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The diagrams whose interference yields the leading contribution to CP- violation 
from Nji-^^ decays in the traditional, Majorana version of leptogenesis. In a su- 
persymmetric theory, the supersymmetrized versions of these diagrams must 

also be included 

A diagrammatical representation of neutrinoless double-/? decay 

A schematic representation, after [361 EH] , of the evolution of baryon number 
B (vertical axis) and lepton number L (horizontal axis) in Dirac and Majorana 
leptogenesis. In Dirac leptogenesis (left panel), the evolution of B and Ltot 
proceeds in three steps: first, two stores of lepton number Li (stored in left- 
handed neutrinos) and Lj,^ (stored in right-handed neutrinos) are produced 
during heavy particle decays; second, sphaleron processes (which act along 
lines of constant B — L) mix Li and B while leaving L^^ alone; third, after 
sphaleron interactions have effectively shut off, equilibration between Li and 
Lj,^ results in a net positive B and L for the universe. This is qualitatively 
quite different from the situation in Majorana leptogenesis (right-hand panel). 
Here, lepton-number-violating heavy particle decays first build up a single, 
nonzero store or lepton number Li = Ltot', second, sphalerons transmute this 
lepton number into a nonzero baryon number for the universe. Only one store 
of lepton number is created, and the result is a universe with negative and 

positive B 

Diagrams that give the leading contribution to the CP asymmetry from decays 
of the scalar fields 0i and (p^. Similar CP asymmetries are generated during 
the decay of the fermionic fields tp^^ and tp- 
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2.3 Diagrams corresponding to the leading order contribution to the equihbra- 
tion of left-handed leptons with right-handed neutrinos. These include, from 
left to right, Higgs decays and inverse decays, scatterings off Standard Model 
fermions and electroweak gauge bosons through a virtual Higgs, and scattering 
off a Higgs boson in conjunction with the emission of an electroweak gauge 
boson. The right two diagrams also have t-channel equivalents. The amplitude 
associated with each diagram is proportional to the (small) effective neutrino 
Yukawa given in equation (12.91) 

3.1 Contours of VLisph"^ |t38j corresponding to the upper and lower bounds from 
WMAP (13. 4p . as well as a simple overclosure bound, as a function of the LSP 
mass irtisp and the reheating temperature Tr associated with cosmic inflation 
for a Wino LSP. The requirement that VLlsp not overdose the universe severely 
constrains Tr, and hence the temperature scale of thermal leptogenesis, for a 
theory like PeV-scale loop-split supersymmetry, in which the LSP is particu- 
larly heavy. For a Higgsino LSP, the results are similar, but the contours move 
slightly to the left 

3.2 Here, we show the variation of (the improved expression for the LSP 
relic density that accounts for the effect of LSP annihilations) with LSP mass 
for several different values of ^3/2 [38]. The horizontal line corresponds to the 
WMAP upper bound on VLcdm- For a given choice of 7/13/2, portions of the 
contour that fall below the horizontal line respect the WMAP constraint and 
are phenomenologically allowed, and the reheating temperature Tr be increased 
beyond the naive upper bound from figure 13. 1[ Portions that lie above it 
overproduce the LSP and are excluded (for example, when 7713/2 = 10^ GeV, 
an LSP mass misp ^ 500 GeV is excluded) 

3.3 A schematic representing the bounds on the leptogenesis scale from grav- 
itino physics as a function of the gravitino mass 7^73/2. When 7^73/2 < 10^ GeV, 
constraints from BBN are quite severe. In the range 10^ GeV < 7r73/2 ^ 
10^ GeV, M is still limited by nonthermal decays to the LSP. For 7/73/2 < 
10^° GeV, decays to gluinos are a concern. This leaves a window 10^ GeV < 
777.3/2 < 10^° GeV within which, depending on the mass of the LSP, Dirac 
leptogenesis can be successful 
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Here, the regions of — 63 space (see equation fl3.30p for a description of 
the Yukawa-matrix parametrization used) for which all constraints on neutrino 
masses and mixings fl3.19p are simultaneously satisfied for some combination of 
the remaining Oj and hi are shown for two different values of 6 [38] . Additionally, 
contours depicting the ranges for S13 = sin 6*13 (which depends primarily on 03, 
but varies slightly with the remaining Oj and hi) are shown. In the right panel 
5 = 10~^; in the left panel, 6 = nie/m^ = 4.83 x 10^^ The plots reveal that 
while 63 is constrained, along with most of the other and hi, to lie reasonably 
near 1, a™"^(5) can be quite large and increases with decreasing 6 = 10~^, 
making it easier to obtain a realistic baryon-to-photon ratio rj. For 6 = 10~^, 
^max ^ ^ _ ^ gg ^ IQ^^, a™"^ ~ 95. In each panel, the configuration 

that yields the greatest decay asymmetry e is marked with an arrow 

Feynman diagrams giving the two leading order contributions to the flavor- 
changing process £~ £"7 due to sneutrino (left diagram) and charged slepton 

(right diagram) mass mixings 

Feynman diagrams for the leading-order process involving (left diagram), 

and for the leading process involving 5m\^ (right diagram), with sneutrinos 
running in the loop in the mass-insertion approximation. Note that any process 
involving 5m\^ necessarily involves two mass insertions, and any one involving 
5m\^ (given by the diagram on the right with an additional Sm'j^j^ insertion) 

necessarily involves three 

Diagrams for 2 2 s-channel processes which transfer lepton number between 
LaggVis (the aggregate lepton number in the sector comprising the fields i, i, 
ur, cr, and cr, which are assumed to be in chemical equilibrium with one 
another due to rapid gauge and (F^)-term equilibration interactions); and Lhid 
(the lepton number stored in right handed neutrinos) . The two t-channel inter- 
actions associated with each diagram appearing above must also be included 

in calculating the full thermally averaged cross-section 

This figure, taken from |43j, illustrates the effect of second-order processes of 
the form 0i0i —>■ ij and (picpi — »• ij on the exclusion contours from leptogenesis. 
Contours are displayed for the case without annihilation (dashed line) and with 
annihilation (solid line) 
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Bands in / — 5 parameter space for which the final baryon number Bf falls 
within the range permitted by WMAP, for different choices of Here 
/ parameterizes the couplings of the $j fields (see eq. f l3.30p ). The configuration 
of the left (right) panel is the one marked with a dot in the left (right) panel 
of figure 13.41 in which 03 = 95 (03 = 4.5). The shaded vertical lines show 
the constraint on the value of 6 coming from neutrino mixings and masses 
{6 = rrie/m^ ^ 4.83 x 10"^ in the left and 6 = 10~^ in the right). Note that 
there are two points that yield a realistic value of Bf for a given 6 (the points 
at which the grey vertical line intersects a given ribbon). When / is small 
enough, the initial baryon number generated is just enough to be consistent 
with WMAP, while the washout effect of inverse decays and 2^2 processes 
is negligible. In this situation the baryon number generated is independent of 
M$ and its the final value is proportional to the CP violating parameter e (see 
eq. (12.161) ) . The dark dotted curve shows the band of consistent baryon number 
calculated in this "drift and decay" limit. Each "X" marks the point in which 
Td/H = 1 for each different M$j. At these points, the "strong washout" regime 
starts. The now active washout processes reduce an initial surfeit of baryon 
number (due to a larger /) down to an acceptable level (see figure |4]4j). . . . 
These two plots, originally appearing in [38j|, show the evolution of baryon 
number B for = 10^^ GeV and 6 = me/m^ = 4.83 x 10~^, in the two dif- 
ferent regimes that produce a realistic value for the final baryon number of the 
universe. Bp. For the rescaled coupling strength / = VA23/i23 = 1-5 x 10"'^, as 
shown in the left panel, the effects of 2 2 lepton-number-changing processes 
are negligible and the final baryon-to-entropy ratio is the same as that initially 
produced by and decays. For stronger couphng / = 3.8 x 10~^, as shown in 
the right panel, baryon number is initially overproduced, but 2^2 processes, 
which are stronger for stronger coupling, reduce B to an acceptable level by 
the time they freeze out. The left and right panels correspond respectively to 
the lower and upper parts of the 'ribbon' in figure 14.31 
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Exclusion plots |13] combining constraints from both leptogenesis and flavor 
violation in the process fi — > 67. The left-hand panel shows exclusion contours 
in M$^-(x) space for a universal scalar soft mass rris = 200 GeV, with tan/? = 
10; the right hand panel shows the variation of the branching ratio BR{fi — >■ 67) 
with respect to using tan /3 = 3, 10 and 30. In both plots, we have taken 
Ml = 160 GeV, M2 = 220 GeV, and /i = 260 GeV. We have also assigned the x 
superfield an F-term VEV \/{F\^ = 10^ GeV. Such a large VEV results in large 
trilinear couplings between Higgs fields and sneutrinos and therefore induces 
potentially sizeable mixings between left-handed and right-handed sneutrinos 
after electroweak symmetry breaking. In each plot, the thick solid contours 
represent the current experimental bound on the branching fraction (13.331) : the 
dashed lines represent the expected future experimental bound of 1.5 x 10^^^ 
from MEG. The thin solid contour in the left-hand panel delimits the region 

allowed by leptogenesis constraints 

Here, we illustrate the two representative cases for sneutrino dark matter in 
Dirac neutrinogenesis. In each case, we have set M' = 200 GeV and m/>^ = 
150 GeV and maintained the relations g{x) = 200 GeV, g{^i) = 100 GeV, 
/i = 500 GeV, and (i/g = (2/9 = 1- The left panel depicts the situation where 
the U{l)]\r gauge coupling is reasonably strong {g = 9.5 x 10~2). As the universe 
evolves, T decreases and z —>■ 00. Here, the number-density-to-entropy-density 
ratio Ydm of a right-handed sneutrino LSP reproduces the WMAP prediction 
(the horizontal dashed line), but the lepton number Lj>^ stored in right-handed 
sneutrinos is washed out almost entirely and all connection between the dark 
matter abundance and the observed baryon number of the universe is lost. The 
right panel shows the situation where the gauge coupling is small g = 5x 10~'^. 
Here, washout is insignificant and Lc^ persists essentially unchanged from its 
initial value until present time, but sneutrino annihilations are inefficient in 
reducing the dark matter abundance to an acceptable level. The equilibrium 
right-handed sneutrino abundance Y^\j- is also shown in each panel for refer- 
ence. It should be emphasized that in the left panel, the correct dark matter 
relic density is obtained, and that while ^Icdm and Qb are not tied to one an- 
other, there is no problem for cosmology, whereas the situation depicted in the 
right panel is cosmologically unacceptable, since dark matter is overproduced. 
The additional contribution to (j(z>i>*) that arises in the right-handed sneu- 
trino dark matter scenario in the presence of additional light scalars Ui and uJj 
charged under the U{1)n- gauge group 
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5.3 The effect on tlie dark matter number-density-to-entropy-density ratio Ida/ of 
adding a number A^^ of vector-like pairs of light scalars to the theory. Here, 
all such scalars have U{1)n have charge magnitude \Q\ = 20, and results are 
shown for N„, = 10, 10^, and 10^, as well as for the unmodified case N„, = 0. 



As in the right panel of figure 15.11 M' = 200 GeV and mc,^ = 150 GeV and 
maintained the relations g{x) = 200 GeV, g{^i) = 100 GeV, /t = 500 GeV, 
(i/g = (2/g = 1, and g = 1.5 x 10^^. The addition of light scalars does not 
affect the lepton number Lj^^ stored in right-handed sneutrinos appreciably for 
> 10. Because g is small, an inordinately large number of light fields are 
required in order to achieve a situation in which L^^^ dictates the size of Ida/- 
When iV^ is large enough, L^,^ creates an effective "floor" for the dark matter 
abundance (corresponding to a dark matter abundance consisting only of z/"r 
and none of its antiparticle zz/j*) |96 

6.1 (</('?/')) AO (the line-of-sight integral through the halo density squared, averaged 
over the angular acceptance AQ) vs. AQ for several halo profiles: smooth 
isothermal, clumpy isothermal, Burkert, NFW, and Moore et al. Taken from [6T] . ll01 

6.2 The expected flux from the annihilation of a heavy Wino LSP, as a function 
of Afl, that would be detected by a satellite detector aimed at the galactic 
center [61]. Also included is the anticipated background flux for a 2.3 TeV Wino 
at such a detector for different values of detector energy resolution, ranging from 
AE/E = 50% (top line) to AE/E = 5% (bottom line). It should be noted 
that the spread in the signal is smaller than the width of the energy bin. From 
this, it is apparent that for the NFW and Moore et al. profiles, the prospects 
for detection increase with better angular resolution (decreasing AQ). For a 
Higgsino LSP, the resulting curves are similar, but the signal is two orders of 
magnitude lower. The background flux for a 1.1 TeV Higgsino also is increased 
by a factor of ~ 10 over the 2.3 TeV Wino case, owing to the lower energy of 

the signal photons Il05 



6.3 The ratio of Nsigaai/ a/ A^background (left panels) and total number of photons 
(right panels) collected by a generic space telescope with an effective area of 
2 X 10^ cm^ and an energy resolution of 1%, over a range of AQ, after 10^ s 
(about 1/3 of an active year) of viewing time, and for both Wino (top panels) 
and Higgsino LSP (bottom panels) [61]. See figure 16.41 caption for the halo 
model key. The threshold for 5a discovery has been included for reference in 
the significance graphs, and contours corresponding to 10, 25, and 100 events 
have been included in the event count graphs. It can be seen here that for such a 
space telescope, no halo profile is capable of producing the 25 events necessary 
for detection, and that only for the Moore et al. profile is the significance 
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xi 



6.4 The expected flux from tlie anniliilation of a heavy Wino LSP, as a function of 
AQ, that would be detected by a Cherenkhov detector aimed at the galactic 
center [^1]. Also included is the anticipated background flux from a 2.3 TeV 
Wino at such a detector for different values of detector energy resolution, rang- 
ing from AE/E = 50% (top line) to AE/E = 5% (bottom line). It should 
be noted that the spread in the signal is smaller than the width of the energy 
bin. From this, it is apparent that for the NFW and Moore et al. profiles, 
the prospects for detection increase with better angular resolution (decreasing 
Afl). For a 1.1 TeV Higgsino LSP, the resulting curves are similar, but the 
signal is two orders of magnitude lower. The background flux is also increased 
by a factor of ~ 10 over the 2.3 TeV Wino case, owing to the lower energy of 
the signal photons 

6.5 The ratio of A'signai/ \/ A'background (left panels) and total number of photons 
(right panels) collected by a generic Cherenkhov array with an effective area 
of 1.5 X 10^ cm^ and an energy resolution of 10%, over a range of AQ, after 
10^ s (about 1/3 of an active year) of viewing time, and for both a heavy Wino 
(top panels) and Higgsino LSP (bottom panels) [61]. See figure [6^ caption 
for the halo model key. The threshold for 5a discovery has been included for 
reference in the significance graphs, and contours corresponding to 10, 25, and 
100 events have been included in the event count graphs. It can be seen here 
that there are real prospects for detection with such a Cherenkhov detector, 
provided that the galactic CDM halo density resembles the NFW or Moore et 
al. profiles 

6.6 Detection boundary contours in Aea — AQ parameter space for the 7-ray sig- 
nature of a 2.3 TeV Wino, based on the 5a significance requirement from equa- 
tion EH for VERITAS (top left panel), with = 1 x 10^ cm^ and AE/E, 
and GLAST (bottom left panel), with Aes = 1.5 x 10^ cm^; and based on the 
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scenario 

D.2 Processes contributing to a{i)i)*) in the right-handed sneutrino dark matter 
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ABSTRACT 

Theory and Phenomenology of Dirac Leptogenesis 

by 

Brooks D. Thomas 

Chair: James D. Wells 

Dirac leptogenesis, in which neutrinos are purely Dirac and develop small but nonzero effective 
masses without the aid of the see-saw mechanism, provides an interesting alternative to the 
standard leptogenesis picture. Here we review the theory and phenomenology of Dirac leptoge- 
nesis and show that it is a viable theory capable of simultaneously satisfying all relevant bounds 
from cosmology, neutrino physics, and flavor violation. In addition, we also explore several po- 
tential extensions of the model, such as the possibility of right-handed sneutrino dark matter 
and the potential for relating the leptogenesis mechanism to the origin of the /x-term. Theories 
with a heavy gravitino and gaugino masses generated by anomaly mediation emerge as one 
natural context for Dirac leptogenesis. In such models the lightest neutralino is often expected 
to be predominately wino or Higgsino, and is a viable dark matter candidate. We conclude 
with an examination of the prospects for detecting the effectively monoenergetic photon signal 
that results from the annihilation of such a dark matter particle in the galactic halo. 



CHAPTER I 
BARYOGENESIS AND LEPTOGENESIS 



1.1 The Problem of the Baryon Asymmetry 

The universe we live in is manifestly asymmetric between baryons and antibaryons. This 
statement can be quantified by introducing a parameter rj, defined as r] = UB/n^- Here 
Ub = Uh — n-f^, where rib and are the baryon density and antibaryon density of our universe, 
respectively, and is the present number density of photons. The value of t] has recently been 
measured with great precision by WMAP [Ij to be within the range 



which implies that the relic density of baryonic matter in our universe is somewhere in the 
range 



Since the standard cosmology is manifestly symmetric with respect to matter and antimatter, 
it must be supplemented with some additional physics which can account for the presence of a 
baryon asymmetry on the observed scale. The generic name for such scenarios is baryogenesis. 

1 



rj = (6.1 ±0.3) X 10 



10 




0.021 < Qbh^ < 0.025. 
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Any successful baryogenesis model must satisfy a set of three requirements which were 
originally pointed out by Sakharov [2J. The first and most self-evident of these conditions is 
that there must be baryon number violation; otherwise, a universe with an initial value of = 
could not evolve to one where rj ^ The second is that the charge conjugation symmetry C 
and its composition with parity CP must also be violated. If this were not the case, baryon 
number violation and antibaryon number violation would occur at the same rate, and hence 
no net baryon asymmetry would be created. Finally, one must have either a departure from 
thermal equilibrium or else posit that CPT, the composition of CP with the time evolution 
operator, is somehow violated during the early universe. This condition can be obtained by 
calculating the thermal average (B) of the baryon number, a dimensionless quantity defined by 
the B = {rih — rij;) /s in terms of the entropy density 

s = %9.T\ (1.3) 

where T is the temperature of the thermal bath and is the number of interacting degrees 
of freedom (in the minimal supersymmetric model (MSSM) during the baryogenesis epoch, 
~ 205). The result is [3]: 

{B) = Ti[e-^^B] = Ti[{CPT){CPT)-^e-^^B]. (1.4) 

If CPT commutes with the Hamiltonian, this becomes 

(B) = Ti[e-^^^B] = Tr[e-'^^{CPT)-^B{CPT)] = -Ti[e-^" B], (1.5) 

which implies B = (conversely, if thermal equilibrium is not maintained, all particle species 
will not have a common temperature, and trace cyclicity can no longer be invoked). If any 
of these three conditions is not met, t] will not deviate from zero. In general, since Lorentz 



3 

invariance is believed to be necessary for the formulation of a consistent quantum field theory, 
this condition has been taken to imply a departure from thermal equilibrium rather than the 
presence of CPT-violation. The alternative, in which Lorentz invariance is temporarily, dynam- 
ically violated in the early universe, was originally explored in the spontaneous baryogenesis 
scenario of [1] and has been elaborated on in models like the radion-induced baryogenesis of [5] 
and the gravitational baryogenesis of [6j], in which a net baryon number is produced via a CP- 
violating effective interaction between the baryon number current and the derivative of the Ricci 
Scalar. Indeed Lorentz invariance is violated whenever a tensor field receives a nonzero vacuum 
expectation value (VEV), and indeed may arise spontaneously in string theories, brane- world 
scenarios, and modified gravity models [7]. Still, no experimental signal of CPT violation has 
ever been detected [HI E] , and most theories of baryogenesis (including the one with which this 
work is principally concerned) assume that CPT is a good symmetry of the Hamiltonian and 
that the generation of a baryon number for the universe involves out-of-equilibrium dynamics 
rather than Lorentz violation. 

1.2 Baryogenesis via Leptogenesis 

1.2.1 Leptogenesis and the See-Saw Mechanism 

A variety of viable baryogenesis models exist, including electroweak baryogenesis, in which 
CP- violation occurs at a bubble wall, or phase boundary, and Affleck-Dine baryogenesis [TOj . 
in which the baryon asymmetry is generated by moduli fields charged under B — L, where 
L = {ne — rij) / s denotes the lepton number of the universe. In this paper, we will focus on 
models which achieve baryogenesis through a framework known as leptogenesis [HI [12] , where 
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decays of heavy particles in the early universe which violate both CP and lepton number L 
produce an initial lepton asymmetry, which is then converted to a nonzero baryon asymmetry 
by sphaleron processes associated with the SU{2) electroweak anomaly [13j. Leptogenesis is 
a particularly attractive model because in addition to its ability to yield a realistic value for 
T] [E], it can also explain why the standard model neutrinos have small but nonzero masses. In 
its most common form, which we will call Majorana leptogenesis, the role of the heavy, decaying 
particles is played by the right-handed neutrinos uji required to fill out the multiplets housing 
the Standard Model fields in many grand unified groups. As they are gauge singlets, nothing 
prevents them from obtaining large (GUT-scale) Majorana masses M,^^, and such masses are 
inherently lepton-number- violating. 

In the following brief review of Majorana leptogenesis, and for the remainder of this work, 
we will work under the assumption that the universe we live in is inherently supersymmetric, 
but that supersymmetry is broken at some high scale Msusy- Although leptogenesis certainly 
does not require supersymmetry, we choose to work in a supersymmetric framework for a variety 
of reasons. Supersymmetry naturally explains the stability of the weak scale against quadratic 
divergences [151 [IS] and improves the outlook for the unification of the Standard Model gauge 
couplings pn [T8]. Furthermore, the assumption of conserved (or nearly conserved) i?-parity 
makes the lightest supersymmetric particle (LSP) a convenient dark matter candidate |19j . 

In supersymmetric Majorana leptogenesis, the most general realizable superpotential one 
can write for the relevant fields is 



where and Na represent the left- and right-handed neutrino superfields respectively, Hu is 



W 3 VafsLaNpHu + M.^^NM 




linglll 



the up-type Higgs superfield, and yap is a dimensionless trilinear couplingl^ This superpotential 
contains both Majorana and Dirac masses for neutrinos (once the Higgs field gets a VEV) and 
yields a mass matrix 



( n ^ 



( \ 



:i.7) 



where fiavor indices have been suppressed to simplify notation. Since mo = yvsin^, where 
tan/3 = Vu/vd denotes the ratio of the up- and down- type Higgs VEVs, the off-diagonal terms 
mixing ul and z//j will be small compared to M^^^a- As a result, when this matrix is diagonalized, 
the mass spectrum of the theory contains three light neutrinos with masses 



If / N 

m^ = mDY7—mD, (1-8) 



which are identified with the Standard Model neutrinos, as well as three heavy neutrinos with 
masses 0{M,y^). This method for obtaining small but nonzero neutrino masses is known as the 
see-saw mechanism [201 EI] , and the liaison it forges between the origin of the observed baryon 
asymmetry and the lightness of the Standard Model neutrinos is one of the most compelling 
aspects of leptogenesis scenarios. 

One of the useful consequences of Majorana component to the light neutrino masses is 
that such a source of lepton number violation would be observable in neutrinoless double-/? 
decay experiments ^22j. The contribution to this process from a Majorana neutrino mass is 
shown diagrammatically in figure 11.21 It should be noted that a Majorana neutrino mass is 
not required for it to occur (certain other other beyond-the-Standard-Model sources of lepton 
number violation can also contribute), but if neutrinoless double-/? decay were ever observed, 
it would serve as compelling evidence in favor of the see-saw mechanism. 



"^We have chosen to work in a basis where the M^^^a are diagonal. 
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Figure 1.1: The diagrams whose interference yields the leading contribution to CP- violation 
from Nfi^ decays in the traditional, Majorana version of leptogenesis. In a super- 
sjTiimetric theory, the supersymmetrized versions of these diagrams must also be 
included. 




Figure 1.2: A diagrammatical representation of neutrinoless double-/? decay. 

The production of lepton number in Majorana leptogenesis occurs when the heavy right- 
handed neutrinos decay out of equilibrium. The leading CP- violating contribution comes from 
the interference between the diagrams shown in figure 11.11 (in supersymmetric models, the su- 
persymmetrized versions of these diagrams must also be included), and results in the generation 
of a nonzero lepton number L for the universe in the traditional, Majorana version of leptoge- 
nesis. In the limit where there is a large splitting between the right-handed sneutrino masses 
(Mj^^i <^ M^^2, M^rs) The contribution from a single Nji^ decay is commonly called the decay 



asymmetry and is given by [23] 

3 lv[i[yl^yipylpyia\ ( M^^i 



where a sum over repeated indices is assumed. This translates into a net lepton number pT 



1.9) 



L = Kw— (1.10) 

g* 

for the universe, where is a dimensionless, model-dependent prefactor included to account 
for the washout of lepton number by 2 ^ 2 lepton-number- violating processes, etc. that deplete 
the lepton number genrated during Nr^ decay. Whether Majorana leptogenesis is capable of 
reproducing the observed value of rj depends both on e and on kw- 

1.2.2 Effects of the Electroweak Anomaly 



So far we have seen that the leptogenesis mechanism is capable of generating a net lepton 
number L 7^ for the universe. A mechanism by which this lepton number can be transformed 
into a net baryon number is built into the Standard Model (SM) [24] . The baryon and lepton 
number currents 

i 

= + (1.12) 

i 

where the sum is over generations, are anomalous due to the triangle anomaly, with divergence 

d^J^ = d^Ji = ^ [glW^^^W^^-^ + glF.^F^'^) . (1.13) 

Here, W^^, and F^^ are the respective field-strength tensors of the SU{2) and U{1)y gauge 
fields, g2 and qy are the respective gauge coupling constants, Nf denotes the number of fermion 



generations, and 



= ^e^'^^'Fp, (1.14) 



is the dual of F^^ (the expression for W""^^ is analogous). The change in baryon number over 
some duration t is then given in terms of the Chern-Simons number 

Ncs{t) = j|l / d'xe'^'ea.cW^W^Wm (1-15) 

by 

5 (t) - 5(0) = j Sxd^'J^ 



= Nf[Ncs{t)-Ncsm- (1-16) 

An analogous equation exists for lepton number. 

This theory has an infinite number of quasi-degenerate vacua, in each of which the W"" are 
pure gauge and consequently Ncs becomes an integer. Equation (11.161) tells us that vacuum- 
to- vacuum transitions involve AB = AL = NfANcs, so the minimum change for such a 
transition in the SM, where Nf = 3, is AB = AL = 3. This corresponds to an effective 
12-fermion interactioEC 

i 

where we have included the B + L subscript to draw attention to another important consequence 
of equation (I1.16p : since the same equation describes the evolution of both baryon and lepton 
number, the combination B — L is conserved in any process of the form (11.171) . while B + L 
is violated by at least 6 units. These effective interactions serve as the primary means of 
conversion between B and L in leptogenesis. 



^The physical scale A associated with this effective operator is the inverse of the magnetic screening length 
Rsph ^ l/{a2T). 
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In calculating the transition rate for these interactions, we will make use of the fact that 
the potential energy of the electroweak theory between any pair of adjacent {ANcs = 1) vacua 
contains a saddle point, and that transitions through that point in field space should dominate 
over all others. This field configuration is known as the electroweak sphaleron, and the height 
of the potential barrier separating any two adjacent vacua is known as the sphaleron energy, 
usually denoted by Esph- At zero temperature, the rate for such transitions is determined by 
the instanton action and turns out to be completely negligible |2i 



At finite temperatures, however, the situation is modified by the possibility of thermal exci- 
tations over the barrier. For low temperatures T < Egph, the rate calculation is reasonably 
straightforward and the result [25] is 



sph 



2^6 



-Esph/T 



(1.19) 



where c is an 0(1) constant and Mw{T) is the mass of the W boson as a function of temperature. 

At higher temperatures T > Egph-, the calculation of the sphaleron interaction rate is some- 
what difficult, but it can be estimated on dimensional grounds by examining the scales of 
the processes involved [26] . Non-perturbative fluctuations of the gauge field that mediate the 
sphaleron transition are associated with magnetic fluctuations with characteristic distance scale 
Rsph ~ 1/ {012T) on the order of the magnetic screening length. If one naively assumes that the 
time scale for the process is of the same order, one obtains the result T sph ~ a2^. However, 
when one takes into account damping effects in the plasma [27J [28] one finds that the time scale 
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is slowed to tsph ~ l/lc^i^); which leads to a rate 

r^p/, = (25.4 ± 2.0)a^T, (1.20) 

where a2 = (7|/47r 1/30 and the value given for the proportionality constant is the result of 
numerical calculations [30l [31] . 

Sphaleron interactions will be in equilibrium whenever the sphaleron transition rate Tsph 
exceeds the expansion rate of the universe, expressed by the Hubble parameter H. When the 
universe is radiation-dominated, the expansion rate is given by 

H = IMgl^^T'^/Mp. (1.21) 

From equations fll.l9p and (11.201) we find that sphaleron interactions are in equilibrium when 

100 GeV < T < 10^^ GeV, (1.22) 

which means that the interconversion of B and L can be considered rapid perhaps during and 
certainly soon after the lepton number asymmetry L is built up by the decays of the heavy 
right-handed neutrino fields and will remain rapid down to around the weak scale. 

In order to relate these baryon and lepton number asymmetries to each other in a quan- 
titative way, we can take advantage of the set of conditions among the chemical potentials fii 
implied by equilibrium conditions among the species i present in the thermal bath. If chemical 
equilibrium is established between particle species ai, 02 . . . , cim, &2, • • • &n by sufficiently rapid 
interactions of the form oi, a2 . . . a/ ^ 6162 • • • their chemical potentials obey the relation 



$^/^a, = ^/ib.. (1.23) 



i=l i=l 



■^It has been argued that the characteristic time scale of these processes may also have a logarithmic 
suppression, so that Tsph ^ [0^2'^ ^'^(l/'^)]- Since what is of interest to us is the numerical result for the rate of 
baryon-number-changing transitions, we can remain agnostic on this issue. 
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For example, when sphaleron interactions are in equilibrium, equation ( 11.17P implies that 

^(3/.,,+/i,J = 0. (1.24) 

i 

Similarly, 5'f/(3) QCD instanton processes and global hypercharge conservation respectively 
require that 

5^(2/i,,-//„,-/i,J = (1.25) 

i 

5^ i^^q^ + 2/i«, - ^^d, - - /ie, + ■^^^Hj = o, (i.26) 

and from the fermion Yukawa interactions 

5^(2/x,,-/i^-/i,J = (1.27) 

i 

J2i2fig^ + fiH - fin) = (1.28) 

i 

5^(2/i,,^-/XH-/^eJ=0. (1.29) 

i 

Since particles of different generations will also be in equilibrium with one another at high 

temperatures, we can take fig^ = fig^, fi^, = fJ'u, f^d, = f^d, fJ-e, = IJ-e, fJ-e, = fJ-e and solve this 

system of equations for one of the /i, (we choose fit). The resulting chemical potentials are 

_2Nj + 3 eN,+i _2Nf-l 

^^-6iv7T3^^ /^^--sivFa/^^ ^"-6iv7T3^^ 

/^g = -^f^i /^H = WJ+3^^^^ {1.30) 

and since the values of B and L are related to these chemical potentials (see equation (IB.SP in 
appendix [B|) by 

1 v-^ Nf . , , , 

B = — - gifii = — ^(2^g + /i« + /id) (1.31) 

i=baryon 

L = Yl = -^{2fie + fie), (1-32) 

5'*s-' , 9*3-'- 

i=baryon 
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one finds that 

+ 4 _ 14iV^ + ^ _ (1.33) 

22A^/ + 13^ ^ 22Ary + 13^ ^ ^ ^ 

in the Standard Model. From this, we can state the relationship between baryon and lepton 
number when sphaleron interactions are in equilibrium: 

B = —L. (1.34) 
51 ^ ' 

In the MSSM, things are modified by the presence of a second Higgs doublet, and the result 
becomes result is 

B = |l. (1.35) 

In Majorana leptogenesis, we know that since B — L is conserved by sphalerons and violated 
only by the lepton-number-producing decays of Nfi^^, the initial value {B — L)init = —Linu 
generated during the leptogenesis epoch and the present value will be equal; hence (in the 
MSSM) 

Btoday = ——Linit = — — I^W — (1.36) 

and the universe receives a nonzero baryon number. 
1.2.3 An Alternative Leptogenesis Scenario 



Although it is not our aim to discuss the details of Majorana leptogenesis models, we note 
that they have been shown to be able to yield a realistic value of 77, reproduce the observed light 
neutrino spectrum, and evade problems associated with relevant astrophysical constraints (for 
recent reviews, see [HI [231 E21 ES]). In this form, leptogenesis emerges as a viable phenomeno- 
logical theory. There may be other forms of leptogenesis that are just as successful, however. 
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The aim of this work is to conduct a thorough examination of one such alternative, supersym- 
metric Dirac leptogenesis, and to show that it is also a phenomenologically viable model. We 
will begin by discussing the symmetries and field content which determine the form of super- 
potential and discuss the consequences of that superpotential for baryogenesis and neutrino 
masses in chapter [III In chapter IIIIl we enumerate the phenomenological constraints on Dirac 
leptogenesis and construct a simple, theoretically- motivated model. In chapter IIV| we solve the 
Boltzmann equations for the evolution of baryon and lepton number in the early universe nu- 
merically and discuss the consequences of these results on the model parameters. We investigate 
some potentially promising extensions of the model in chapter |Vl and in chapter IVIIi 



CHAPTER II 



DIRAC LEPTOGENESIS 



2.1 Super potential and Fields 



While Majorana leptogenesis is a certainly successful model, it is not the only way in which 
leptogenesis can be realized. There are several reasons why exploring potential alternatives 
is a worthwhile endeavor, and in particular why it is advantageous to have a viable leptoge- 
nesis mechanism in a model without Majorana neutrinos. One is that the non-observation 
of neutrinoless double-/? decay at future experiments [221 El] could significantly constrain the 
parameter space of Majorana leptogenesis to the point where severe model tensions might arise 
between such constraints and others from astrophysics, flavor-physics, etc. Another is that 
it is not at all obvious that massive right-handed neutrinos emerge naturally from a string- 
theory context. While landscape surveys, even for particular classes of models and specified 
orbifold compactifications, are limited by computational complexity, preliminary results [55] 
performed for heterotic BSL^ models on the Z3 orbifold suggest that a see-saw mechanism in 
the traditional sense is not a generic feature of otherwise phenomenologically promising string 
models. 
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For all these reasons, it would be useful if one could find a way to link the smallness of the 
physically-observed neutrino masses to a successful baryogenesis mechanism without having to 
introduce singlet neutrinos with Majorana masses. It is indeed possible to do this in the context 
of a scenario that has come to be known as Dirac leptogenesis or Dirac neutrinogenesis [36], [37] . 
In this scenario, an additional symmetry (the precise form of which is not terribly important or 
stringently constrained by the model framework) is introduced, and charges are assigned under 
this new symmetry in a manner which forbids, at tree level, both the Majorana and Dirac 
neutrino mass terms appearing in (11. 6p . A set of heavy, vector-like pairs of fields introduced 
whose couplings to the standard model fields contain nontrivial, CP-violating phases are also 
introduced. These fields will play the role that heavy right-handed neutrinos play in Majorana 
leptogenesis, and their decays during the early universe will lead to the buildup of equal and 
opposite lepton asymmetries and Lj,^ in the left-handed lepton and right-handed neutrino 
sectors, while conserving the overall lepton number for the universe Ltot = Li + L^^ = 0. The 
equilibration rate between these stores is suppressed by the smallness of the effective neutrino 
Dirac mass term; as a result, the electroweak sphaleron processes which convert Li into a baryon 
asymmetry B effectively shut off before Li and L,^^ have a chance to equilibrate. Consequently, 
unlike in Majorana leptogenesis, the universe ends up with a net positive lepton number as well 
as a net positive baryon number — a result which is depicted schematically in figure 12. 1[ 

We will begin by writing down the superpotential for Dirac leptogenesisjll which is modified 
from that of the MSSM only in the lepton sector. The field content and charge assignments of 
the model are essentially the same as those presented in [37\. Of the usual quark and lepton 

supermultiplets of the MSSM, the left-handed lepton multiplets La (a is a family index) and 

^Once again, we remark that Dirac leptogenesis does not require supersymmetry and functions perfectly well 
without it [361 • 
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Figure 2.1: A schematic representation, after [36l [38] , of the evolution of baryon number B 
(vertical axis) and lepton number L (horizontal axis) in Dirac and Majorana lepto- 
genesis. In Dirac leptogenesis (left panel), the evolution of B and L^t proceeds in 
three steps: first, two stores of lepton number Li (stored in left-handed neutrinos) 
and Lj,^ (stored in right-handed neutrinos) are produced during heavy particle de- 
cays; second, sphaleron processes (which act along lines of constant B ~ L) mix Le 
and B while leaving L,y^ alone; third, after sphaleron interactions have effectively 
shut off, equilibration between Li and L,^^ results in a net positive B and L for the 
universe. This is qualitatively quite different from the situation in Majorana lep- 
togenesis (right-hand panel). Here, lepton-number- violating heavy particle decays 
first build up a single, nonzero store or lepton number Le = Ltot', second, sphalerons 
transmute this lepton number into a nonzero baryon number for the universe. Only 
one store of lepton number is created, and the result is a universe with negative 
Ltot and positive B. 



the Higgs multiplets Hu and Hd will be pertinent to leptogenesis. We also include a right- 
handed neutrino superfield for each family, an exotic chiral multiplet aiid a number N<p 
of vector-like pairs of chiral multiplets and $j. While at least two such pairs are required 
for leptogenesis, any iV$ > 2 is in principle allowed from a baryogenesis standpoint. We also 
need to introduce an additional symmetry to forbid Majorana masses for the Na, which we 
choose to be an additional global U{1) we will call U{1)n, due to the fact that the right-handed 
neutrinos are charged under it. The charge configurations of these fields under this additional 



f/(l), as well as the rest of the relevant symmetries of the theory, are shown in table I^TTI 



The 



^With the field content given in table \2A] the global U{1)n is anomalous. To deal with this issue, one 
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Field 


U{1)l 


U{1)n 


SUi2) 




Pm 


N 


-1 


+ 1 


1 





-1 


L 


+1 





2 


1 

2 


-1 


Hu 








2 


1 
2 


+1 


Hd 








2 


1 
2 


+ 1 





+1 


-1 


2 


1 
2 


-1 





-1 


+1 


2 


1 
2 


-1 


X 





-1 


1 





+1 



Table 2.1: One possible set of charge assignments, taken from [37], that leads to the Dirac 
leptogenesis superpotential given in (12.11) . Here the additional symmetry employed 
is a U{1) (which may in principle be either global or local). Only the charges of 
the fields relevant to leptogenesis, which include the Higgs doublets Hu and Hd, 
the left-handed lepton superfield L, the right-handed neutrino superfield A^, the 
heavy fields $ and $, and the additional field Xi have been included. Here, U{1)l, 
SU{2), and U{1)y respectively denote lepton number, SU{2), and U{1) hypercharge 
quantum numbers, Pm denotes matter parity, and U{1)n denotes the charge under 
the additional U{1). 

most general superpotential that can be constructed out of this set of fields is 

W 3 Xic^N^^.Hu + hiaLa^iX + Mij^i'^j + ^HuHd, (2.1) 

where a is a family index, Aj^ and hia are Yukawa couplings, Mij is a matrix of supersymmetry- 
respecting mass terms coupling the $j and $j fields, and fi is the usual Higgs mass parameter. 
We note that we can always choose to work in a basis where the the mass are diagonal and 
real, and hence from this point onward, without loss of generality, we will adopt the convention 

that Mjj$j$j = M$^$j$j and absorb all phases into the coupling matrices Xia and hia. 

can modify the theory at high energies by introducing additional heavy fields, appealing to the Green-Schwarz 
mechanism [39 , etc. Alternatively, one may modify the field content of the low-energy effective theory and 
introduce a set of fields with the appropriate charges to cancel the U{1)n anomalies. We will explore one such 
choice in section 15.31 
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2.2 Dirac Leptogenesis, Baryogenesis, and Neutrino Masses 

In order for successful baryogenesis to occur, we must satisfy the Sakharov criteria discussed 
in section [m The out-of-equihbrium condition is satisfied during the decays of the component 
fields in $ and $ (both the scalar and fermionic components of the $ and $ supermultiplets, 
which we denote by 0, 0, and will play the role that the z//j play in Majorana leptogen- 
esis). The mass of the lightest $ — $ pair, which we will denote M$^, defines the leptogenesis 
scale. Since the parameters Xia, and hia may in general be complex (as may the elements of 
(Mjj) in a general basis), they will in general contain nontrivial CP- violating phases that cannot 
be rotated away. This allows us to satisfy the CP-violation criterion. As in Majorana leptoge- 
nesis, a source of baryon number violation is provided by electroweak sphaleron processes. The 
primary difference between the two models is that explicit lepton-number-violating terms are 
present in equation (II. 6p but absent in equation fl2.ip . Unlike in Majorana leptogenesis, where 
L is explicitly violated by the decays of heavy fields, here the overall values of B and the total 
lepton number L^t are altered only by sphaleron processes, and B — L is never violated. 

When the temperature of the thermal bath drops below the leptogenesis scale , both 
the scalar and fermionic components of the $i and $1 superfields will decay, generating a net 
CP asymmetry and building up stores of lepton number in the lepton fields z/^j and i, and 
in the associated slepton fields upt, i. The leading contribution to CP- violation arises due to 
the interference of tree-level and one-loop-level diagrams. Those relevant to (p and (p decay 
are shown in figure 12. 2j the fermion fields and ip^_^ undergo similar decays, but in the 
approximation of unbroken supersymmetry the amplitudes (and resultant CP-asymmetries) in 



19 




Figure 2.2: Diagrams that give the leading contribution to the CP asymmetry from decays of 
the scalar fields (pi and 0^. Similar CP asymmetries are generated during the decay 
of the fermionic fields ih^ and th- . 

the fermion case will be the same as those for the scalar case, so the rates need not be separately 
evaluated. For completeness, we have included contributions involving the scalar component 
of Hd, which start to become important when the supersymmetric Higgs mass parameter fi is 
of the same order as the M$ . , but in what follows we will assume that /x -C Mi and hence the 
contribution from these diagrams is negligibly small. 

For purposes of illustration, let us begin by examining a simple toy model. We will consider 
the case where there are only two sets of $ and $, the minimum number required for CP- 
violation, in which case [lO], one may parameterize the associated lepton number violation by 
defining a single decay asymmetry e, which represents the amount of lepton number generated 
in any particular lepton-number-carrying species by the decay of a single heavy particle. This 
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implies the relations 



r($i ^ iV^i/^) - r(<l>5 — ^ iV,ff„) = eVo (2.2) 

r(<i>i ^ - r($i ^ L^x') = -er^ (2.3) 

r($i ^ L^X') - r($i ^ = eTD (2.4) 

r($, iv,i^„) - r($i iv^ff^) = -er^ (2.5) 



among the rates for the processes depicted in figure [2T2| where F^i is the total decay width of 
any of the heavy fields in the $i or $i supermultiplets, and we have used the superfield notation 
for $1, N, etc. because in the assumption of unbroken supersymmetry, the super symmetrized 
versions of the diagrams appearing in fig. 12.21 yield the same result as the unsupersymmetrized 
ones. Explicit calculation of Yd and e yields 

= -^M^, V (|Ai„p + . (2.6) 

OL 

and 



47r(|M$,|2 - |M<i,J2)(|Ai^|2 + |/ii^|2)' 
where in both equations, a sum over the repeated indices a, /3 and 7 is implied. It will be 
convenient to define a parameter 5 = |M$J/|M$2|, and in terms of 5, 

Im(ALA2a/it^/^2/3e^^) ( b , g) 



47r(|Ai^|2 + |/ii^|2) Vl-^sy 

where il) is the relative phase between and M^^. This tells us that for small values of 5, the 
final baryon-to-photon ratio will be approximately proportional to 5. Since and have 
equal and opposite charges under the global U{X)l symmetry, the individual Li and Ly^ lepton 
numbers respectively stored in left-handed leptons and right-handed neutrinos will likewise be 
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equal and opposite. Supersymmetry enforces a similar condition = —Lc,^^ in the sneutrino 
sector, and consequently no net lepton number is produced by the decays of $1 and $1. 

In order to discuss the subsequent evolution of these stores of lepton number, we must now 
take a moment to examine how things look at temperatures far below the leptogenesis scale 
Ml. Here, the theory can be described by an effective superpotential Wgjf in which the heavy 
$i and $i have been integrated out: 

^-ff 3 ^^xLpHuNa^ + i^H^Hd. (2.9) 

If we arrange for the scalar component of the x superfield to acquire a VEV (x), an effective 
Yukawa matrix for neutrinos, proportional to the ratio (x)/M$^ will result; then, when 
acquires its VEV during electroweak symmetry breaking, this will translate into a neutrino 
mass matrix with entries given by 

m,^p={x)vsm(3 (2.10) 

If (x) ^ Ml, this setup ostensibly yields small but nonzero masses for neutrinos, even when the 
elements in Ajq, and hi^ are 0{1), and thus stands as an alternative to the traditional see-saw 
mechanism. Additionally, since the mass matrix in fl2.10p has a reasonably simple structure, it 
can yield interesting predictions about the mass hierarchy among the standard model neutrinos, 
especially when certain additional, well-motivated constraints are applied, as we shall see in 
section 13.21 

We have so far said nothing about how x receives its requisite VEV, but there are a va- 
riety of ways of engineering such a thing. One workable example is the O'Raifeartaigh-type 
model employed in [?T], in which the F-term of x acquires a large VEV (F) ~ 'm-^/2Mp, and 
supergravity effects give rise to a nonzero VEV (x) — 1611771^/21^^^ for the scalar component of 
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Xi where k is an undetermined dimensionless coupling constant. Another, in which the VEV 
is induced by introducing new superpotential couphngs involving x with additional exotic su- 
perfields charged under Uil)]^ and no F-term VEV develops for x, is presented in chapter IVl 
For the moment, we will not concern ourselves with the precise manner in which a x VEV 
comes about, but it will be important to distinguish theories based on the presence or absence 
of a large F-term VEV (F^)- '^^^ reason for this is that when (F^) is nonzero, the interaction 
lagrangian resulting from (12.91) includes an effective A-term 

(F ) 

^eff 3 \^\iahipHJpV^ + c.c. (2.11) 

through which left- and right-handed sneutrinos can equilibrate. When (F^) is large, this 
interaction will result in the lepton number asymmetries L^- and Lj>^ stored in the slepton 
sector being rapidly equilibrated away; when it is small or vanishing, these two stores of lepton 
number do not equilibrate until late times. This effective A-term can induce potentially large 
flavor- violating effects, as will be addressed in chapter lllli 

As discussed in chapter one, B and L are not separately conserved in the early universe 
due to the electroweak anomaly, and as we saw in section ll.2.2[ sphaleron processes associated 
with that anomaly will intermix the two, conserving B — L and violating B + L. Since the 
net lepton number for the universe is zero in Dirac leptogenesis, we would normally expect 
these processes to wash out B, but this conclusion can be avoided if the individual stores of 
lepton number we have generated in charged (s)leptons and (s) neutrinos do not have a chance 
to equilibrate among themselves until well below the scale of electroweak phase transition T^, at 
which point sphalerons have effectively shut off. For simplicity's sake, let us begin by assuming 
that any lepton number stored in sleptons and sneutrinos is rapidly annihilated through an 
effective A-term of the sort described above. In this case, there are only two stores of lepton 
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Figure 2.3: Diagrams corresponding to the leading order contribution to the equihbration of 
left-handed leptons with right-handed neutrinos. These include, from left to right, 
Higgs decays and inverse decays, scatterings off Standard Model fermions and elec- 
troweak gauge bosons through a virtual Higgs, and scattering off a Higgs boson in 
conjunction with the emission of an electroweak gauge boson. The right two dia- 
grams also have t-channel equivalents. The amplitude associated with each diagram 
is proportional to the (small) effective neutrino Yukawa given in equation (12. 9p . 

number and Lj,^ present. Since the £ fields are charged under SU{2) x U{1)y whereas vr 
are not, sphaleron effects will act on Li only, transmuting this store of lepton number into a 
nonzero baryon number along lines of constant B — L^in the manner described in section [1.2.2t 
while leaving L^^ untouched. As discussed above, when Li and L^^^ finally do equilibrate, a net 
baryon number for the universe will already have frozen in, and the universe will end up with 
net positive B and L, as shown in figure 12. 1[ 

The leading contributions to the equilibration rate between left- and right-handed neutrinos 
are shown in figure [231 The amplitude for each of these processes is proportional to the effective 
Higgs couphng term in equation (12.91) . and hence suppressed by {x)/^^i- order for Dirac 
leptogenesis to work, we must ensure that the net equilibration rate does not become significant 
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compared to the expansion rate of the universe until well after the electroweak phase transition. 
The equilibration rate may be estimated on dimensional grounds to be 



|A|W(X)' 
Ml 



re, ~ ' ' ^2 aiJ. (2.12) 



where Qsc is an 0{1) gauge or top Yukawa coupling (see figure 12. 3p and T is temperature. 
Requiring that V ^.^ < H ioi T > T^. leads to the condition 



Mi 



< 10"^ (2.13) 



which is easily satisfied, as it is essentially a recapitulation of the statement that neutrino 
masses must be small. This limit can be translated into a bound on the neutrino mass using 
equation (12.10p 

< 1.74 (sin/?) keV (2.14) 

and thus the stipulation that left-right equilibration happen only well below Tc is satisfied 
automatically when neutrinos are given realistic masses. 

In addition to this constraint, there are several other consistency checks which Dirac lep- 
togenesis must pass in order to be considered a legitimate baryogenesis model. One of these 
is that the Sakharov criterion that the abundances of 0i and 0^ depart from their equihbrium 
values must be satisfied. In order for this to occur, the decay rate To (12.61) must be slower than 
the rate of the expansion of the universe H{T) (11.211) evaluated at the temperature T ~ 
when 01 can no longer be treated as effectively massless and its abundance begins to fall off, 
or in other words 
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This constraint also favors small couplings and large M$^. However, in addition to these two 
requirements, we must ensure that the present value of rj satisfies the bounds in equation (11. ip . 
To do this properly, one must solve the full system of Boltzmann equations, which we do in 
chapter HVl though a rough estimate can be made using the "drift-and-decay" approxima- 
tion [12], in which we assume that the heavy particle decays occur well out of equilibrium and 
that the effects of inverse decays and 2^2 processes where AL^ ^ are negligible. Including 
contributions from both scalar and fermion decays, this gives the result 

Le = ^ = ^K2il) = 7.32 x lO^^e, (2.16) 

s n^g^ 

where nj^™* is the initial number density of 0, K2{x) denotes the modified Bessel function of 
the second kind, evaluated at x, and e is the decay asymmetry given in (12.71) . Since the final 
baryon-to-entropy ratio B (related to 1] hj B = 7]/ 7. OA) generated by sphaleron processes will 
be on the same order {B ^ 0.35L^), this can serve as a rough estimate for B. Thus even if 
equation (12.151) is satisfied, the final baryon-to-entropy ratio of the universe will be proportional 
to e; and from equation (12.71) . we see that for e to be large, either the couplings must be large 
or the splitting between M<i,j and M^^ must be small. 

The upshot of all this is that while there are tensions among the model parameters in 
Dirac leptogenesis, they are not difficult to reconcile — in part because relevant physical scales 
in the theory are determined by the interplay of a large number of model parameters: (M$^, 
S, (x)), and the elements of the complex coupling matrices Ajq, and hi^. Unlike in Majorana 
leptogenesis, where neutrino masses are determined solely from the neutrino Yukawa y^, matrix 
and the masses of the heavy fields according to the see-saw mechanism, neutrino masses in 
Dirac leptogenesis (I2.1UI) depend not only on the superpotential couplings A and h and the 
masses M$., but also on (x). Consequently, the model parameters of Dirac leptogenesis are 
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far less constrained. Of course this versatility comes at the price of introducing an additional 
intermediate scale, though we will show that it is possible to relate (%) to other physical scales, 
for example the Higgs /i-term, in chapter |Vl 

So far the constraints we have discussed have been limited to those which function as 
consistency checks on the model. The real tensions among M<i,j, S, Xia, hia, and (x) are not 
those inherent in the Dirac leptogenesis framework, however, but those that arise when we 
demand that the model respect the full battery of additional constraints from neutrino physics, 
flavor physics, and cosmology. We now turn to address these constraints and their implications 
for Dirac leptogenesis. 



CHAPTER III 



CONSTRAINTS ON THE MODEL 

3.1 Astrophysical Constraints 

3.1.1 Baryogenesis and the Gravitino Problem 

So far, we have seen that a simple toy model of Dirac leptogenesis is capable of yielding a 
nonzero baryon-to-photon ratio for the universe and explaining the observed scale of neutrino 
masses, but we have not yet shown that this scenario is a viable phenomenological model. To do 
this satisfactorily, we must first ensure that none of the modifications we have made disrupt the 
standard cosmology, i.e. that it is compatible with Big Bang nucleosynthesis (BBN), cosmic 
inflation, etc. [38j. Furthermore, we must investigate whether the theory is simultaneously 
capable of yielding a realistic neutrino spectrum and avoiding all current bounds on flavor- 
violation in the lepton sector [13]. At the same time, we must show that the theory is capable 
of reproducing the value of r] observed by WMAP (11. II) once all these constraints are applied. 

There are two primary ways supersymmetric Dirac leptogenesis could potentially disrupt 
BBN. First, constraints on the light element abundances place a bound on the number of 
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additional light neutrino species [44] 



AA^^ < 0.3 



(3.1) 



might cause one to worry that the presence of three additional light, sterile neutrino fields 
in the theory would violate this bound. However, the u^ia are not in thermal equilibrium with 
the bath during the BBN epoch; hence their contribution to AA^ is suppressed by and entropy 
factor 



and thus the bound in (13.11) is respected. The second way in which supersymmetric Dirac 
leptogenesis can impact BBN is that, as in any supersymmetric theory, energy and entropy 
released during late decays of heavy sparticles (gravitinos will be of particular concern) could 
potentially distort the light element abundances from their observed values. This merits an in- 
depth discussion, as gravitino physics can place quite stringent constraints on the leptogenesis 
scale M$^. 

The connection between leptogenesis and gravitino physics, which might not at first seem 
intimately interrelated, occurs through the reheating temperature associated with cosmic 
infiation. Since leptogenesis requires that a thermal population of $i and $i be generated from 
the thermal bath during reheating, we require that < Tr; hence an upper bound on Tr 
translates into an upper bound on M<j,^. On the other hand, if ^3/2 ^ heavy gravitinos 
will also be generated from the thermal bath, and one must take care that they do not cause 
problems for the standard cosmology. There are two distinct varieties of gravitino problem 
that must be addressed. First, as alluded to above, late gravitino decays can disrupt BBN by 
releasing energy in the form of photons and other energetic particles into the system; second, 
as the assumption of i?-parity conservation implies that at least one LSP will be produced at 




(3.2) 
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the end of the decay chain resulting from each late gravitino decay, the potentially large, non- 
thermal population of stable particles produced in this manner could overdose the universe — or 
if the right amount is produced, could make up the majority of cold dark matter (CDM). Both 
of these issues are contingent on the gravitino lifetime r3/2, and for cases where 7713/2 ^ 
(where is the rough scale of the MSSM sparticle masses), this lifetime is approximated 

by m 

Careful analysis of the BBN constraints (see, for example, |16] and references therein) reveals 
that unless 1713/2 > 10^ GeV, Tr cannot be greater than around 10^ GeV, and for 7/13/2 < 
5 X 10'^ GeV, cannot exceed 10^ GeV. In models where 777.3/2 is at the PeV scale, however, 
^3/2 ~ 10~^ s, which implies that gravitinos produced in the thermal bath decay long before the 
BBN epoch (at tuniverse ~ 1 s), and thus there is no gravitino problem of the former type for 
models with 7773/2 at or above these scales. However, since a weakly interacting LSP decouples 
on a timescale t/ ~ 10^^^ s, it will have long frozen out by the time gravitino decay occurs unless 
7713/2 is larger than around 10* GeV; hence in theories with smaller gravitino masses (including 
simple PeV-scale supersymmetry with anomaly-mediated gaugino masses), LSPs produced by 
gravitino decay will be unable to thermalize and the latter type of gravitino issue cannot be 
ignored. 

In order to avoid any complications from late gravitino decay, we require not only that the 
LSP not overdose the universe, but that its surviving relic density Qlsp must be less than (or 
ideally, if the LSP is to constitute the majority of cold dark matter, equal to) the relic density 
of CDM as measured by WMAP [47j . 

ncDuh^ = 0.11 ± 0.01 (WMAP 68% C.L.). (3.4) 
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In general, ^Ilsp will have both a thermal and a non-thermal component, so that ^Ilsp = 
^LSP + ^LSP- '^^^ thermal component filsp niay be ascertained by solving the relevant set of 
Boltzmann equations for the LSP abundance at freeze-out. The results, for the case where the 
LSP is essentially either a pure Wino or Higgsino, are [H] 



Of course it is also possible that the gravitino itself is the LSP, in which case the primary 
concerns are that next-to-lightest supersymmetric particle (NLSP) decays to gravitino do not 
disrupt BBN and that the gravitino abundance does not exceed the WMAP bound (13. 4p . In 
Gauge- Mediated supersymmetry breaking and other theories where the gravitino is exceedingly 
light (with a mass on the order of a few keV or less), the gravitino contribution to the energy 
density of the universe during the BBN epoch [49j is also a concern. It has been shown [HD] 
that these constrains can be translated to a reheating temperature limit Tr < 10"^ GeV, which, 
as we shall see in sect ion |4 . 2 . 2 1 1 ur ns out to be incompatible with Dirac leptogenesis. When the 
gravitino is reasonably heavy, with a mass of 0(100 GeV), the gluino mass is around 500 GeV, 
and the NLSP is Higgsino-like (the best-case scenario) the constraints from NLSP decay and 
f^CDM conspire to produce a reheating temperature bound Tr < 10^ GeV [5T], which will also 
turn out to be problematic for the theory. It thus appears that a gravitino LSP is essentially 
incompatible with Dirac leptogenesis; thus we shall henceforth focus our efforts on models in 
which the gravitino is not the LSP. 

Now we turn to evaluating ^igp- We begin by addressing the regime in which there is 
no significant reduction in Q^gp from LSP annihilations. Assuming for the moment that the 




for Higgsino LSP. 



for Wino LSP 



(3.6) 



(3.5) 
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dominant contribution to the non-thermal rehc abundance comes from late gravitino decays 
and that all the LSPs produced from such decays survive until present day, flisp is given by 

fiS = ""-'f'^^' nMT,/.). (3.7) 

Pcrit 

where pcrit is the critical density of the universe, Tq is the present temperature of the uni- 
verse, and 5^3/2(73/2) is the number of gravitinos per co-moving volume at the characteristic 
temperature T3/2 at which the gravitino decays, which is given by ^52j 

^3/2(^3/2) = 0.856 X lO-^M—^^ — ) ( 1 - 0.0232 In (—^^ — ) ). (3.8) 

Substituting this into equation fl3.7l) yields 



^ 2.96 X 10- (^) (^) (1- 0.0232 m 



(3.9) 



lOiOQeV^ 

In figure [3?n we plot the contours corresponding to the WMAP upper and lower bounds from 
equation (13. 4p on the total relic abundance of a Wino LSP, as well as the simple overclosure 
bound Vtisph"^ = 1, taking into account both thermal and non-thermal contributions, as a 
function of itllsp and Tr. The gently- sloping portion of each contour corresponds to the 
nonthermal abundance in equation 13.91 which becomes significant when Tr is large; the nearly 
vertical portion on the right side of the graph corresponds to the thermal abundance given in l3.5[ 
In the region above and to the right of the WMAP upper bound given in (13. 4p . dark matter is 
overproduced and hence excluded. In the narrow strip between the upper and lower WMAP 
bounds, the thermal LSP abundance and the nonthermal LSP abundance resulting from late 
gravitino decay conspire to reproduce the observed value for Qcdm- In the region below and to 
the left of the lower WMAP contour, the experimental upper bound is not violated and hence 
this region of parameter space is phenomenologically allowed, provided there are other sources 
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Figure 3.1: Contours of Qisph"^ [38j corresponding to the upper and lower bounds from WMAP 
(13.41) . as well as a simple overclosure bound, as a function of the LSP mass itllsp 
and the reheating temperature Tr associated with cosmic inflation for a Wino LSP. 
The requirement that Qlsp not overdose the universe severely constrains Tr, and 
hence the temperature scale of thermal leptogenesis, for a theory like PeV-scale 
loop-split supersymmetry, in which the LSP is particularly heavy. For a Higgsino 
LSP, the results are similar, but the contours move slightly to the left. 

of dark matter to make up the deficit between Qlsp and ft com- These sources could include 
other non-thermally-generated contributions to Qlsp ot contributions, contributions from other 
exotic particles, or some combination of both. It can be seen from figure [3TT] that when niisp 
is large enough that fij^p ~ ft com and CDM is essentially thermal in origin, the ceiling on 
Tr is quite low-around 10^ GeV. In regions of parameter space where rriisp is smaller and the 
majority of CDM is generated non-thermally, Tr may be raised a bit, but is still constrained 
to be below ~ 5 x 10^° GeV. As we shall see in chapter HV] Tr > 10^ GeV turns out to be 
problematic for Dirac leptogenesis (in terms of the final baryon-to-photon ratio generated), 
largely due to the out-of equilibrium condition in (12.151) . This implies that if we want to raise 
Tr above 10^ GeV and still have the LSP relic density dominate ftcDM-, the majority of the 
dark matter abundance (be it from late gravitino decays or something else) must be essentially 
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non-thermal in origin. 

We now turn to address the regime where LSP annihilations do play a role in reducing 
^igp, and thus the upper bound on Tr may be raised. This effect becomes important when 
m3/2 ^ misp- When it is taken into account [53], the non-thermal LSP relic density is modified 
to 

=min((]S°),(]S-")), (3.10) 

where ^isp'^ ^^^^^ density given in equation (13. 9p . and fi^^ is the relic density 

obtained by solving the full system of Boltzmann equations for the LSP. For a Wino LSP, 
^^L5P IS given by 



where xy^ = my/ /rriisp] for a Higgsino LSP, which annihilates far less efficiently, Vt^gp"'^^^ will 
be even higher. 



In figure [3^ we show the relationship between Q^gp'^^^^ and rriLSP for several values of 
m3/2- The horizontal line corresponds to the WMAP upper bound on Q^^^ . From this plot 
it is evident that annihilations are only effective in reducing the LSP relic abundance below 
this bound when 7713/2 is much larger than itilsp- However, when 7/13/2 is increased beyond 
around 10^ GeV, r3/2 becomes short enough that gravitino decay occurs before LSP freeze-out, 
and ^^gp drops to zero regardless of what the ratio of rriisp to 7773/2 is, and nonthermal LSP 
overproduction from late gravitino decays no longer remains a concern. 

While the problems that can arise for small gravitino masses have now been thoroughly 
addressed, the caveats associated with extremely large 7773/2 should also be mentioned. As has 
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Figure 3.2: Here, we show the variation of VL^gp"'^^'^ (the improved expression for the LSP 
rehc density that accounts for the effect of LSP annihilations) with LSP mass for 
several different values of m^i2 [38] • The horizontal line corresponds to the WMAP 
upper bound on VLcdm- For a given choice of 7/13/2, portions of the contour that fall 
below the horizontal line respect the WMAP constraint and are phenomenologically 
allowed, and the reheating temperature Tr be increased beyond the naive upper 
bound from figure I3.1[ Portions that lie above it overproduce the LSP and are 
excluded (for example, when 7713/2 = 10^ GeV, an LSP mass misp > 500 GeV is 
excluded) . 

been shown in [5l] , split supersymmetry models with a large hierarchy between the gravitino and 
gaugino masses can suffer from phenomenological problems associated with the overproduction 
of gluinos, including the distortion of both the CMB and the light element abundances through 
their late decays. Since gluinos will be bound into i?-hadrons at temperatures below the scale 
Aqcd associated with the QCD phase transition, a precise analysis of their decay rate at late 
times has not yet been performed. Still, while a precise ceiling for m3/2 must wait until the 
decay of i?-hadrons is better understood, it is known that this ceiling falls somewhere in the 
7713/2 — 10^° — 10^^ GeV range. For this reason, one should be wary about making the gravitino 
mass arbitrarily large. There are also caveats associated with the gravitino-producing decays 
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of scalar sparticles in models where one or more scalars has a mass larger than m3/2 |55j . 

As discussed above, bounds on the reheating temperature in supersymmetric theories can 
be viewed as bounds on the leptogenesis scale M$^. This is not in any way peculiar to Dirac 
leptogenesis either: in a Majorana leptogenesis model, the constraints still apply with the 
lightest right-handed neutrino mass M^^ in place of M$j. In figure 13.31 we display these 
constraints graphically. For light gravitinos (^3/2 ^ 10^ GeV), BBN limits severely constrain 
Tr, and hence M$^. For slightly heavier gravitinos (10^ GeV < ^3/2 ^ 10^ GeV), there are still 
constraints on Tr from nonthermal decays which are only alleviated (via LSP annihilations) 
when rriLSP ^ ^3/2- For extremely heavy gravitinos, with mass 7713/2 ^ 10^° GeV, decays 
to gluinos become worrisome, and when 7713/2 ^ 10^^ GeV they will almost certainly become 
problematic. 




10- 10*' 10'' 10** 10'" lo'- 
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Figure 3.3: A schematic representing the bounds on the leptogenesis scale from gravitino 
physics as a function of the gravitino mass 7713/2. When 7713/2 ^ 10^ GeV, constraints 
from BBN are quite severe. In the range 10^ GeV < 7773/2 < 10^ GeV, M is still 
limited by nonthermal decays to the LSP. For 7773/2 < 10^° GeV, decays to gluinos 
are a concern. This leaves a window 10^ GeV < 7773/2 < 10^° GeV within which, 
depending on the mass of the LSP, Dirac leptogenesis can be successful. 

While we have not said much about other constraints on M<i>j, but it is not difficult to see 
that substantial model tension arises when the reheating temperature is constrained to be below 
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around 10^° GeV. The out-of-equilibrium condition ( 12.15P demands that Xia and hia be at most 
0(10-^) for ~ 10^° GeV; on the other hand, (E^D and (12161) imply that unless there is a 
large hierarchy among the Yukawa couplings to different sets of $ and $ (so that Aiq, -C A20 for 
all a), decreasing the Aj^ and hia below C(10~^) will yield insufficient baryon number — and we 
have not yet taken into account the effects of 2 <-> 2 processes and inverse decayso Furthermore, 
engineering Xia and hia to be extremely small while the rest of the trilinear couplings in the 
superpotential are C(l) seems to defeat the purpose of leptogenesis, the advantage of which 
was its ability to explain the smallness of neutrino masses without resorting to arbitrary fine- 
tuning. While any precise statement about the value of must wait until after we solve 
the Boltzmann equations governing the development of B and L during the leptogenesis epoch, 
leptogenesis seems to become difficult or contrived when < 10^'' GeV. This means that 
when 771^/2 ^ 10^ GeV and LSP annihilations are ineffective, substantial tensions arise among 
the out-of-equilibrium decay criterion, overclosure bounds related to the reheating temperature 
Tr, the equation that determines the decay asymmetry e, etc., and problems are likely to arise. 
We shall confirm these suspicions in chapter llVi 

The indication that 7713/2 must be quite large in order for Dirac leptogenesis to work suggests 
that the model could be quite successful in the context of split supersymmetry [561 SI]- In 
this scenario, a hierarchy generated between the gaugino masses and the masses of the scalar 
sparticles in which the latter are elevated to a high scale in order to evade unwanted flavor 
violation effects, while the former are kept at the TeV scale or below to constitute dark matter. 
The gravitino mass also tends to be quite large in split supersymmetry, hence it would provide 
a solution to the model tensions that gravitino cosmology presents Dirac leptogenesis. It will 



^One can also get around this by making and M^^ essentially degenerate. We will elaborate on this 
possibihty in section l44l 



37 

soon be made clear that Dirac leptogenesis and split supersymmetry are indeed compatible 
models. 

As a case in point, a particularly simple and phenomenologically interesting split supersym- 
metry scenario in which it is difficult (though not impossible) to get Dirac leptogenesis to work 
is the PeV-scale supersymmetry of [57] , sometimes also referred to as loop-split supersymmetry, 
in which anomaly mediation is invoked in the gaugino sector, but not in the scalar sector. One 
assumes that the messenger fields X responsible for transmitting the effects of supersymmetry- 
breaking to the visible scale are charged under some symmetry, and consequently, while the 
scalar masses are still given by 



mi = c^gi. (3,12) 



,2 _ JX^X 



where c is an 0{1) constant, the term which normally yields dominant contribution to the 
gaugino masses 

! (fe—WW —X^X" (3.13) 
J Mp Mp ^ ^ 

is not longer gauge-invariant. The leading contribution arises at one loop, via the anomaly- 
mediated expression 

If we assume that supersymmetry is broken at an intermediate scale, around ~ 10^ — 10'' 
GeV, then all scalars in the theory, (with the exception of one light Higgs particle) receive 
masses around the PeV scale while the gauginos acquire masses at the TeV scale |58]. This 
model is attractive in its simplicity, and furthermore it is connected to the rich phenomenology 
associated with anomaly-mediated models which includes possibilities for the detection of dark 
matter via observations at the next generation of 7-ray telescopes |59l[60l[6T] and characteristic 
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gluino decay signatures that could be observed at the LHC [62l [63] . 

The problem with getting Dirac leptogenesis to work in loop-split supersymmetry is that 
not only does the theory require m^j2 to be around the PeV scale (which is in itself not 
particularly worrisome), but also mandates a particular relationship between 7713/2 and the LSP 
mass through equation (13 .14^ and the AMSB gravitino mass relation [581 El] 

"I/. = (3.15) 

Figure [3]2] indicates that Dirac leptogenesis prefers a splitting m-ij^lrnLSP substantially larger 
than that dictated by equation 13.141 which means that i^^|p will generally exceed the WMAP 
bound. Loop-split supersymmetry can still be made to work — for example in a situation where 
lepton number production is amplified by resonance effects — but tuning the M$. to this degree 
begs some sort of motivation or additional theoretical machinery. It should be emphasized, 
however, that the splitting between rriLSP and 7773/2 can be much larger [HSIES] in more general 
split supersymmetry scenarios, making Dirac leptogenesis far easier to realize. 

3.1.2 Goldstone Bosons and Symmetry Breaking 



Finally, in addition to the battery of constraints outlined above, Dirac leptogenesis must 
respect cosmological bounds associated with the production of light scalars. As was mentioned 
in chapter [III some new symmetry must be posited in order to construct the Dirac leptogenesis 
superpotential and forbid Majorana masses for the right-handed neutrinos. In Dirac lepto- 
genesis, neutrino masses are the result of the scalar component of the x superfield acquiring 
a VEV which breaks this new symmetry, producing a Goldstone boson or pseudo- Goldstone 
boson, depending on the way the symmetry is broken (which we have not yet specified). Let 
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us first examine the case where the breaking has an explicit component and the relevant scalar 
is a pseudo-Goldstone boson. Constraints on such particles arise from both BBN and cosmic 
microwave background (CMB) considerations [66] as well as from the detection of abnormalities 
in the neutrino flux associated with supernova events [67j, and they can become problematic 
(depending on the mass of the Goldstone boson) when the symmetry-breaking VEV is less than 
around 1 GeV. 

In the case where symmetry breaking is completely spontaneous, a true Goldstone boson 
will result, which these constraints seem to rule out. Of course this assumes that the Goldstone 
boson is a physical state: in the case where the additional symmetry which forbids neutrino 
Majorana masses is a gauge symmetry, the Goldstone boson is "eaten" by the gauge field and 
the relevant constraints become those associated with extensions of the Standard Model gauge 
structure. If it is an Abelian gauge symmetry, for example, one must take care that bounds 
associated with U{1) mixing are not violated and that the gauge theory is free of anomalies 
(which would not be the case given the field content in table 12. ip . We will discuss these 
requirements and investigate their implications further in section 15.31 

3.2 Neutrino Physics 
3.2.1 Experimental Constraints 

In addition to respecting constraints arising from cosmological considerations, in order to 
be phenomenologically viable, a given Dirac leptogenesis model must yield a neutrino spectrum 
that accords with current experimental constraints. The most stringent such constraints come 
from solar and atmospheric neutrino oscillation experiments [68l 169] , and place limits both on 



40 



the mass splittings 



m„ 



m 



Vf,1 



(3.16) 



where the indices a and h label the different neutrino mass eigenstates, and on the mixing 
angles Oab between these eigenstates. The primary connection between the latter and observable 
physics occurs through the leptonic mixing matrix 



U 



MNS 



(3.17) 



where U^'^^ is the neutrino mixing matrix and U^^"^ is the charged lepton mixing matrix. In the 
basis where the charged lepton mass matrix is diagonal, Umns is the neutrino mixing matrix 
and may be expressed in terms of the neutrino mixing angles 6ab as 



U 



MNS 



C12C13 



S12C13 



-S12C23 - C12 5235136*''^^ C12C23 - Si2S23Sl3e "^^^ 



5136^'^^^ 
S23C13 



v 



■S12S23 ~ Ci2C23Si3e ^^^'^ — C12S23 — Si2C23Si3e *''<^^ C23C13 



(3.18) 



where Cab = cos 9ab, Sab = sin 9ab, and 6cp is a CP- violating phase. The present limit cl on the 
Am^^ and 9ab are [TOj 



sin^ 012 = 0.30+^^, sin^ ^23 = 0.50lo.i2, 



sin2 0i3< 0.031, 
Ami = X lO'W, \Aml\ = (2.2+°;^) x lO'W. 



(3.19) 



The smaller of the two mass splittings, Arn^i, is to be identified with the AttIq obtained 
from solar neutrino data (the MSW-LMA solution); the larger, Am|^, with the Am\ from 
atmospheric neutrino data. When we take these constraints and substitute them into the U mns 



^We do not take into account the LSND result which would require an extra neutrino mass eigenstate. In 
the case that forthcoming data from experiments such as MiniBooNE corroborate the LSND signal, it will be 
necessary to extend the neutrino content of our model. 
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matrix, we arrive a set of bounds 



^ .79 -.86 .49 -.58 - .18 ^ 



v 



\Umns\= .30 -.58 .40 -.68 .61 - .80 , (3-20) 
.19 -.46 .50 -.77 .59 -.79 ^ 

However, the constraints in (13.191) say nothing about the sign of the largest mass squared 
difference Am|^. As as result, we are left with two possibilities: the physical neutrino z/3 can 
be either the heaviest of the three mass eigenstates, i.e. mi < m2 <^ (a situation dubbed 
the "normal hierarchy") or the lightest, i.e. <^ mi < 1712 (the "inverted hierarchy"). 

The reason Umns is of particular importance is that in the basis where the charged lepton 
mass matrix is diagonal, it becomes the unitary matrix responsible for diagonalizing the squared 
neutrino mass matrix: 

i'^Ddtag = ^MNS ^^^1 ^MNS- (3.21) 

We can therefore estimate [70] the generic form of the neutrino mass matrix squared, since it 
has to be diagonalized by Umns- In the normal hierarchy scenario, we find that 



{ 1 1 



(3.22) 



where the ^ are small compared to the 0{1) entries and not necessarily the equal to one another. 
Similarly, in the inverted hierarchy scenario, we find 



^^^^ 

e 1 1 



(3.23) 
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In either case, the ratio of the 0{1) entries to the small ^ is constrained to be at least of order 
P32 = Amgj^/Amgi, which, according to f l3.19p . must respect the bounds 



In order for Dirac leptogenesis to be phenomenologically viable leptogenesis model we need 
simultaneously to be able to satisfy the above constraints and reproduce the form of the neutrino 
mass matrix given in f l3.22p or (13.231) . Ideally, we should like to find a simple way of arriving 
at one of these matrix structures that is motivated by theoretical considerations as well as 
data pressures. Of course this solution must also be compatible with successful baryon number 
generation, the astrophysical constraints discussed in section \3A\ etc. 

3.2.2 The Flavor Structure of the Trilinear Couplings 

As we saw in chapter [Tll the neutrino mass-squared matrix in Dirac leptogenesis is given by 



As long as A and h are completely generic and there are at least three sets of $ and $, it is 
apparent that a matrix of this form can yield an arbitrary neutrino mass spectrum. On the one 
hand this is good, for it means that the theory is perfectly viable, in the sense that there exists 
some set of parameters that will satisfy the battery of constraints given in equation (I3.19P : 
on the other hand, this arbitrariness comes at the price of introducing many additional free 
parameters, whose relative values must be determined by some additional underlying physics. 
As a first step toward understanding what that additional physics ought to involve, let us 
examine the spectrum of a simplified model containing only a single pair of heavy fields $i and 



20.0 < P23 < 39.7 . 



(3.24) 




2 2 (or 3) 3 



(3.25) 
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$1 of mass M$j, or alternatively, a theory in which <C M<j,. for alH > 1. In such cases the 
mass-squared matrix is proportional to an outer product of the family-space vectors Aiq, and 
thus its eigenvalues are 

3 



m 



(3.26) 



Here, two of the physical neutrinos are massless. For each additional set of $ and $ with a mass 
similar to , an additional neutrino acquires a nonzero mass. Thus in "short-suited" models 
where <^ Mq>^ and M$2 -C M$. for alH > 2 (i.e. where there are effectively only two sets 
of (f) and (f) involved in determining the neutrino spectrum), one neutrino mass eigenstate is 
massless and a hierarchy will exist between the other two, determined by 5 = M^^/M^j. One 
can then arrange for the other two neutrino masses to take the experimentally observed values 



V2 



Am2i and m^, ^ Am|]^ by tinkering with 5 and the structures of A and h. Conveniently, 



leptogenesis in such models is well-approximated by the toy model considered in section [TTl in 
which there were only two sets of $ and 

Let us now examine the effect of A and h on the neutrino spectrum in models with this 
sort of hierarchy among the heavy particle masses. It is apparent from equation (13.251) that the 
matrix structure of ml is primarily determined by h rather than A. In the limit where 5 -C 1, 
this matrix becomes 

l^llP ^11^12 ^11^13 ^ 

= sin/3) 



m,. 



\ 



\hn\' 


huhi2 




^12^11 


\hi2? 


/l*2^13 


^13^11 




\hi3\' 



(3.27) 



If for some reason the off-diagonal elements of h are much smaller than the diagonal elements 
(i.e. hi2, hi3 ^ hii), this matrix will take form (13.221) that reproduces the normal hierarchy 
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solution. If hi2 and hi^ are roughly of the same order, so that hu ~ /113 = h, equation (13.27P 
becomes 



\h\ 



17I 



Ml 



\ 



e 




A 




1 


1 




1 





(3.28) 



where ^ = hu/h. Here we have implicitly assumed that contributions of 0{hii/h) dominate 
over the 0{6) contributions. If this is not the case, and in particular if hn = 0, the matrix 
structure becomes 



^5^ S 



(ml) oc 



v 



(3.29) 



6 1 1 
6 1 1 

Again the generic structure yielding the normal hierarchy scenario is obtained, but in this case 
the value of 6 fixes the ratio p23- In any case, what is significant here is that imposing the two 
simple hierarchical requirements, <^ ^$2,3 ^11 ^ ^12 ~ ^13 generically gives rise to 
the correct neutrino phenomenology. 

These considerations suggest that we ought to look for some setup that will ensure that the 
diagonal element hu is small compared to the off diagonal terms hi2 and hi^. One particularly 
simple way to do this is to posit an antisymmetry condition on the matrix h, in which case the 
neutrino mass matrix will take the form given in (13.291) . with small entries of 0{6). This option 
is of particular interest because there is also a theoretical motivation for it: antisymmetric 
Yukawa matrices emerge quite naturally in certain grand unified theories and in models with 
non-Abelian flavor symmetries [TH [731 [TH [72] . We shall put off discussion of such matters until 
section 15. 2^ where we provide a theoretical motivation for the hierarchies among the M$- and 
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elements of A and h. 

3.2.3 Constrained Hierarchical Dirac Leptogenesis 



Motivated by the preceding remarks, we now define a particular model that ought to be 
able to yield a phenomenologically acceptable neutrino spectrum. We define constrained 
hierarchical Dirac leptogenesis (CHDL) [38] as the setup in which: 

1. The mass matrix (M$) is real and diagonal. 

2. The coupling matrices A and h are antisymmetric 

3. The large mixing angles in the neutrino sector are result from the smallness oi 6 = 
In this model, the antisymmetry of A and h allows us to parameterize them in the manner 



A = / 



V 



1 a2 
-1 as 
-a2 —as 



h 



h = f 



J 



-h 6s 
-62 -b3 



(3.30) 



which is convenient when the a^ and bi are all roughly 0{1). Since the assumption of a hierarchy 
among the M$. leads to the neutrino mass-squared matrix of the form fl3.22p . we expect that as 
and 63, the effects of which show up only at the 0{6) level, will be less tightly constrained than 
the rest of the a^ and 6j, which contribute to the leading term. This is in fact the case: if ai, a2, 
bi, or 62 deviates significantly from one, the neutrino spectrum cannot satisfy the constraints 
in (13.191) . It is therefore appropriate, since the value of / is unimportant as far as this set of 
constrains are concerned (any rescaling of / can be compensated for by a similar rescaling of 
(x)), to analyze constrained hierarchical models as functions of as and 6s alone. 
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In figure 13.41 we show the region of viability in 03 — 63 space for two different values of 6: in 
the left-hand panel, we set S = me/m^ (and Mq,^/Mq,^ = 171^/111^-) as required in the minimal 
version of CHDL discussed above; in the right-hand panel, we set 5 = 10^^ We consider a given 
combination of 03 and 63 to be phenomenologically viable if there is any combination of the 
remaining and bi for which the combination simultaneously obeys all the neutrino oscillation 
constraints in (13.191) . The entirety of the shaded region shown in each panel, including all 
differently-shaded bands, is permitted by these constraints: the bands represent of contours 
sin ^13 the value of which will be measured or constrained in future neutrino experiments. This 
plot demonstrates two important features of the Yukawa matrices in CHDL: first, it is indeed 
possible to satisfy the neutrino oscillation constraints for 6 = rrie/mn] second, while 63, like 
most of the other and bi, is constrained to lie fairly close to 1, is permitted to be quite 
large when 6 is small. In the same figure, we also show contours for the value of sin6'i3, the 
value of which will be measured or constrained in future neutrino experiments. It is seen that 
the value of sin ^13 increases with increased 03 until reaching its maximum experimental bound. 

Let us now take a moment to address how these results affect leptogenesis. Since we are 
assuming that M$3 ^ , M$2 , the formula (12.81) for the decay asymmetry e tells us that 
lm{Xl^X2ahlph2i3Mq>-^M^^) will vanish (when it involves diagonal elements of A or h) unless 
(3 = a = 3. This means that 

e oc Im (020362^3) r (3.31) 

1 — 

in CHDL; the two panels in figure [3^ then show that for a given 6, the largest amount of left- 
handed lepton number Ll (and therefore the largest baryon asymmetry) will be obtained when 
as = a'^°'^{S), where a^'^^{6) is the maximum possible value of as for a given 6 consistent with 
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Figure 3.4: Here, the regions of 03— 63 space (see equation fl3.30p for a description of the Yukawa- 
matrix parametrization used) for which all constraints on neutrino masses and 
mixings fl3.19p are simultaneously satisfied for some combination of the remaining 
aj and 6j are shown for two different values of 5 [3B] ■ Additionally, contours depicting 
the ranges for S13 = sin ^13 (which depends primarily on 03, but varies slightly with 
the remaining and 6j) are shown. In the right panel 6 = 10~^; in the left panel, 
6 = irie/m^ = 4.83 x 10~^ The plots reveal that while 63 is constrained, along with 
most of the other and bi, to lie reasonably near 1, a™'^^(5) can be quite large and 
increases with decreasing 6 = 10^^, making it easier to obtain a realistic baryon- 
to-photon ratio t]. For S = 10~^, a'^'^ ~ 4.5; for b = 4.83 x IQ-^, a^'^^ ~ 95. In 
each panel, the configuration that yields the greatest decay asymmetry e is marked 
with an arrow. 

neutrino masses and mixings. It is interesting to note that since the maximum experimental 
value of sin6'i3 sets the value of a^""^ , then the maximum baryon asymmetry will be obtained 
when sin ^13 acquires its maximal experimental value, which we take to be sin^^^ ^13 = 0.031. 
For 6 = 10~^, the result is a™'^'^(l/10) ~ 4.5, as indicated in the right panel of figure 13.41 
It is also of note that a™'^^(5) rises sharply as 6 is decreased, as indicated by the result for 
6 = 4.83 X 10^'^ (a value which, as we shall see in section is predicted by certain theoretically- 
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motivated extensions of CHDL) shown in the left panel, in which case a'^"'^ ~ 95. 

One of the assumptions we have been implicitly making here is that we have maximal CP 
violation in the decays of the fields $ and $ (i.e. the overall phase in the product 020362^3 must 
be 7r/2). In this case one can obtain a value for the effective CP violation in the lepton sector. 
This can be defined in a phase invariant way in terms of the quantity J = Im(f/i 2^/22 ^23^^*3) 
[75| [76] . where Uij = {UMNs)ij- Taking for example the two points that will generate the most 
baryon number in the left and right handed panels of figure 13.41 (marked with an arrow) we find 
J ~ 0.034 for the point in the left panel and J ~ 0.030 for the point in the right panel (since 
the maximal value for J is |f/i2f^22^23f^r3l, which for these points yields a value of J ~ .038, 
the CP- violation here is close to maximal). It could be interesting to do a more detailed study 
of the issue of linking more generally the effective CP violation in the lepton sector to the CP 
violation in the interactions of the heavy fields $ and Nevertheless, we will not pursue this 
issue further in the present work. 

We have shown here that it is possible to construct Dirac leptogenesis scenarios that are 
capable of satisfying current constraints on the light neutrino spectrum. In CHDL, the set 
of parameters relevant to neutrino physics and to leptogenesis is quite small, given the anti- 
symmetry condition on the two coupling matrices A and h: the leptogenesis scale M<i,j, the 
intermediate scale (x), the trilinear matrix entries 03 and 63, the mass ratio 5, and the overall 
coupling strength / are essentially the only free parameters in the model. In fact, even this 
represents an overcounting: equations fl3.25p and fl3.19p imply an additional constraint on these 
parameters. For example, take the case where 5 = 10^^ and 03 and 63 have been set to the 
values most advantageous for baryogenesis, 03 = 4.5 and 63 = 2.2 (see figure . In this case, 
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the constraint becomes 



Pix) 



sin/3 = 1.009 X 10 



-13 



(3.32) 



and the model contains only five free parameters (and that's including 63, which is rather 
stringently constrained). We will generally use this constraint to eliminate (x) or /, depending 
on which is most convenient for the purpose at hand. 

While we again emphasize that it is certainly not the only scenario for obtaining a realistic 
neutrino spectrum in Dirac leptogenesis, CHDL is a particularly simple and predictive model — 
and it will be the one on which we will focus our attentions from this point forward. It 
still remains to be seen whether CHDL (or indeed any Dirac leptogenesis model) can also 
simultaneously satisfy the constraints from gravitino physics and yield a realistic baryon number 
for the universe, however, though as we shall see in chapter HVl thermal Dirac leptogenesis will 
indeed turn out to be workable in a variety of models. 

3.3 Soft Masses and Flavor Violation 
3.3.1 Flavor Violation in Supersymmetric Models 

Another requirement which Dirac leptogenesis must satisfy in order to be considered viable 
is that the theory must be compatible with present constraints on fiavor violation in the lepton 
sector. The most stringent bounds come from measurements of the branching ratios for fiavor- 
violating decays and conversions of heavy leptons, such as /i ^ 67, r — ^ /X7, fj, — > eee and 
fiA eA. The current experimental limits on the 2-body decay processes are [77] 



BR{n 67) < 1.2 X 10 



,-11 



(3.33) 



BR{r 



fi-f) < 1.1 X 10 



6 



(3.34) 
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The MEG experiment [78] is expected to improve on bound on — 67 by several orders of 
magnitude in the near future, bringing it to 0(10^^^ — 10~^^) or lower. Other related projects, 
such as PRIME [79] (sensitive to fiA eA conversion), are expected to go online over the 
next few years, and projects have also been proposed [80] that would lower the bound on 
BR{t /i7) to 0(10-9). 

Supersymmetric theories with generic soft parameters tend to result in unacceptable levels 
of flavor violation due to flavor misalignment between quark and lepton mass eigenstates on 
the one hand, and squark and slepton eigenstates on the other. There are a variety of ways 
to address this problem, the most common one being the assumption of soft mass universality. 
Here, one assumes that whatever mechanism gives rise to supersymmetry breaking at some 
high scale M (the Planck scale, the GUT scale, etc.) results in a set of squark and slepton 
soft mass-squared matrices that are flavor-blind and diagonal and A-terms proportional to the 
Standard Model Yukawa couphngs, i.e. that 

niQ = m| = mil, m| = mil, 

m| = m|l, m| = m|l, m^=m^l, 
= AuYuy = AuYd, = AuYe, = AuYu, (3.35) 

then all squarks and sleptons become degenerate in mass and can be freely rotated into one 
another up to A-term-induced mixings, which will only be large for sfermions of the third 
generation. Flavor- violating effects would then be expected to be very small. 

In Dirac leptogenesis, even if soft masses are universal at scale M, as they are run down 
from M to the leptogenesis scale fiavor-off-diagonal terms will be generated by quantum 
corrections, due to the nontrivial matrix structure of the trilinear coupling matrices A and h. 
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It can be shown that, under assumption of universal supersymmetry breaking, this is not a 
damning problem for the MSSM (for a brief review, see [Slj). However, since the lepton sector 
of the Dirac leptogenesis superpotential (12.11) is modified from that of the MSSM and includes 
new coupling matrices with nontrivial flavor structure, we must ensure that any off-diagonal 
contributions to the slepton mass matrices mj^i^, ^"rr and rn\j^ generated by these modifications 
do not result in an unacceptable level of flavor violation. The renormalization group evolution 
(RGE) equations for parameters appearing in a general superpotential soft supersymmetry 
breaking Lagrangian are well known [82], and the RGE equation for a soft mass is given by 



J 167r^ 



+al^a^^'^ - SgiCa{maV5l + 2^7^(T«)lTr[T'^m^]J , (3.36) 

where the ajj^ are soft A-terms, Ma are gaugino masses, Uijk are trilinear couplings appearing 
in the superpotential, and a sum over gauge groups is implied. The Caij) are quadratic Casimir 
group invariants, defined in terms of the generators T° by the relation 

C^m = (T'^Ty,. (3.37) 

Assuming the Universality condition at the high scale and that all soft A-terms are equal to 
the relevant Yukawa coupling multiplied by the universal soft super symmetric mass m^, one 
can simply estimate the flavor violating corrections to the mass matrix by integrating the RGE 
equations iteratively [83l [M] , and using this method, one obtains off-diagonal contributions to 
the slepton masses 5m|^ and 5m\j^: 



StjiIl ~ "2^^"^ (^M^J ^^"^^^"^^ "^^-^^^ 
SmlR ^ -^\n(^) X*^Xifsml (3.39) 
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Another off-diagonal scalar mass term arises from the effective A-teim in (12. lip once electroweak 
symmetry breaking occurs in models where the F-term of x acquires a VEV (F^) 7^ (a 
corollary in many mechanisms where its scalar component obtains its VEV). After electroweak 
symmetry breaking, this term results in a contribution 

6ml^ = hlX,p^vsm/3 (3.40) 

to the sneutrino mass matrix which mixes left-handed and right-handed sneutrinoqf . 

In order to examine the effect of these mixings, the full mass matrices for both the charged 
sleptons and sneutrinos must be taken into account. For simplicity, we will continue to assume 
that the leading soft breaking sector is flavor diagonal and universal with a common scalar 
mass rris- The resulting additional contributions to the slepton mass squared matrices, given 
by equations (13.381) . (13.391) . and (I3.40p . can thus be expressed in terms of the 3x3 submatrices 



5m|^, Srnj^j^, and Smj^j^ as 



6ml^ ; \ / 6ml^ 

drni, = I 1 5m 



(3.41) 

The only contribution to the charged slepton mass squared matrix comes from 5m\j^ (no off- 



0:0/ V (5mi^)t 



^^RR 



diagonal terms will be generated radiatively in the RR part of charged slepton matrix mj^ since 
the Yukawa coupling matrix Yg can always be chosen diagonal at the high scale), while the 
sneutrino mass squared matrix receives not only additional flavor mixings among left-handed 
and among right-handed sneutrinos, but also an effective A-term from 5m\ji which intermixes 
left-handed and right-handed sneutrinos. 



^The CP-violating phases in A and h required for leptogenesis can induce new phases in the slepton mass 
matrices during the RGE running or through the mixing term of Eq. (|3.40p . We have checked that the effect of 
these phases is too small to reach the experimental bounds on lepton EDM's, within the region of parameter 
space considered here in which leptogenesis is successful. As for the rest of phases of the MSSM, they are 
assumed to be small enough to avoid violating these same experimental bounds. 
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In CHDL, where we have a specific fiavor structure for the matrices A, h and M$, one 
has a specific prediction for fiavor mixing among sleptons once the electroweak and hidden 
symmetries are broken. Using the neutrino mass constraint (13.321) . we can express the overall 
dependence of the slepton mass terms in fl3.38p , fl3.39p , and fl3.40p on the relevant mass scales 
in the theory: 

5ml^ oc oc ^ (3.42) 



(5m? 7? oc TT- 



It should be noted here that the proportionality constants for the bottom two equations are 
not dimensionless: the ones associated with 5m\i^ and 6m\^ each contain a factor of m^ and 
have mass dimension [m]^, while the one associated with 5m|^ contains a factor of {F^)v and 
has mass dimension [m]^. 

3.3.2 Lepton-Sector Flavor Violation in Dirac Leptogenesis 



1 




Figure 3.5: Feynman diagrams giving the two leading order contributions to the fiavor-changing 
process i~ ij'j due to sneutrino (left diagram) and charged slepton (right dia- 
gram) mass mixings. 



The effective interaction leading to lepton fiavor violating decays of the form £j ij^y, 
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where li and ij are charged leptons, can be written as 

J = iem,^Ui{q-p)a^pq'' {A^^Pl + A'^Pr) Uj{p)e*{q), (3.43) 

where q and p are the momenta of the photon and the outgoing lepton £j respectively, and mi. 
is the outgoing lepton mass. The resulting decay rate is 

r(/7 _ l~ ^) = ^mlilA^l' + \A^\'). (3.44) 

The leading contributions to the amplitudes A^ and A^ appear at one loop level and are shown 
in figure [3751 They involve both a sneutrino (and chargino) mass eigenstate and charged slepton 
(and neutralino) mass eigenstate running in the loop. These amplitudes were computed in [83] 
for a general MSSM scenario and we have included the modified expressions which account for 
the presence of light right-handed neutrino and sneutrino fields in appendix 1X1 




Figure 3.6: Feynman diagrams for the leading-order process involving 5m\^ (left diagram), and 
for the leading process involving Sm\^ (right diagram), with sneutrinos running in 
the loop in the mass-insertion approximation. Note that any process involving 
5m|^ necessarily involves two mass insertions, and any one involving Smj^j^ (given 
by the diagram on the right with an additional Sm^j^ insertion) necessarily involves 
three. 



Since the relationship between the Dirac leptogenesis model parameters /, (x), and M$^, 



which enter into equation fl3.44p via the masses and mixings defined in appendix |Al and the 
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flavor-violation rates ( 13.44p is somewhat obscure, let us take a moment to make it a bit more 
transparent. Consider an exclusion contour in parameter space corresponding to the /i — 67 
bound (which will turn out to be the more stringent of the two) in equation fl3.33p . Along 
this bound, the branching ratio for BR{fi 67) is by definition quite low and flavor-violating 
effects will be small. It is therefore valid to use the mass-insertion approximation there and 
treat 5m\j^, ^^"iw ^^rr small corrections to the slepton propagators. Let us focus on 
corrections to the sneutrino propagator, which enters in the left diagram in figure 13. 5[ since 
it can receive corrections from all three of these. The leading contributions involving each of 
6m\]^ (left panel) and 5m\^ (right panel) are shown in figure 13.61 (the leading contribution to 
5m\^, which would require three mass insertions, is the least important of the three, and is not 
shown). Since there is no coupling between leptons and right-handed sneutrinos, corrections 
from 5m\j^ and 5m\j^ only appear at second and third order in the mass insertion expansion, 
respectively. Therefore, if there is no substantial hierarchy among these three sets of mixing 
terms, mixings between left-handed sleptons provide the primary source of flavor violation. In 
the approximation that all slepton, chargino, and neutralino masses are roughly on the same 
scale msoft, the branching ratio may be estimated as 



In interpreting this result, let us treat and (x) as free parameters, treat S, 03 and 63 as fixed, 
and eliminate / using the constraint (13.321) . Equation (13.451) tells us that contours of branching 
ratio in the M,j,^ — (x) plane correspond to contours of 5m\j^. According to equation (I3.42p . 
5m\^ = ciM$j/ (x), where ci is a dimensionless proportionality constant with dimension [m]^. 




(3.45) 
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so the exclusion contour associated with left-left mixing takes the form 

lnM^ = ln(x) + CiL, (3.46) 

where Cn = ~" IiiI'^^^ll/'^i) an as yet undetermined constant. 

In the absence of any large hierarchy between 5m\i^, ^'^\r) ^"^^ ^^"rr the oblique contour 
defined by equation (13.461) is the only important one. However, there is no a priori reason why 
such a hierarchy should not exist. The 5Tn\^ contribution (I3.40p is proportional to {F-^, which 
has little relevance to baryogenesis other than that it serves to equilibrate the left- and right- 



handed sneutrino fields and is essentially unconstrained. As was pointed out in section 12.21 
a/ {F) can potentially be quite large (10^ GeV or higher), and if this is the case, contribution 
from 5m\j^ could be as important as those from 5m\i^. Let us assume for a moment that this 
is the case and examine the constraints related to Smj^j^ and Smj^j^ together. In regions of 
parameter space where (x) is small, we now have 

BRif, ^ 67) oc , s r = constant (3.47) 

along any exclusion contour. Equation (I3.42p tells us that, 

Smlji = C2-, (3.48) 

where C2 has mass dimension [m]^. The associated contour is therefore given by 

ln(x) = Clr, (3.49) 

where Clu = ln(c2/5m|^). As for equation (I3.42p implies that it becomes impor- 

tant (and in fact the dominant contribution to the right-handed sneutrino mass) in regions 
of parameter space where (x) is small and M<i>j is large. In this regime the mass insertion 
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approximation can no longer be used, but we can use the approximation m?^ ~ ^^'rr the 
sneutrino propagator in equation (13 .47^ . The amphtude for the process in the left panel of 13.51 
then implies 

where C3 is a proportionality constant with mass dimension [m]^. This implies yet another 
oblique exclusion contour corresponding to the line 

\nM^, = -\n{x) + CRR, (3.51) 

where 

This line runs parallel to the one from left-left mixing (I3.46p . As discussed above, it will 
generally be the case that the contour provides a more stringent bound than 5m|.^, and 
hence this contour can generally be ignored. 

Taken together, the contours determined by equations (I3.46P and (I3.49p suggest that unac- 
ceptable amounts of flavor violation will occur in regions with large and small x- It is not 
yet obvious exactly how large and how small x must be before problems arise (or to put it 
another way, what the precise values of Cll and Cm are) in a given Dirac leptogenesis model, 
however, nor is it obvious that these constraints are compatible with those from gravitino cos- 
mology, neutrino physics, etc. These issues will be addressed in the next chapter through a 
careful numerical calculation which takes into account the full formulae for the masses, mixings, 
et al. given in appendix lAl 



CHAPTER IV 



EVOLUTION OF THE BARYON ASYMMETRY 



4.1 Boltzmann Equations 



Let us now turn to the numerical calculation of rj and the solution of the full Boltzmann 
equations. This is slightly more comphcated for Dirac than for Majorana leptogenesis because 
in the latter, it is only necessary to keep track of the overall lepton number. In contrast, the 
former involves the creation of several distinct stores of lepton number with different properties. 
The heavy fields aside, in Dirac leptogenesis there are six particle species charged under lepton 
number {ur, Dr, £, £, and the right-handed charged lepton and slepton fields cr and en), and 
thus six individual stores of lepton number to keep track of: L^, Lj,^, L^, L^^, Lg^, and Lg^. 
Some of these stores are positive and some negative, some of the particle species are involved 
in sphaleron interactions while others are not, and so forth, so it might appear necessary that 
we keep track of each field and each store individually. Including an equation for the overall 
baryon number of the universe B and accounting for the dynamics of the heavy fields in the 
$1 and $1 supermultiplets raises the total number of equations in the Boltzmann system to 
twenty-one. The situation can be greatly simplified, however, by noting that the fields £, £, 
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ept, and cr participate in SU{2) and/or f/(l)y gauge interactions, which should be sufficiently 
rapid (compared to other processes relevant to leptogenesis) that these species will always be 
in chemical equilibrium with one another. As a result, any lepton number stored in any one 
of them will be rapidly distributed among L^, L^, Lg^, and Lg^ in proportion to the relative 
number of degrees of freedom of each respective field. These fields then compose a distinct, 
"visible" sector of the theory with an aggregate lepton number Lmg. 

As for the remaining two species charged under lepton number, there are two possible 
choices. While the vr, by construction, have no interactions with the visible sector and thus 
can be seen as forming a "hidden" sector which only equilibrates with the visible sector fields 
at late times through the effective neutrino Dirac Yukawa, no such requirement exists for the 
right-handed sneutrinos vr. When [F^ is large, the vr fields have additional, rapid interactions 
with the visible sector fields (and in particular t) through effective A-terms (12. lip and should 
thus be considered part of the visible sector. As shown in section [3.3.2l an F-term VEV of order 
(F^) ~ 10^ GeV is still permitted by flavor violation constraints in theories with weak-scale 
slepton squared masses, and in split supersymmetry scenarios with much heavier sfermions, 
even this bound no longer applies, so rapid sneutrino equilibration is far from excluded. On 
the other hand, when (F^) is small, the vr fields decouple and become part of the hidden 
sector. In either case, since none the hidden sector fields interacts via electroweak sphalerons, 
an aggregate lepton L^ij^ can be defined for the hidden sector, which is just the sum of the 
lepton numbers individually stored in each of its constituent fields. We will henceforth refer to 
these two situations as the \axge-{F^ and small- (F^) scenarios. 

The definitions of L^s and L^ij^ simphfy our task considerably; in the limit of rapid equili- 
bration within the visible sector, only three Boltzmann equations are required to describe the 
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evolution of lepton and baryon number asymmetries of the light fields in Dirac leptogenesis: 
one for Lyis, one for Lhid, and one for the overall baryon number B. In addition to these, it 
turns out (see appendix [B]) that only three additional equations are needed to describe the 
dynamics of the heavy fields: two to describe the evolution of the individual lepton numbers 
L^^ and L^_ stored respectively in (pi and (pi, and one to track the abundance of one of the 
heavy fields (we choose ). Our Boltzmann system has now been reduced from twenty-one 
equations down to a far more manageable six. 

In addition to the rapid gauge interactions that equilibrate the fields in Lyis, there are 
a variety of additional processes which we must take into account. First, we must include 
the decays and inverse decays of the fields in $i and $i. This introduces a pair of rates 
Tl = r(0 —^i + x) and Tfi = r(0 — > z^j^ + -f^^), which are related to the overall (p decay rate 
To given in (12. 6p by the defining relation 

rD = rL + rn. (4.i) 

Second of all, we must include interactions involving virtual (pi and (p^ fields (and their fermionic 
superpartners) that can transfer lepton number between L^i^ and Lhid- The dominant contri- 
bution to the transfer rate comes from 2^2 processes, and since this rate, which we will 
call T2^2, will be suppressed by inverse factors of the leptogenesis scale for T ^ M^^, these 
processes cannot be considered rapid compared to the gauge interactions. Third, since annihi- 
lation processes second-order in the heavy fields can serve to reduce the abundances of these 
particles during the leptogenesis epoch [HSIES], we should include them as well, with rate Ta- 
Fourth and finally, we must include sphaleron processes with rate T^ph, which will interconvert 
B and L^ig. 

We are now ready to write the full set of Boltzmann equations governing the evolution of 



61 



baryon number B in the early universe, a full derivation for which is provided in appendix [HI 
along with explicit definitions for all quantities used therein. Since these equations will differ 
slightly between the large-(F^) and small-(F^) scenarios, due to the differing constitution of 
Lyis and Lhid, they need to be written down separately for each case. We will concentrate on 
the large-(F-^) scenario, in which they are 
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in terms of the variable z = M^^/T. Here, the inverse decay rates (r^)^^, (Ti) and (F^) 



(4.5) 
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are defined by 



(r«)„ = ^ fl^l r«, (4,10) 



where F^) is tlie total decay widtli of (pi given in equation (12.61) . n^^ is tfie equilibrium number 
density of 0i (and 0^, etc.), and the quantities Ti and represent the partial decay widths 
for (j) — > vjiH (or cf) i^rH) and (p —>■ ix (or — * ix)- The ratio of modified Bessel functions 
Ki{z) and K2{z) appearing in these rates is a result of averaging over time-dilation factors: 
see ( IB.23P in Appendix [Bl We will examine this scenario from this point forward. 

4.2 Numerical Analysis of the Boltzmann System 
4.2.1 Calculation of Rates 



From this point forward we will focus our attention on the large-(F-^) case and proceed to 
solve the Boltzmann equations (14.21 - 14.7( 1 numerically. In order to do this we must first calculate 
the relevant rates appearing in these equations in terms of the model parameters M$^, 6, a^, 
63, and / (or alternatively (x))- We already have expressions for Ygph (11.201) and Yd (12. 6p . 
from the latter of which and Tn can be obtained trivially. This leaves only the rate F2^2 
for processes that shuffle lepton number between L^s and Lhid and the heavy field annihilation 
rate F^. We now turn to address each of these in turn. 

While there are a large number of 2 2 processes which shuffle lepton number between 
different particle species, only a few will transfer it between L^s and Lhid- The rest, which 
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collectively serve to assist the rapid SU{2) x U{1)y gauge interactions in equilibrating lepton 
number among the fields in the Lms sector, can be ignored. The relevant s-channel diagrams 
are pictured in figure 14.11 In addition to these, there are contributions from the t-channel 
transforms (two per diagram) of these diagrams. In order to evaluate diagrams containing 



virtual heavy fermions, we define the Dirac spinor 

[All] 



( (I ^ \ 



Di 



where (^/'^.Ja and are the Weyl spinor components of the $j and $j superfields. Numerical 

calculation of the thermally averaged cross-sections for the diagram pictured on the left in the 
top row of figure 14. H which involves two Yukawa-type couplings of a scalar to two fermions, 
yields {(Jv)iap ~ af^\\iXWi3?, where [SSj 

{2Y) ^ ^Q_2 ^2) 



For the diagram on the right in the first row of figure 14. ![ which includes one Yukawa- type cou- 
pling and one trilinear scalar coupling proportional to M$ . , the result is very nearly temperature 
independent and well approximated by {av)iai3 — o^i^'^'^^^ \\a\'^\hii3\'^ , where 

-."™'-0-5x-r^. (4-13) 

The contribution from each of the two diagrams in the second row of figure BTTl which involve two 
Yukawa couplings and a mass insertion from the heavy fermions, is equal to that in (14.131) . These 



interactions dominate among 2^2 processes. The diagram in the bottom row of figure 14. 
which involves a trilinear scalar coupling to the down-type Higgs, may be approximated by 

{crv)iafs ^ a-™^|Aiap|Ai^p, where 

« -A., (4,14) 
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The constant of proportionality in this equation is and is thus suppressed relative to the 

rate given in (14. 13p by yu/M$^. Here, we will assume that fi is several orders of magnitude 
smaller than all the M$^, and therefore the effect of these processes can be neglected. 



/ \ / 



T>R ^R X^ ^ ^R 



Figure 4.1: Diagrams for 2^2 s-channel processes which transfer lepton number between 
LaggVis (the aggregate lepton number in the sector comprising the fields £, £, vr, cr, 
and cr, which are assumed to be in chemical equilibrium with one another due to 
rapid gauge and (F^)-term equilibration interactions); and Lhid (the lepton number 
stored in right handed neutrinos). The two t-channel interactions associated with 
each diagram appearing above must also be included in calculating the full thermally 
averaged cross- sect ion. 



Taking into account contributions involving virtual fields in the $2 and $2 supermultiplets, 
as well as those in $1 and $1, we find the total interconversion rate between L^ig and L/^j^ in 
the large- (F-^) scenario to be 



r2^2 ^ 3n,5^5^(ar)|Ai.nA,,r + af^)|A2.r|A2, 

a (3 



a p 



+ 3n,5^5^ (ar^|Ai„nAi,r + af^)|A2„r|A2,r) , (4.15) 

a p 

where we have assumed an equilibrium number density for all non-leptonic light species (e.g. 
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x) involved. 

The annihilation rate for the heavy fields Ta, which is associated with second order processes 
of the form — > ij and 010^ — > ij (and miscellaneous supersymmetrizations thereof), is 
most readily expressed in terms of the reaction density, given by the general expression 

where T is temperature, s is the usual Mandelstam variable, and a{s) is the total reduced cross 
section for annihilations of etc. into light fields. This is defined by the formula 

^(«) = ^/*'EI-^^WI'^^' (4-17) 

where both t and s denote the Mandelstam variables. The hmits of integration are given by 
t± — — s(l =pr)/2, with r defined below. The relationship between the reaction density for 
the annihilation of two particles i and j with number densities ni and rij, the thermally-averaged 
cross-section {(t\v\) for annihilation, and the annihilation rate F via 

j^ninj{a\v\) = -^r. (4.18) 

Hence, in the case under consideration here, the annihilation reaction density7A between 0i 

and (pl is given by 

7^ ~ -r|e (4.19) 

in the (very good) approximation than L<^^ <^ F^c . 

In supersymmetric Dirac leptogenesis, the total reduced cross-section 7^, including all rel- 
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evant decay processes, is [43 



±.tot 1 

'^SUSY - 



16n 



4 / / 66\ ^ / 16 16 A,fl + r 
+(72^ 32 + - r + 3 -- - - + 9 In ' 



where x = s/M'^^, r = -y^l — 4/x, and (72 and ^fy are the SU{2) and t/(l)y couphng constants. 
The effect of such second order annihilation processes on the parameter space of Dirac leptoge- 
nesis is shown in figure 14. 2[ where, for comparison, we show two sets of leptogenesis exclusion 
contours: one representing no second-order processes and one representing annihilation in a 
supersjTiimetric model. It is evident from this graph that second order processes do indeed 
lower the upper exclusion contour, though the effect is not a dramatic one. 

4.2.2 Numerical Results 



Now that we have calculated the requisite rates, we proceed to a numerical evaluation 
of the Boltzmann system. As discussed above, CHDL contains five free parameters, so in 
presenting our results, we will begin by selecting a particular value of 6 small enough to be 
consistent with the neutrino physics constraints in fl3.19l) (we choose 6 = 10~^) and tune 03 
and 63 to the values most advantageous for leptogenesis (in this case 03 = 4.5 and 63 = 2.2). 
Treating the reparameterized coupling strength /, the overall scaling factor in both A and h 
in the parametrization defined in equation fl3.30p . as a free parameter, we display the results 
corresponding to several different choices of in / — 5 parameter space in the right panel of 
figure 14.31 Here, the regions of / — 5 parameter space in which the final value of rj generated 
falls within the WMAP-allowed range given in (II. ip . appear as thin 'ribbons' corresponding 
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Figure 4.2: This figure, taken from [l3], illustrates the effect of second-order processes of the 
form 0101 ij and 0i0i — >■ ij on the exclusion contours from leptogenesis. Con- 
tours are displayed for the case without annihilation (dashed line) and with anni- 
hilation (solid line). 



to each value of M$^. We should reemphasize that S specifies the values of 03 and 63 and 
hence the results displayed are precisely valid only for points lying along the grey, vertical line 
corresponding to the specific value of 6 chosen. To illustrate the effect of changing 6 on the 
shape of the ribbons, we have included a second plot for which 6 = 4.83 x 10~'^ in the left panel 
of figure 14.31 The motivation for selecting this particular value of 6 for contrast will be made 
apparent in section [521 

In figure |473| the effects of the processes detailed in section [4.2. II are apparent, and certainly 
nontrivial. Physically, these effects can be interpreted as follows: increasing the strength of 
the neutrino-sector couplings (here parameterized by /) increases Td, which in turn increases 
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the initial value of B; however, from equation f l4.15p . increasing / also increases the rates 
for the 2 <-i> 2 processes which shuffle lepton number back and forth between Lyis and Lhid- 
Furthermore, it increases the rate for inverse decays. This allows two possibilities for generating 
a realistic final value for B. In the first case, where / is small, the initial baryon number 
produced by (f) and (j) decays is approximately within the range allowed by WMAP, and 2^2 
and inverse decay processes are so slow as to be negligible; this is the "drift-and-decay limit" of 
equation (12.161) . In the second case, where / is large, a surfeit of baryon number will initially be 
produced, but these processes, which occur more rapidly for larger /, subsequently reduce B to 
a phenomenologically acceptable level; this we refer to as the "strong- washout regime" . These 
two possibilities are shown in the two panels of figure 14. 4[ in which the dynamical evolution of 
B, Lyis and other relevant quantities has been plotted, for 6 = 4.83 x 10^^ and M<j, = 10^^ GeV. 
In figure 14. 3[ the two regimes are represented respectively by the lower and upper portions of 
each ribbon — or more properly, by the two points at which these two portions of the ribbon 
intersect the grey line corresponding to the chosen value of 5 in each graph. In the 6 = 10~^ case, 
the strong washout corresponds to / = 8.3 x 10^^ and the drift-and-decay limit corresponds 
to / = 1.5 X 10^^. Alternatively, one can use equation (13.321) to express things in terms of 
(x) once is specified: for = 10^^ GeV, (x) ~ 5 GeV in the strong washout case and 
(x) ~ 50 TeV in the drift-and-decay limit; for = 10^° GeV the two cases converge and 
(X) ~ 1 TeV. 

Having discussed the general effects of the washout processes described above as a group, it is 
also important to address their characteristics relative to one another. Inverse decays dominate 
over 2^2 processes only for a brief period, where 1 z < 50, but during this period they are 
extremely effective in reducing lepton number, and in fact are the primary factor in determining 
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Figure 4.3: Bands in f —5 parameter space for which the final baryon number Bf falls within the 
range permitted by WMAP, for different choices of [38j. Here / parameterizes 
the couplings of the fields (see eq. (13.301) ). The configuration of the left (right) 
panel is the one marked with a dot in the left (right) panel of figure 13.41 in which 
03 = 95 (03 = 4.5). The shaded vertical lines show the constraint on the value of 
5 coming from neutrino mixings and masses {5 = mf./m^ ^ 4.83 x 10~^ in the left 
and S = 10~^ in the right). Note that there are two points that yield a realistic 
value of Bf for a given 6 (the points at which the grey vertical line intersects a 
given ribbon). When / is small enough, the initial baryon number generated is just 
enough to be consistent with WMAP, while the washout effect of inverse decays 
and 2^2 processes is negligible. In this situation the baryon number generated 
is independent of M$ and its the final value is proportional to the CP violating 
parameter e (see eq. (12.161) ). The dark dotted curve shows the band of consistent 
baryon number calculated in this "drift and decay" limit. Each "X" marks the 
point in which Td/H = 1 for each different M$j. At these points, the "strong 
washout" regime starts. The now active washout processes reduce an initial surfeit 
of baryon number (due to a larger /) down to an acceptable level (see figure S3]). 



the final value of rj. For larger z, until they freeze out, the 2^2 interactions dominate and 



further reduce L^s and Lhid (and consequently B). It should be noted, however, that the total 
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B — Ltot number of the universe is manifestly conserved by the Boltzmann equations ( 14.2p - 
(14. 7p (the sum of the rates for the various lepton numbers involved is zero), and since we began 
with B — Ltot = 0, we end up with B — Ltot = in any case, as expected. 




Figure 4.4: These two plots, originally appearing in [38], show the evolution of baryon number 
B for M$j = 10^^ GeV and 6 = irie/m^ = 4.83 x 10~^, in the two different regimes 
that produce a realistic value for the final baryon number of the universe. Bp- For 
the rescaled coupling strength / = ^J\22,h2z = 1-5 x 10~^, as shown in the left panel, 
the effects of 2 -s-^ 2 lepton-number-changing processes are negligible and the final 
baryon-to-entropy ratio is the same as that initially produced by and decays. 
For stronger coupling / = 3.8 x 10~^, as shown in the right panel, baryon number 
is initially overproduced, but 2 2 processes, which are stronger for stronger 
coupling, reduce B to an acceptable level by the time they freeze out. The left and 
right panels correspond respectively to the lower and upper parts of the 'ribbon' in 
figure 14.31 

The upshot of figure 1473) is that Dirac leptogenesis can work even when the constraints from 
chapter UTTl are taken into consideration, but that these constraints have important implications 
for the theory. Combined constraints from baryogenesis and neutrino physics make it extremely 
difficult simultaneously to obtain a realistic neutrino spectrum and obtain the correct baryon 
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number when < 10^° GeV, which in turn requires a reheating temperature Tr > 10^° GeV. 
This means that the constraints from gravitino cosmology discussed in section [3?T1 and especially 
the particularly stringent BBN constraints that arise when m^j2 ^10^ GeV, are of genuine 
concern. Even for heavier heavy gravitino masses 10^ GeV < m^j2 ^ 10^ GeV, we must require 
either that rriLSP is light enough that the naive reheating temperature bound permits a reheating 
temperature above 10^° GeV (see figure |3TT]) . or else that the ratio m3/2/mLSP be large enough 
that LSP annihilations are effective (figure I3.2p . Barring some auxiliary mechanism (such 
as resonant leptogenesis) for increasing baryon number generation, then Dirac leptogenesis 
indeed appears to require heavy gravitinos. However, when 7713/2 is sufficiently large that 
> 10^° GeV is permitted, CHDL succeeds in providing an explanation for the origin of 
the observed baryon asymmetry and neutrino mixings. 

4.3 Satisfying the Flavor Constraints 

The results in the previous section are enough to support the claim that Dirac leptogenesis 
works in theories with a heavy gravitino and sfermion masses heavy enough to evade the flavor 
violation constraints in 13.331 We now expand our analysis to the case where the sparticle 
spectrum is light. Doing so greatly increases the number of parameters in our theory, however: 
the flavor violation rates in (13.441) depend not only on the CHDL model parameters but also on 
the gaugino masses Mi and M2, the Higgs mass parameter fj,, the ratio of Higgs VEVs tan/?, 
and the soft masses for the sleptons in the manner discussed in appendix \M In order to obtain 
precise predictions for the rates, these need to be specified. In our analysis, we choose the 
values Ml = 160 GeV, M2 = 220 GeV, /i = 260 GeV, and tan/3 = 10. As for the slepton soft 
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masses we will assume a common scale = 200 GeV for them. We will also examine what 
effect varying and tan/5 from these chosen values has on BR{fi 67) and BR{t /X7). 
We will assume that the high scale at which soft masses are universal is M = 2 x 10^^ GeV, 
though the results are not particularly sensitive to this choice. We will continue to work in the 
large- (F-^) regime, and since the flavor violation rate is sensitive to {F^, through the left-right 
slepton mixing term fl3.40p (whereas baryon number is not), we need to specify it exphcitly; we 
choose ^/JF^ = W GeV. 

The results of our calculation are displayed in figure In the left panel, we show exclusion 
contours in M$^-(x) space for = 200 GeV. The areas below and to the right of the lower 
contour (the white region) are excluded by the experimental bound from /i — 67. The contours 
associated with r — /i7 are far weaker and do not eliminate any otherwise acceptable region 
of parameter space, and as such we do not display them here. We also include contours 
demarcating the region wherein baryogenesis can succeed. The exclusion contours take the form 
anticipated by equations (13.461) and fl3.49p . and things become problematic when is large 
and when {x) is small, as predicted. As discussed in section [3. 3. 2[ the lower contour represents 
the result of left-right sneutrino mixing contributions Smj^j^ and is primarily controlled by 
(F^). When (F^) is increased, the lower contour in figure 143] is raised, and we find that when 
\/1F]J ^10^ GeV, the entirety of parameter space is excluded. 

In the right panel, we also show how varying the universal scalar mass affects the branching 
ratio for /i 67, which reaches a maximum when is around the weak scale. This is 
to be expected: when is much larger than the weak scale both the slepton mass-squared 
eigenvalues and the flavor-violating terms scale like ml and the sneutrino and charged slepton 
mixing matrices asymptote to a constant value, while the branching ratio is still suppressed by 
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Figure 4.5: Exclusion plots [13] combining constraints from both leptogenesis and flavor vio- 
lation in the process fi — > ej. The left-hand panel shows exclusion contours in 
M$^-(x) space for a universal scalar soft mass = 200 GeV, with tan (3 = 10; the 
right hand panel shows the variation of the branching ratio BR{fi — > 67) with re- 
spect to rus using tan /5 = 3, 10 and 30. In both plots, we have taken Mi = 160 GeV, 
M2 = 220 GeV, and /i = 260 GeV. We have also assigned the x superfield an F- 
term VEV a/ {Fy^) = 10^ GeV. Such a large VEV results in large trilinear couplings 
between Higgs flelds and sneutrinos and therefore induces potentially sizeable mix- 
ings between left-handed and right-handed sneutrinos after electroweak symmetry 
breaking. In each plot, the thick solid contours represent the current experimental 
bound on the branching fraction fl3.33p : the dashed lines represent the expected 
future experimental bound of 1.5 x 10^^^ from MEG. The thin solid contour in the 
left-hand panel delimits the region allowed by leptogenesis constraints. 

the masses running in the loop; as decreases below the weak scale, 5m|^ and 5m\^ go to 
zero and the slepton masses are dominated by flavor diagonal electroweak contributions. We 
also observe that, as in the supersymmetric see-saw case [Ml IH1], the flavor violation rate is 
quite sensitive to tan/5. The lepton-flavor- violating processes depicted in flgure [375] involve a 
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tan /^-dependent chirality flip, and hence are augmented when tan/3 is large. 

The important result we learn from figure 1^5) is that Dirac leptogenesis and the satisfaction 
of flavor violation constraints associated with weak-scale soft masses are indeed compatible — a 
conclusion that would not have been obvious a priori. There is substantial overlap between the 
region of parameter space within which CHDL is capable of reproducing the observed baryon 
asymmetry of the universe and the region in which lepton-sector flavor violation accords with 
experimental bounds. The improvement of these bounds from the next generation of lepton 
flavor- violation experiments will either confirm that some beyond-the-standard-model source 
of flavor violation exists or else rule out a substantial amount of the available parameter space. 
Thus these experiments will serve as an important check on Dirac leptogenesis models with 
slepton masses near the weak scale. We again emphasize that these tensions can be relaxed 
by elevating the universal sfermion mass parameter m^. In theories like split supersymmetry 
where is extremely large, flavor violation constraints can be ignored entirely. 

4.4 The Resonant Escape 

It is clear at this point that the primary model tension in Dirac leptogenesis is a conflict 
between the ceiling gravitino cosmology places on the reheating temperature after inflation and 
the condition that baryogenesis reproduce the WMAP 77. In CHDL, this tension is overcome 
by arranging a large hierarchy among the Yukawa couplings to different sets of $ and $. 
Couplings to the second lightest set, which are crucial to baryogenesis but otherwise have 
little physical effect, are large enough to produce substantial baryon number while couplings 
to the first remain small so as to avoid violating the out-of-equilibrium condition (12.131) . which 
forces Ala and hia down when M^^ is itself required to be small. There is, however, a second 
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possibility for avoiding this bound. It can be seen from the definition of the decay asymmetry 
e in equation (12.71) that if M^^ ^ Mq,^ (in other words, if 5 ~ 1), a resonance condition in is 
obtained. The perturbation theory we have used in calculating equation (12. 7p is good as long 
as the separation of and M^^ is substantially greater than the value of the off-diagonal 
elements Mij in the $ — $ mass mixing matrix induced at the one-loop level: 

~ g*g'{M^, + M$J/,„ (4.21) 

where gi and g'^ represent the appropriate A and summed over the fermion family index, 
and lij is a numerical factor on the order of l/lGvr^ from the loop integral. Thus for small A 
and 5 can be set very close to one and M^^ and M^^ may be very nearly degenerate. This 
allows for the possibihty of resonant leptogenesis, as has been done in the Majorana leptogenesis 
case 1^ [871 [881 El!, which can be invoked to generate a large baryon number in cases where 
the bounds on Tr are more severe. This possibility could potentially allow Dirac leptogenesis 
to work in scenarios such as PeV-scale supersymmetry, where is is otherwise all but ruled out. 



CHAPTER V 
EXTENSIONS OF THE MODEL 



Thus far, we have shown that Dirac leptogenesis is a self-consistent model and that there 
exists a simple version of that model, CHDL, which can simultaneously reproduce both the 
observed baryon asymmetry of the universe and the observed neutrino spectrum without vio- 
lating any additional constraints from astrophysical observations, flavor-violation experiments, 
etc. At this point, we turn to outline some of the theoretical motivations for CHDL and to 
examine some interesting possible extensions of the model, including a potential connection 
between the Dirac leptogenesis superpotential and the origin of the Higgs /i-term, a method for 
canceling U{1)n anomalies, and the possibility for right-handed sneutrino dark matter. 

5.1 Origin of the //-term 

One of the advantages of Dirac leptogenesis is that unlike in Majorana leptogenesis, where 
they are essentially determined by the leptogenesis scale, neutrino masses are set not only by the 
masses M$ . of the heavy fields, but also by the scalar VEV (x) . This leaves us more freedom to 
choose the mass scale of the heavy fields, but can also be thought of as a drawback since we have 
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lost some predictiveness. For this reason, it would be interesting if some of these additional 
intermediate scales could be tied to other intermediate scales that arise in supersymmetric 
models. One interesting possibility involves making use of the U{1)n symmetry to link the 
dynamics of the x fi^ld to the origin of the supersymmetric Higgs mass parameter /x. A simple 
way of doing this is to posit an interaction coupling the up- and down-type Higgs fields with x 
and modifying equation (I2.9p to 

Weff 3 ^^xLpH^N^ + xHuHd, (5.1) 

where is an 0{1) coupling constant. The Higgs superfields are now necessarily charged under 
the new hidden sector symmetry, and therefore the charge assignments given in table 12.11 must 
be revised. Furthermore, the rest of the standard model quark and lepton fields which couple 
to the Higgs fields must also be assigned nontrivial charges under this symmetry. This makes 
the task of arranging anomaly cancelation, arranging for gauge couplings to unify at some high 
scale, etc. vastly more difficult. In this situation, we have ^ = Y^{x) and any observation or 
limits on Higgsino dark matter would directly constrain the VEV (x). (Split) supersymmetry 
might also give us some ideas as to how to relate the heavy fields $ to the SUSY breaking scale. 

5.2 Theoretical Motivations for CHDL 



While the combination of small 6 and antisymmetric h and A which constitutes CHDL may 
yield the correct neutrino phenomenology, we have not offered any explanation of how such a 
situation might arise. Ideally, we would like to have a theoretical framework that would provide 
both the form of the superpotential required for Dirac leptogenesis f l2.ip and the necessary flavor 
structure. We would like to take a moment and indicate one interesting possibility in which 
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simple assumptions enable us to reproduce the needed conditions, while at the same time 
reducing the number of new free parameters that must be introduced. 

In many grand unified theories and models with non-Abelian flavor symmetries, the Yukawa 
matrices (and in general the flavor interactions) can be symmetric, antisymmetric, or both (see 
for example [71], [73l [71] and references therein). Operating in this paradigm, let us assume 
that the SM left-handed leptons have the same flavor charge as the heavy fields $, and the SM 
right-handed charged leptons have same flavor charge as the fields $. Upon breaking of the 
flavor symmetry, let us assume that the charged lepton Yukawa matrix is symmetric due to a 
symmetric flavon VEV configuration (Sajj)- The corresponding effective superpotential is 

WdYi ^^H^L^ep + F$ {S^f,)^^<l>f,+h.c. (5.2) 
IVI p 

where Yi and Yq> are dimensionless couplings, and a and are flavor indices. The charged 
lepton Yukawa matrix (yja/j and the mass matrix (M$)q,/3 of the heavy fields $ would be both 
symmetric and proportional 

{M^)^p = Mp^ (yj„^ (5.3) 

This specific structure predicts exactly the mass spectrum for the fields in terms of the 
flavor scale Mp, which is at the origin of the intermediate scale required for successful thermal 
leptogenesis. 

Mf Mp ^ Mp 
= rrie , iW$2 = M^,, = 711^ (5.4) 

V V V 

If the flavor scale is of the order of some GUT scale Mp ~ 10^^ GeV, we would then expect 
Mq,^ ~ 10"*^^ GeV, with 5 = M<^^/Mq,^ = nie/m^ = 4.83 x 10"'^ being a small parameter. 
Because this ratio is quite small and e is approximately proportional to 6 for small 6, one 
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might worry that this might doom baryogenesis. However, as indicated in the left panel of 
figure [374l a™'^^(me/m^) ^ 95, so that a hierarchy between couplings to different sets of $ and 
$ is permitted, and the result is that e is only suppressed by an 0{1) numerical factor. This is 
corroborated by the results of our numerical analysis, displayed in the left panel of figure 14.31 
and explains why we chose to examine the 5 = 4.83 x 10~^ case. 

Another related possibility that can help reduce the proliferation of scales (M$^, Mp, etc.) 
in the model is to link the supersymmetry-breaking scale Msusy to other physical scales in the 
theory by coupling the spurion field X responsible for supersymmetry breaking to the and 
$j fields. If the spurion field is charged under the same symmetry (e.g. the U{1)n introduced in 
section I^TTI) . responsible for dictating the form of the superpotential in fl2.ll) . one could arrange 
that 

W^^{Yi HaL^ep + X l>,$^) +h.c. (5.5) 
Ivi p 

where again Yj and 1$ are dimensionless couplings. Supersymmetry breaking effects can provide 
the field X with a VEV {Ax) ~ (m3/2Mp^)^^^ and upon flavor symmetry breaking we could 
get the effective superpotential 

Weff D (l/Ja/3 {HdL^ep + ^m^/2Mpi <l„$^) +h.c. (5.6) 

where we have assumed that the original constants Yi ~ Y$. In this situation, for example 
with m3/2 ~ 10^*^ GeV (a value sufficiently high as to avoid any model tensions associated with 
gravitino cosmology), the mass of the lightest $ field would be ~ 10^^ GeV. 

Now, let us also assume that the fiavon VEV {Sai3) is symmetric and that the coupling 
becomes antisymmetric upon flavor breaking, i.e. the superpotential can be written as 

W D {y,y^^H,L^ep + {M^Y^^^.^p + KpN^^pH^ + C/'^^l-.L^x (5.7) 
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where we have Aq/j = {A^^) and h'^^ = (^q^), with the flavon acquiring an antisymmetric 
VEV configuration and with the VEV configuration of A^ being arbitrary in fiavor space. 
As alluded to above, this assumption is well-motivated from GUT considerations. While the 
fiavor structure of A^ is not important for our present purposes, in a specific model of fiavor 
it would likely end up being some linear combination of symmetric and antisymmetric VEV 
configurations. 

Let us see what the implications of these ingredients are for our model. In general, the 
charged lepton Yukawa matrix can be diagonalized by a biunitary transformation of the 

type 

U^yy=yl'^^'^s)^ (5.8) 
but when yi is symmetric, this biunitary transformation takes the simpler form 

U'^y,U = yl'''^\ (5.9) 

If (M$) oc (yj, as in the setup described above (see equation (15.31) ). then (M$) will be diago- 
nalized by the same transformation. When the mass matrices for the charged leptons and the 
— system are simultaneously diagonalized (that is, when we go to the charged lepton basis), 
A and h transform as 

A' = f/^Af/ h' = U^h'^^'^U. (5.10) 

Transformations of this type preserve the antisymmetry of h (and A if also antisymmetric); 
thus in the charged lepton basis the matrix h remains antisymmetric and (M$) is real and 
diagonal; this reproduces the neutrino mass matrix structure given in (13.291) . in which the 
diagonal elements of h are zero and it is the smallness of 6 that is responsible for the large 
mixing angles observed in the lepton mixing matrix Umns- 
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From a purely structural point of view, it will be noted that any matrix of the form 

(M$) = + B/3X3, (5.11) 

where A and B are arbitrary constants and /3X3 it the 3x3 identity matrix, can be diagonalized 
along with yi. Thus (M$) oc yi is not required for the transformation rules in (15.101) to hold. It 
follows from this that M$ may receive arbitrary diagonal contributions and 6 = Mi^^/Mq>^ = 
me/rrif^ is not required: 6 can in principle take any value (as long as it is consistent with the 
observational bounds on the neutrino spectrum). This gives us a great deal more freedom 
to adjust model parameters without sacrificing the theoretical motivations we have developed 
here. 

5.3 Sneutrino Dark Matter 

5.3.1 Dark Matter in Dirac Leptogenesis 

In order to be a viable dark matter candidate, the LSP must be stable (which follows 
automatically from i?-parity conservation) and interact only weakly with the standard model 
fields. The most obvious choice is the lightest neutralino, A^^i, and this choice is perfectly 
compatible with Dirac leptogenesis. The choice of model parameters outlined in section 14. 3[ 
for example, yields a stable, predominately bino LSP with a mass of around 160 GeV. In the 
context of PeV-scale supersymmetry or split supersymmetry, little changes except for that the 
LSP mass may be quite large, potentially facilitating the detection of energetic photons from 
cold dark matter annihilation at the galactic center by the next generation of Cherenkhov 
telescopes [201 ESI EDI EI] . We explore this possibility further in chapter |VT1 
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Dirac leptogenesis does, however, offer a new potential candidate for LSP: the hghtest right- 
handed sneutrino Unhke in Majorana leptogenesis scenarios, where right-handed sneutrino 
masses receive a supersymmetry-respecting contribution on the order of the leptogenesis scale, 
here they acquire mass only through soft terms, and since the i7r are still singlets under the 
Standard Model S'f/(3) x SU{2) x f/(l)y gauge group, it is certainly possible that their soft 
masses are significantly smaller than — and perhaps of a different origin entirely from — those 
of the SM squarks and sleptons. Obtaining the observed value of ^cdm (13.41) with a pure 
right-handed sneutrino LSP is somewhat difficult, however, owing to the fact that the particle's 
interactions are so weak that annihilations are insufficient in reducing its relic density to the 
appropriate level. One solution to this problem is to make the LSP a mixture of left- and 
right-handed sneutrinos [911 ES] with mixing angle 6. The SU (2) gauge interactions of the ul 
then permit the particles to annihilate efficiently enough to reduce Qlsp to an acceptable level, 
but (since the annihilation rate is suppressed by sin^ 9) not too efficiently. Here we will take 
another approach: we will gauge the global U{1)n symmetry we introduced in order to forbid 
Majorana masses for neutrinos and allow up^^ to annihilate down to an acceptable level via their 
interactions with [/(I) at gauge bosons and gauginos. 

Since neutralino dark matter still works in the context of Dirac leptogenesis, investigat- 
ing the viability of another, more exotic CDM candidate might appear a somewhat frivolous 
endeavor. However, there is a compelling reason why sneutrino dark matter is of particular 
interest. It is a puzzling but well-documented fact about our universe that the dark matter 
abundance flcDM (13. 4p and the abundance of baryonic matter Qb (11. 2p are roughly the same 
order, and in the standard picture of neutralino dark matter there is no known reason why 
this should be the case. However, in situations where an asjnnmetry develops in some globally 
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conserved quantum number under which the LSP is charged [93], it is possible to hnk these 
two quantities. Consider for a moment an effective theory containing both visible and hidden 
sectors, coupled only by dynamics at some high scale whose precise nature is unimportant. 
As the universe evolves down from the scale of the unifying dynamics, an asymmetry in this 
conserved quantum number develops between the two sectors. Suppose the hidden sector con- 
tains a number of light fields charged under the relevant global symmetry, each with charge 
Qi and number density rii, and their antiparticles, with charge —Qi and number density n,. 
The respective particle asymmetries A^s and Ahid in the visible and hidden sectors, normalized 
with respect to the entropy density of the universe s, are then 

Ahid = ^ = ~Ams. (5.12) 

s 

If the hidden sector contains a massive, stable (or at least extraordinarily long-lived) particle 
a with mass m^, this particle becomes a potential dark matter candidate. Let us assume that 
it also carries a sizeable fraction of this asymmetry Aa = cxAhid at the time of the freeze-out 
epoch and consider the limit where annihilation processes are efficient and conserve Aa- In this 
case, a and its conjugate a* will annihilate rapidly with each other until only one or the other 
remains (even if the annihilation rate is quite large, a and a* cannot annihilate further via such 
processes once either a or a* is entirely depleted). The particles left over, having no conjugates 
with which to annihilate, become the dark matter, with a relic abundance 

rtcDM = "^"^"^"^^^ = 2.236 X 10^°ma\Aa\ (5.13) 

Pcrit,0 

proportional to the asymmetry. An asymmetry on the order \Aa\ ~ 10^^^ — 10^^° (depending 
on nia) will thus replicate adequately the observed dark matter abundance. Since Majorana 
neutralinos are their own antiparticles, they must have vanishing U{1) charge, and thus cannot 
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be used to link Qb and ^cdm in this manner. In contrast, when effective A-terms mixing 
left- and right-handed sneutrinos (such as those engendered by (F^)) are small, a right-handed 
sneutrino LSP in Dirac leptogenesis fits the bill perfectly: it is stable; it resides in a hidden 
sector essentially decoupled from the visible sector fields until well after the freeze-out epoch; 
and it is charged under the globally conserved quantum number B — L, under which the visible 
sector fields are also charged. The B — L asymmetry between the hidden and visible sectors 
created during the decay of the heavy $i and $i superfields will also determine the relic density 
of uri, and consequently the mechanism responsible for baryogenesis will also be responsible 
for generating VLcdm- 

The simplified scenario resulting in fl5.13l) is somewhat difficult to realize in practice, how- 
ever, due the fact that it simultaneously requires the net rate T{a*a X) (where X represents 
some unspecified final state comprising light, hidden-sector particles) for A^-conserving anni- 
hilation processes to be efficient in reducing the dark matter abundance to a cosmologically 
acceptable level and the rate T{aa X) for A^- violating processe^jto be negligibly small. The 
correct hierarchy of rates can be achieved in certain circumstances — in the presence of large 
s-channel contributions to r{a*a — > X) that have no t-channel equivalent, for example — and 
we will discuss these possibilities when we solve the Boltzmann equations for a right-handed 
sneutrino LSP. Before we do this, however, we turn briefly to address the modifications neces- 
sary for gauging the U{1)n symmetry, as this is the simplest and most convenient method of 
making L-conserving sneutrino annihilations efficient. 



""^We assume that CP-violation can be neglected and that the rates for the conjugate processes aa — > X and 
a*a* X* are equivalent, hence Taa = ^a*a*- 
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5.3.2 Extending the Hidden Sector 



In order for our theory to be consistent, we must require that our new U{l)i\f symmetry 
is non-anomalous. With the field content listed in table 12.11 we find that there are no mixed 
SU{2)IU{1)n, U{1)I.U{1)n or U{1)yU{1)1 anomalies, but that 

{grav.fU{l)N: 3(+l) + 3(-l) + 3(+l) + (-1) = 2 (5.14) 
U{1)%: 3(-l)3 + 3(-l)=^ + 3(+l)3+ (-1)3 = 2 (5.15) 



are non-vanishing. It has been shown [91] that in the case where the only non-vanishing 
anomalies for a hidden sector U{1) are U{lY and {grav.)'^U{l) (i.e. mixed anomalies vanish), it 
is possible to cancel those anomalies by introducing new fields charged only under that U{1). 
We therefore look for an appropriate expansion of the field content in the hidden sector that will 
cancel the above anomalies. The simplest way to cancel these anomalies would be to introduce 
two additional fields with the same charges as but mixing between the fermionic components 
of these fields will ensure that at least one linear combination of them will be massless. Since 
these fields are necessarily odd under i?-parity, this choice is problematic for dark matter and 
therefore must be excluded. The next simplest optioiu requires four new fields (a = 1, 2), yj, 
and p, the charge assignments for which are recorded in table 15.11 (the reason for the i?-parity 
assignments therein will become evident in a moment). Adding these fields permits us to write 
down two new superpotential terms 

W D pXP + CaX^a^, (5.16) 

where /i is a supersymmetric mass term and (a (c^ = 1, 2) are dimensionless couplings. 



^Here we have assumed that aU new superfields have U{1)n charges (not necessarily integral) Qi < 20. 
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Field 


U{1)l 


f/(l)7V 


SU(2) 


U{l)y 


Pm 







-4 


1 





+ 1 


V 





5 


1 





+1 


p 





1 


1 





+1 



Table 5.1: This table shows the charge assignments for minimal set of additional superfields 
necessary to cancel the U{1)% and {grav.YU{l)isi anomalies of the fields in table l2.1l 

Since our aim is to arrange for P/j^ to be the LSP, we must make sure that none of the new 
i?-parity-odd fields in the hidden sector ends up massless and that the scalar component of x 
can still acquire a VEV on the appropriate order to yield a realistic set of neutrino masses. 

The scalar potential in the hidden sector, which contains both an F-term contribution 
from (15.161) 

^ \PP + Ca^a^l^ + l/iXr + ICaXU^ + (101' + IC2niXV^|' 
^ IP'W + l/^rixl' + {pXa^a^P* + h.C.) 

+ ICam^l' + ICaHeaXp + (IClP + IC2niX¥'P (5.17) 



and a D-term contribution from the new U{1)n gauge interaction given by 



Vn = '-{\n\' 



2N2 



Xp + |pP-4|6p-4|6P + 5|^n 



(5.1^ 



both of which serve to stabilize the fields. In addition, a contribution Vs will arise from soft 
supersymmetry-breaking Lagrangian. As for the fermions (which, as they are singlets under 
the Standard Model SU{3) x SU{2) x U{1)y, will be referred to as "neutralinos" from this 
point forward), their mass matrix, in the {Z', p, ^i, basis, is 



M 



( M' -V2g ix) V2g (p) -AV2g (6) -4V2g (6) 5V2g (if) \ 

M Ci{^) C2{v) Ci(ei) + C2(6) 



Ci(x) 

C2(x) 

V • • . . . / 



(5.19) 
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As long as the scalar components of x ^"^^ receive nonzero VEVs, all the physical mass 
eigenstates resulting from this matrix will be massive, as desired. This can be engineered by 
introducing tachyonic soft squared-masses for these fields in V5. Assuming the other scalars 
^2, V^, and p receive large, positive squared-masses the effective (tree- level) potential for the 
remaining fields becomes 

v = {\p?+ mi)\x? + (i/ir + mim^ + icipieirixr + y (ixr + m'Y- (5.20) 

It should be noted that since the U{1)n charges of x and are of the same sign, there are no 
problems with Z^-flat directions. 

Minimizing the potential in equation fl5.20p . we find that in terms of the quantities 

ml ^ -i\fx\' + ml), ml ^ -(|/i|2 + m|j, (5.21) 

the expectation values for and x are given by 

(4 ml — ) 



ixr> = 4 

9 



3|C: 



|2 



~ 2 vy ' J / A ~ 2 ~ 2 



(5.22) 



3|Ci 

Equat ion 15 . 2 1 1 implies that in order to avoid having to introduce fine-tunings among the hidden 
sector mass parameters fi, m^, and m^^, these parameters ought to be of roughly the same 
magnitude as the physical scales (x) and (^i) (if fl were small or vanishing, some of the hidden- 
sector neutralinos would be too light). In order for both VEVs to be positive, we must have 

< ^ < 4, 5.23 

a J 

where a = (f?^ + ICiH/SlCiP ^ 1/3- Both of these inequalities must be satisfied in order for 
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and X to receive VEVslfl As for Z' vector boson associated with the gauged U{1)n symmetry, 
when the symmetry is broken by the VEVs of and it acquires a mass 

M|, = 2/((x)' + 16(6)'). (5.24) 

Thus all fields in the hidden sector acquire masses on the order of the scales (x) and {^a) 
or larger, and the lightest right-handed sneutrino remains a viable LSP candidate. In fact, it 
should be pointed out that gauging U{1)n introduces new D-teim. contributions to the physical 
mass of the right-handed sneutrino which serve to drive down ml^. The total result, including 
soft terms, but assuming that left-right mixing through effective A-terms is negligible, is 

< = - / {{\x?) + 4 im + 4 (i6r> - (ipr> - 5 , (5.25) 

and since m?^ ^^^^ is essentially a free parameter, it can be tuned freely so that the mass of 
the lightest right-handed sneutrino is small, as desired. It should also be pointed out that this 
setup results in a VEV for the scalar field x without necessarily producing an F-term VEV 
(x), and hence the right-handed fields are indeed part of the hidden sector, as is required for 
relating Qcdm and Qb in the manner discussed above. 

5.3.3 Evolution of the Dark Matter Abundance 

In order to ascertain the dark matter abundance in any nontrivial model, it is necessary 
to solve the coupled system of Boltzmann equations that govern the evolution of ur^ and z/|j.^ 

during the freeze-out epoch. The two-to-two processes which yield the leading contribution to 

^Equation (|5.23|) also indicates why the matter parity assignments given in table 15.11 were chosen as they 
were. If and tp were taken to be odd under matter parity, the coupling constant in (|5.16|) would vanish, 
sending a ^ oo. 
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dark matter annihilation include both processes of the form vr^v*^^ X, which contribute to 
washout of the dark matter density but do not alter Lj>^ (which is equivalent to {B — L)yj^, 
since electroweak sphalerons do not affect the fields in the hidden sector and hence B and L 
are separately conserved there), and those of the form i'r^i'r^ X, which contribute to the 
depletion of both Qcdm and {B — L)i)^. Let us call the net rate for the former processes Tj^o* 
and the rate for the latter rj>,>. The Boltzmann equations can be written either in terms of 
t'R^ and z/^^ or, equivalently, for the dark matter abundance Ydm = + Yj^*^ and the lepton 
number asymmetry La^. We find these equations to be 



dz H (1711)^)3 



dYnM z 



dz 2H{mj;,^)s 



Taa {l,,Ydm - ^-^{U., - L,,)YS'm) (5.26) 
r(+) (r-/ - Y^,,) + r(-) L^J , (5.27) 



where z = rriy^/T, and H{m^^) is the Hubble parameter at scale T = m^^j. The effective rates 

T,, = ^,,/Y^i and = T,,, ±2T,, = A{^,,. ±^,,)/Y^i (5.28) 

are determined from the reaction densities 'yj^c and 7i>i>*, defined as in (14.161) . which are in 
turn determined from the reduced cross sections 6"(z/z>) and (j(z/z/*), a calculation of which is 
included in appendix [Dl Given in terms of r = ^1 — 4m?^/s, = s/m^^, a;„ = s/m?^, and 
A = {2x^^ — 2x^^ — 1), the results are 

+!!l(^)(2Ar+^J^ln(^Y\. (5.30) 

87r Vs-M| M 2 V^-^, 
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We note that the equations ( 15.26l - [5l27l) are actually quite general and can be applied in 
a variety of scenarios in which an asymmetry in some globally conserved quantum number 
between visible and hidden sectors (or between a pair of distinct hidden sectors) develops 
during the evolution of the universe and the dark matter particle is located in the hidden sector 
and charged under that quantum numbero Any model-dependent aspects of the Boltzmann 
evolution enter only in the relevant annihilation rates, in the mass of the dark matter particle, 
and in the input value for L^i^. 

Before we move on to address the Boltzmann system f l5.26l - [5l27j) numerically, we note that 
in certain limiting cases it reduces to a particularly simple form. One of these cases is the 
trivial case, where Lj^^ = and we are left with the single equation 



{yZs'-yE,,) (5.31) 



dt 

whose solution is the standard one for a dark matter abundance. Perhaps more interesting is 
the case where ~ and L,>^ ^ 0. In this case, Lj^^ is fixed and the equation for the dark 
matter abundance becomes 

dYDM 1 



— r 



{YdV - yh.) + Ll^] ■ (5.32) 



dt 2s 

This represents the case described in [93]. As mentioned above, this situation is somewhat 
difficult to engineer for a right-handed sneutrino in Dirac leptogenesis, but in the next section 
we will provide one or two possible schemes in which it can be realized. 

5.3.4 Numerical Results and Discussion 



■^In special cases, such as that of a nearly-degenerate LSP-NLSP pair in supersymmetric models, the Boltz- 
mann equations must be modified |95| . but we do not consider such exceptions here. 
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We are now ready to solve the Boltzmann system given in ( ]5.26p - ( 15.27P for right-handed 
sneutrino dark matter in supersymmetric Dirac leptogenesis. Obtaining the correct dark matter 
abundance in this scenario is not difficult, provided the f/(l)iv gauge coupling constant is on the 
appropriate order, but generating Qcdm via the mechanism described in equation (15.131) and 
thereby forging a link between this quantity the observed baryon asymmetry of the universe is 
not. The reason for this is that the Li>-conserving and L^^-violating rates determined respec- 
tively by fl5.29p and fl5.30p are naturally on the same order unless a hierarchy exists among the 
masses of the virtual particles exchanged. The Z' mass given in f l5.24p depends primarily on the 
hidden-sector gauge coupling g, the soft squared masses and m|^, and the supersymmetric 
mass parameter /2 (which ought to be of roughly the same order as the aforementioned soft 
masses if we wish to avoid fine-tuning); while the neutralino mass and U{1)n gaugino content 
depend primarily on M'. Thus unless a hierarchy between scalar and gaugino masses emerges 
naturally from supersymmetry breaking in which gaugino masses are far heavier than scalar 
masses, (5"(z/z/) and (j(z/z/*) will be of roughly the same order. 

This leads us to choose between the two less-than-ideal outcomes depicted in figure 15.11 
Here, we show the evolution of Ydm and Lj^^ in two representative regimes: one in which the 
U{1)n gauge coupling g is comparatively strong {g = 9.5 x 10^^), and one in which it is weaker 
{g = 5 X 10^^). In each case, we have maintained the relations g{x) = 200 GeV, g{C,i) = 
100 GeV, fi = 500 GeV, and Ci/g = (2/9 = 1 and set M' = 200 GeV. By examining figure SSI 
one can verify that these parameter choices are both compatible with successful leptogenesis, 
given an appropriate choice of M$^. The mass of the lightest right-handed sneutrino is taken 
to be m^)^ = 150 GeV, which makes it the LSP. In the strong coupling case presented in the left 
panel, r(z>z/) and r(z/z/*) are large enough to reduce the sneutrino abundance to a level in accord 
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z=m|_gp/T z=mLsp/T 



Figure 5.1: Here, we illustrate the two representative cases for sneutrino dark matter in Dirac 
neutrinogenesis. In each case, we have set M' = 200 GeV and m^^ = 150 GeV 
and maintained the relations g{x) = 200 GeV, g{C,i) = 100 GeV, ft = 500 GeV, 
and (i/g = (2/ 9 = 1- The left panel depicts the situation where the U{1)n gauge 
coupling is reasonably strong {g = 9.5 x 10~2). As the universe evolves, T decreases 
and 2 — > 00. Here, the number-density-to-entropy-density ratio Ydm of a right- 
handed sneutrino LSP reproduces the WMAP prediction (the horizontal dashed 
line), but the lepton number stored in right-handed sneutrinos is washed out 
almost entirely and all connection between the dark matter abundance and the 
observed baryon number of the universe is lost. The right panel shows the situation 
where the gauge coupling is small g = 5 x 10~^. Here, washout is insignificant and 
L{^^ persists essentially unchanged from its initial value until present time, but 
sneutrino annihilations are inefficient in reducing the dark matter abundance to an 
acceptable level. The equilibrium right-handed sneutrino abundance Y^\j is also 
shown in each panel for reference. It should be emphasized that in the left panel, 
the correct dark matter relic density is obtained, and that while ^Icdm and Qf, are 
not tied to one another, there is no problem for cosmology, whereas the situation 
depicted in the right panel is cosmologically unacceptable, since dark matter is 
overproduced. 

with the WMAP bound, but r(z/i>) is so large that Lj,^ is washed out rapidly and essentially 
completely. Here, the model is successful in yielding the correct dark matter abundance, but 
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all connection between L^^ (and therefore B) and flcDM is lost, and the evolution of Ydm 
reduces to the standard situation described in equation f l5.3ip . In the right panel, for weak 
gauge coupling, things go even more awry: T^uu) is small enough to prevent the washout of 
L^^, but since r(i>z/*) is likewise small, sneutrino annihilations are inefficient in reducing the 
right-handed sneutrino abundance to an acceptable level — in essence, we recover (literally, in 
the g —>■ limit) the situation we had before gauging the U{1)n symmetry. The message here 
is that the right-handed sneutrino in Dirac leptogenesis is a viable dark matter candidate, but 
that a relationship between Qcdm and fib does not emerge naturally from the Dirac leptogenesis 
framework, as least in the simple model presented here. 

Arranging a hierarchy between r,>i>. and r,>i>, though difficult, is not altogether impossible, 
however. As mentioned previously, one way this could occur is if a hierarchy between gaugino 
and scalar masses, in which the former were much higher than the latter, emerged naturally from 
the supersymmetry-breaking mechanism, but this is not easy to achieve. The most promising 
model would appear to be gaugino mediation [SSI EZl ES, EH], in which the supersymmetry- 
braking sector is confined to a "source" brane parallel to our own, visible brane in a five- 
dimensional bulk. The matter fields of the theory (the MSSM fields, right-handed neutrino 
superfields, etc.) are confined to the visible brane, while gauge supermultiplets are allowed to 
propagate in the bulk. As a result, one can arrange for a loop-splitting between gaugino and 
scalar soft masses, but it is difficult to arrange for this splitting to be as pronounced as we 
require. We thus turn to alternative methods for engineering a scale-separation between Fj^ijc 
and r£;j>. 

A second way of arranging a hierarchy between these rates is to further expand the field 
content of the theory and introduce additional light scalars into which the ur^ can decay. This 
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can be arranged by introducing pairs of vector-like superfields Ui and uji that are even under 
matter parity (so ur^ is still the LSP) and charged only under U{1)n- Charges Qi are assigned 
them so that the form of the Dirac leptogenesis superpotential (12.11) and its hidden sector 
extension (15.161) are unchangedp save for the addition of mass terms 

W D M^^UiuJi. (5.33) 

Since the new superfields come in vector-like pairs, they will not spoil anomaly cancelation. 
The soft Lagrangian will also contain new soft mass terms for the scalar fields in uJi and cJj, 
and a partial cancelation between a tachyonic soft squared mass and the contribution M^. from 
the superpotential can arrange for the scalar components of Ui and Ui (also denoted Ui and cJj) 
to remain light compared to ur-^ , while their fermionic superpartner a)j and aJj acquire masses 
> fTT-UR- In this case, an additional s-channel contribution to the right-handed sneutrino 
annihilation cross-section via the diagrams in figure 15.21 will arise, which has no t-channel 
equivalent. The contribution to o"(i>z>*) from each pair of uJi and oJj is 

g'Qlf s_ 



where r is defined as above equation (15.29p . and if the number of additional scalars is large 
(or the charges assigned them are large), the L^^^-conserving rate can be increased without 
simultaneously increasing the L,y^-violating one. 

While this idea seems promising, it turns out to be difficult to realize in practice. The 
reason is that the reduced cross-section (15.341) for sneutrino annihilations, like the usual expres- 
sions (I5.29P and (15.301) . is proportional to g'^, and g must be kept small in order to prevent L,>^ 



^This is trivial and can be easily accomplished, for example, by giving assigning fractional charges to all the 
LOj and cJ,- . 
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Figure 5.2: The additional contribution to alui)*) that arises in the right-handed sneutrino 
dark matter scenario in the presence of additional light scalars Ui and ZJj charged 
under the U{1)n gauge group. 

from being washed out. This substantially hinders the effectiveness of this new contribution at 
reducing Ydm- one must either add a very large number of vector-like superfield pairs to 
the theory — or else make the charges Qi of the new superfields inordinately large — in order to 



achieve the desired hierarchy between F^i) and F/^i^*. In figure [573| we show the effect on Ydm 
of adding light scalars to the theory. Here, we have chosen the (quite large) value |Q| = 20 for 
the magnitude of the charge assigned to all Ui — ZJi pairs and set (7 = 1.5 x 10~^ (the input 
values for all other relevant model parameters are the same as in figure ISTTl) . To compensate 
for the smallness of the ~ 10~^^ prefactor in fl5.34p . an additional ~ 10^ field pairs are 
needed. This is highly unmotivated, to say nothing of the havoc such fields would wreak on the 
RGE running of the U{1)n couphng in terms of the prospects for gauge coupling unification, 
etc., and therefore not a particularly compelling solution to the rate correlation problem. 

The possibility that Dirac leptogenesis presents for right-handed sneutrino dark matter is 



indeed an interesting one, and as the left panel of figure 15711 one that can reproduce the observed 
value of VLqdm- The disappointment is that it it a challenge to tie the size of the observed dark 
matter matter abundance to the baryon number of the universe through their mutual connection 
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Figure 5.3: The effect on the dark matter number-density-to-entropy-density ratio Yd a/ of 
adding a number of vector-hke pairs of hght scalars to the theory. Here, all 
such scalars have U{l)]sf have charge magnitude \Q\ = 20, and results are shown 
for A*"^ = 10, 10^, and 10^, as well as for the unmodified case A^^ = 0. As in the 
right panel of figure 15.11 M' = 200 GeV and m,>^ = 150 GeV and maintained the 
relations g{x) = 200 GeV, ^(^i) = 100 GeV, jl = 500 GeV, (i/g = C2/9 = 1, and 
g = 1.5 X 10^'^. The addition of light scalars does not affect the lepton number 
stored in right-handed sneutrinos appreciably for > 10. Because g is small, an 
inordinately large number of light fields are required in order to achieve a situation 
in which dictates the size of Yd a/- When A^^ is large enough, L^,^ creates an 
effective "floor" for the dark matter abundance (corresponding to a dark matter 
abundance consisting only of and none of its antiparticle i^r*). 



to L,y^ in a way where the former is driven automatically driven by the latter to the correct 
value. It should again be emphasized that the link between L^^ the baryon number B of the 
universe is forged at the leptogenesis scale > 10^° GeV. The two asymmetries do not 
communicate with one another during the CDM freeze-out epoch, and the the rate Fp,) does 
not alter the total lepton number Lhi^i stored in the hidden sector, but rather transfers the 
portion of that asymmetry L,^^ stored in right-handed sneutrinos to right-handed neutrinos 
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(or perhaps light scalars). Thus the dynamics associated with Tj^o have no effect on B, and 
baryogenesis proceeds unhindered no matter to what degree L^^^ is washed out. Thus sneutrino 
dark matter is perfectly compatible with baryogenesis in Dirac leptogenesis scenarios. However, 
it is generally not feasible to forge a connection between Qcdm and the rehc abundance Qb of 
baryonic matter through a mutual connection to Lp: a correlation between the the rates F^jj and 
ri>j> for Lj>-conserving and L^^-violating processes makes it difficult simultaneously to preserve 
while reducing to an acceptable level. 



CHAPTER VI 



INDIRECT DARK MATTER DETECTION 

6.1 Detecting Heavy Neutralino Annihilations 

In section I5.3[ it was mentioned that Dirac leptogenesis is compatible with neutrahno dark 
matter. Furthermore, it strongly prefers a heavy gravitino and hence is quite naturally com- 
patible with split supersymmetry, which is associated with a particular, characteristic type of 
neutralino dark matter. In theories where anomaly mediation is responsible for the generation 
of gaugino masses Mi, M2, and M3 (as is the case in split supersymmetry), the neutralino spec- 
trum is dictated by the /^-functions of the Standard Model gauge groups by equation fl3.14p . 
In the MSSM, to lowest order, the relationship is 

M3 ~ 3Mi ~ 9M2, (6.1) 

which implies that if a neutralino is in fact the LSP, it will be predominately either Wino or 
(depending on the value of /x), Higgsino. We now turn to investigate the possibility for detecting 
such an LSP experimentally. 

In section 13.11 we discussed the relic abundance of such an LSP, which had both a thermal 
component ^Isp given by equation (13.51) or (13. 6p for a majority Wino or Higgsino LSP, re- 
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spectively, and in theories with a heavy gravitino, an additional, nonthermal component fi^sp 
given by equations (13. 9p - (13. lip , resulting from gravitino decay. As is evident from figure [3T 
the thermal component places an upper bound on misPi which turns out in each case to be [61 



rriLSP < 2.5 TeV (for Wino) (6.2) 
rriLSP < 1-2 TeV (for Higgsino). (6.3) 

The direct detection of a heavy Wino or Higgsino LSP with a mass misp Mw GeV is 
all but precluded at present and planned facilities [59] for all but the most unnaturally peaked 
halo models, and the cross-section for the annihilation of such an LSP into positrons is too 
small to yield a detectable signal. The only promising detection method available is to search 
for high-energy photons resulting from LSP annihilation processes in the galactic halo, and in 
particular NiNi 77 and NiNi —>■ 7Z. Expressions for the cross-sections for these processes 
have been computed [lUUj . and in the limit where misp is much larger than the weak scale, 
av{NiNi — > 77) and av{NiNi Z7) take the asymptotic forms 

av(WW 77) ~ 4.0 X 10"^^ cm^s"^ (6.4) 

av{WW Z-f) ~ 9.0 X 10"^^ cm^s-\ (6.5) 

for a Wino LSP, and 

av{HH 77) ^ 9.0 X 10"^^ cm^s^^ (6.6) 

(rv{HH Z-f) ^ 2.0 x 10"^^ cm^s'^ (6.7) 

for a Higgsino LSP. The reason these processes are of particular interest is that the photons 
produced thereby are effectively monoenergetic, with respective energies 

1 - J ■ (6-8) 
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The detection of either signal would be compelling evidence for heavy Wino or Higgsino dark 
matter, since effectively monoenergetic signals at TeV-scale energies tend to originate in par- 
ticle physics processes (rather than astrophysical ones); taken together, they would serve as a 
distinctive signal — were it possible to resolve them. Our aim here is to analyze the prospects 
for the detection of this signal at present and future 7-ray telescopes. 

Calculation of the observed integral flux $ of 7-rays from the galactic center is complicated 
somewhat by our ignorance of the precise density distribution profile plip, s) of dark matter in 
the galactic halo (expressed here as a function of line-of-sight distance s and the angle ip away 
from the galactic center). The density profile enters into the integral flux (usually expressed in 
cm~^s~^) through the expression 



where N^{X) is the number of photons in the final state X, and the integral is evaluated along 
the line of sight. It should be noted that p(^/', s) is proportional to the LSP mass and the 
asymptotic expressions for av{NiNi X) given in 06.41 - 16. 71) are independent of misp, so 
$ is actually insensitive to misp when misp ^ M]y. Since the density integral is essentially 
independent of the particle physics, it is common practice to abstract it by defining the quantity 



While for an ideal detector, one would be interested only in J(0), the value for the line passing 
directly through the galactic center, a real detector will receive incoming photons originating 
over some finite slice of solid angle. The relevant quantity is angular acceptance AQ, the solid 
angle from which the detector actively receives light, which for a given detector may be adjusted 
as desired over a window ranging from a detector's angular resolution to its field of view. Thus 




(6.9) 




(6.10) 
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a P 'J R 


Isothermal profiles 
NFW 

Moore et al. 


2.0 2.0 3.5 
1.0 3.0 1.0 20.0 
1.5 3.5 1.5 28.0 



Table 6.1: The defining parameters a, (3, 7 and R (see equation 16. lip for the halo models we 
examine. R is given in kpc. 

instead of J(0), the quantity of interest is {J{ip))An, the average of J{ip) over AQ. 

-i9 



— 1 1 — n 1 1 — n 1 1 — n 1 1 — r 



T — TT 1 1 — n 1 1 — n r- 



C. ISO. 
Burkert 
S. ISO 
NFW 
Moore 




Figure 6.1: {J{ip)) An (the line-of-sight integral through the halo density squared, averaged over 
the angular acceptance AQ) vs. AQ for several halo profiles: smooth isothermal, 
clumpy isothermal, Burkert, NFW, and Moore et al. Taken from j6T] . 



We have yet to discuss the dark matter density profile itself, and since its precise shape is 
not well known, we will survey a representative sample of profiles rather than focussing on any 
single model. In general, these profiles are derived from numerical simulations |101] . and most 
take the form 



p(r) 



Po 



(6.11) 



(r/i?)7(l + (r/i?)")(/5-7)/"' 
where the density function p has been expressed in terms of the radius r away from the galactic 
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center and spherical symmetry is assumed. The three power-law indices a, /3, and 7, along 
with the characteristic radius i?, serve to define a given model. The models we examine here 
include the relatively cuspy Moore et al. profile [102j . the widely used Navarro- Frenk- White 
p.ome . ana a pa. ot .ot.... .ode. „. one a ..ooth ae..t. d..*u«o.. t.e 
other with a greater level of dark matter clumping |j We list the a, j3, 7, and R assignments 
which define each of these models in table 16. 1[ We also include in our survey one model in 
which the dark matter profile is not defined by equation (16. lip , but by the relation 

where ro and po are fiducial distance and density parameters, respecively (for our galaxy, 
ro ~ 11.7 kpc and po ~ 0.34 GeV cm~'^ [l05j ). This model was originally proposed by Burkert 
et al. |106j and is also in use. Figure lOl shows the relationship between {J{ip))An and AQ that 
results from plugging each of the distributions we consider into equation (16.101) . As one would 
expect, the more sharply peaked distributions (the Moore and Burkert profiles), in which the 
dark matter is concentrated near the galactic center, yield a higher value of {J{'ip))An for a given 
AQ, and therefore to a higher 7-ray fiuxj^ It should be noted that different profiles can lead to 
drastically different values for $. The situation is further complicated by the possibility that 
the presence of a massive black hole at the center of the galaxy could significantly modify the 
dark matter abundance near the galactic center and lead to a pronounced density spike |108j . 



^The dark matter density given by the Moore et al. profile is divergent for ?■ <C -R and is assumed to be 
truncated at some small radius r^. Here, we choose rt = 10~^ kpc. 

^It is worth noting that the models surveyed here assume a spherically symmetric halo, and it has been ob- 
served |107j that triaxial models provide a better fit to the cosmological density profiles obtained in cosmological 
simulations than do spherical ones. The corrections introduced by taking such considerations into account would 
amount to an 0(1) factor in the overall result, and since it will soon become evident that detection prospects 
will be an issue of orders of magnitude rather than numerical prefactors, such prefactors do not significantly 
compromise our results. 
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though there is some debate over the precise effect such a black hole would have, and for this 
reason we will not consider such modifications here. 

6.2 Outlook at Present and Planned Facilities 

Now that we have calculated $ (up to considerations involving the shape of the dark matter 
profile), let us turn to a discussion of the performance of present and future 7-ray telescopes. 
These can be divided into two main types: satellite facilities, which detect incoming photons 
directly, and atmospheric Cherenkhov detectors (ACTs), which operate by observing showers 
of Cherenkhov light that result when high-energy 7-rays enter the Earth's upper atmosphere. 
In assessing performance and the ability to register a discovery at the 5a level, the relevant 
detector attributes to consider are angular acceptance AQ, energy resolution AE/E, and ef- 
fective collection area Aeff. In table 16.21 we provide a list of these performance specifics for 
a variety of present and planned experiments, including both satellite facilities and ACTs, as 
well as those for two hypothetical facilities (one a satellite detector, one an ACT) slightly more 
advanced than any currently planned facility of its kind, which we will use when examining 
the potential for observation of heavy Wino or Higgsino dark matter annihilation at the next 
generation of 7-ray experiments. 

The advantage of satellite facilities, examples of which include EGRET and the soon-to- 
be-launched pair-production telescope GLAST jll3j , is that they can have excellent angular 
resolution and energy resolution. The field of view for GLAST, for example, will be on the order 
of a steradian, its angular resolution in the TeV range will be ~ 0.1°, and its energy resolution 
will be on the order of 4%. The background seen by such facilities is — as one might expect — 



104 



Facility 



Aeff (cm^) AE/E Afimin (sr) ehad 



HESS (Namibia) 
VERITAS (Arizona) 
EGRET (Satellite) 
GLAST (Satellite) 



WHIPPLE (Arizona) 
GRANITE II (Arizona) 



3.5 X 10« 30% 1.88 X 10"^ 1.0 

5 X 10« 20% 9.56 X 10'^ 1.0 

7 X 10^ 15% 9.56 X 10-6 0.25 

1 X 10^ 15% 3.83 X 10-^ 0.25 

1 X 10^ 15% 3.22 X 10~2 

1.5 X 10^ 4% 9.56 X 10"6 



Next generation ACT 
Next generation PPT 



1.5 X 10^ 10% 1.00 X 10-^ 0.25 
2 X 10^ 1% 1.00 X 10-^ 



Table 6.2: The performance parameters |109l IllOj for current and planned 7-ray telescopes, 



including both ACTs (WHIPPLE, GRANITE II, HESS [m], and VERITAS [TT2] ) 
and space telescopes (EGRET and GLAST |113j ). Also included are the parameters 
corresponding to the hypothetical "next generation" atmospheric Cerenkhov tele- 
scope (ACT) and space-based pair production telescope (PPT) we have used in our 
analysis. 



the actual diffuse gamma-ray background, which is not currently well known for energies in 
the TeV range. The best that can currently be done is to make the assumption that the 
power law spectrum from EGRET data (good up to ~ 100 GeV) can be extrapolated to higher 
energies |114[ I115i 1116] . In doing this, one obtains a power-law of the form 



with a numerical prefactor A^^o on the order of 10 ^ cm^ and an exponent a somewhere between 
2.0 and 2.5 (following [nS], we take a to be 2.1 and Nq to be 7.32 x 10^). Since GLAST 
will provide a great deal of information about the diffuse 7-ray background at high energies, 
we can expect the uncertainties in a and Nq to be dramatically reduced once it begins taking 
data. In figure [621 we plot the 7-ray flux that a satellite detector would observe as a function of 
angular acceptance along with the 7-ray background flux (for several different values of AE/E) 
determined from equation (16.21) . It should be noted that while the signal contours are largely 




(6.13) 
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independent of misp, the background contours are not: the power-law relation ( ]6.13p indicates 
that with positive a, background flux will decrease with increasing m^^p. Thus the results given 
for misp near the bounds (16. 2p - (16.31) represent the best-case scenario. One can see that for the 
sharply peaked Moore et al. and NFW profiles, the signal exceeds the background significantly 
when AQ is small. 

10-^° 

'e 

LL 

C 

10"'' 10"^ 10"^ lO"* 10"^ 10"^ 10"'' 10° 10^ 

AO. 




Figure 6.2: The expected flux from the annihilation of a heavy Wino LSP, as a function of AQ, 
that would be detected by a satellite detector aimed at the galactic center [61j. Also 
included is the anticipated background flux for a 2.3 TeV Wino at such a detector 
for different values of detector energy resolution, ranging from AE/E = 50% (top 
line) to AE/E = 5% (bottom line). It should be noted that the spread in the signal 
is smaller than the width of the energy bin. From this, it is apparent that for the 
NFW and Moore et al. profiles, the prospects for detection increase with better 
angular resolution (decreasing AQ). For a Higgsino LSP, the resulting curves are 
similar, but the signal is two orders of magnitude lower. The background flux for a 
1.1 TeV Higgsino also is increased by a factor of ~ 10 over the 2.3 TeV Wino case, 
owing to the lower energy of the signal photons. 

Satellite detectors do have a major drawback, however, which is that they are limited by 
collection area constraints. In order for the signal registered at any detector to be interpreted 
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as a discovery, not only must the significance level (the ratio of Ns, the total number of 



signal photons registered, to V^bg, where Nbg is the total number of background photons 
registered) exceed 5a, but the total number of detected photons must exceed 25, the threshold 
below which Poisson statistics give an equivalent confidence limic|. These requirements, when 
written explicitly in terms of A^s, AE/E, AQ, and observation time, are 



^V<l>BG(A^], AE/E, Chad) ' " 

<^{AQ)A^st > 25. (6.15) 

Since the effective collection area is constrained by the size of the telescope itself, present 
satellite facilities tend to have an Aes on the order of 10"^ cm^, and dramatic (i.e. order of 
magnitude) improvements on this in the future are unlikely. The consequences of this are 
shown in figure 16.31 Here, contours of total event count at an advanced, hypothetical space 
telescope (the "Next generation PPT" from table 16.21) are presented alongside significance 
plots corresponding to each of the halo models discussed above for both a Wino and Higgsino 



LSP, assuming an exposure time of 10'' s. While obtaining the necessary Ns / \/ Nbg to claim 
a discovery is difficult in itself (and only achieved for the Moore et al. profile), the most 
troublesome issue is that none of the halo profiles surveyed yield enough events to register a 5a 
discovery. Space telescopes are thus "event-count-limited," and as a result, the prospects for 
detecting monoenergetic photons from heavy LSP decay at such facilities (including GLAST) 
are quite dim. 

For ACTs, the situation is quite different. This class of detector (examples of which include 
the HESS |111] and VERITAS |112] arrays) operates by observing showers of Cherenkhov light 



■^While the likelihood of random statistical fluctuations at the 5a level increases with improved energy 
resolution, these can be differentiated from a true signal by requiring the signal to be consistent over multiple 
trials. 
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Figure 6.3: The ratio of A^signai/ a/ A^background (left panels) and total number of photons (right 
panels) collected by a generic space telescope with an effective area of 2 x 10^ cm^ 
and an energy resolution of 1%, over a range of AQ, after 10^ s (about 1/3 of an 
active year) of viewing time, and for both Wino (top panels) and Higgsino LSP 
(bottom panels) [61]. See figure [6^ caption for the halo model key. The threshold 
for 5(7 discovery has been included for reference in the significance graphs, and 
contours corresponding to 10, 25, and 100 events have been included in the event 
count graphs. It can be seen here that for such a space telescope, no halo profile 
is capable of producing the 25 events necessary for detection, and that only for the 
Moore et al. profile is the significance criterion even achieved. 



produced by the incidence of high-energy 7-rays on the Earth's upper atmosphere. The general 
scale of the effective area associated with an ACT is determined by the area of the Cherenkhov 
light pool on the ground, which is ~ 5 x 10^ cm^. Compared to the aforementioned A^s 
scale associated with satellite facilities, this is enormous. The trade-off is that uncertainties in 
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reconstructing the energy of the primary photon from the properties of the radiation shower 
place hmits on the energy resolution, which is generally quite poor: a single imaging detector 
can achieve AE/E ~ 30 — 40%; an array of parallel detectors, 10 — 15%). Furthermore, 
the background "seen" by an ACT is not simply the diffuse 7-ray background. Any event 
that precipitates a similar Cherenkhov cascade is a source of background, and the primary 
background at such detectors actually comes from cosmic-ray protons, electrons, etc., which 
dominate over the diffuse gamma-ray background by an order or two of magnitude. It is possible 
to discriminate between hadronic showers and those initiated by 7-rays to a degree due to the 
shape of the cascade and to the time spread of the light pulse, but showers initiated by leptons 
(predominately electrons) are indistinguishable from 7-ray cascades. These backgrounds are 
higher than those seen by satellite detectors, though their spectra are reliably known up to 
5 TeV . The power-law behavior |117j for hadronic and leptonic background events is given by 



where we have replaced A'^o and the power-law index a with their explicit numerical values. 
Here, ehad is a dimensionless coefficient which represents the detector's ability to reject hadronic 
events based on the criteria discussed above, normalized to the performance of the WHIPPLE 
telescope (for which ehad = 1)13 The total observed background is the sum of the hadronic and 
leptonic backgrounds given in equations (16.161) and (16.171) . Improvements in hadronic rejection 
techniques have since lowered ehad to around 0.25 at instruments such as HESS and VERITAS. 

Once again, both hardonic and leptonic contributions to the overall Cherenkhov background 

^In other words, Chad = 0.5 indicates a hadronic rejection twice as good as WHIPPLE'S and Chad = indicates 
perfect hadronic rejection and a background consisting purely of leptonic and photonic events. 




(6.16) 



(6.17) 
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scale with photon energy by a power law with a negative index; hence the best-case scenario 
will involve an LSP mass near the upper bounds given in equations (16. 2p - (16. 3p . 




10"'' 10"^ 10"^ lO"* 10"^ 10"^ 10"'' 10° 10^ 
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Figure 6.4: The expected flux from the annihilation of a heavy Wino LSP, as a function of AQ, 
that would be detected by a Cherenkhov detector aimed at the galactic center [6T]. 
Also included is the anticipated background flux from a 2.3 TeV Wino at such a 
detector for different values of detector energy resolution, ranging from AE/E = 
50% (top line) to AE/E = 5% (bottom line). It should be noted that the spread 
in the signal is smaller than the width of the energy bin. From this, it is apparent 
that for the NFW and Moore et al. profiles, the prospects for detection increase 
with better angular resolution (decreasing AQ). For a 1.1 TeV Higgsino LSP, the 
resulting curves are similar, but the signal is two orders of magnitude lower. The 
background flux is also increased by a factor of ~ 10 over the 2.3 TeV Wino case, 
owing to the lower energy of the signal photons 

In figure 16.41 we plot the 7-ray flux observed at an ACT a function of angular acceptance 
for the different halo profiles surveyed along with background flux contours corresponding to 
different values of AE/E. The background is larger than in the space telescope case (compare 
figure 16.21) by about a factor of ten. In figure 16.51 we plot both the signal significance and total 
number of signal photons recorded at an advanced, hypothetical Cherenkhov array (the "Next 
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generation ACT" from table [6l2|) for a range of AQ and an exposure time of 10''' s, in both the 
Wino and Higgsino LSP cases. It is evident from this plot that despite the larger background 
flux and poorer energy resolution, ACTs are far more likely to detect a statistically significant 
7-ray signal from heavy LSP annihilation than are satellite detectors, due to their far greater 
effective area (on the order of 10^ — 10^ cm^). Since HESS and VERITAS have an angular 
resolution (and hence minimum AQ) on the order of 10~^ steradians, detection would not be 
difficult if the dark matter density conformed to one of the more cuspy profiles. However, since 
the field of view (and hence maximum AQ) for an ACT is generally around 10~^ steradians, 
the less sharply peaked dark matter distributions (the Burkert profile and the two isothermal 
models) will still be out of reach. It is an interesting coincidence that the NEW profile nearly 
demarcates the line between detection and non-detection for presently operational facilities: if 
the actual dark matter distribution is cuspier than that given by the NEW profile, the 7-ray 
signature of Wino dark matter in PeV-scale split supersymmetry should be detectable at the 
next generation of ACTs; if the actual profile is much less sharply peaked, it is unlikely that 
such a signal would ever be detectable at an ACT. 

As mentioned previously, the energy resolution for ACTs is limited by considerations related 
to event reconstruction from the Cherenkhov cascade [109] , and hence comparatively poor — 
around 10 — 15%. Due to this intrinsic limitation, it is unlikely that future facilities will offer 
significant (i.e. order of magnitude) improvements in AE/E over present facilities. Since the 
splitting between the 77 and Z7 lines (16. 8 p requires AE/E = 0.8% when tjilsp = 500 GeV, and 
even more precision for LSP masses near the upper bound (16.21) - (16. 3p from thermal generation, 
the prospects for differentiating these lines at an ACT are effectively nil. Still, even if these 
lines cannot be resolved, the detection of any effectively monoenergetic photon signal in the 
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TeV range would serve as provocative evidence for heavy neutralino dark matter. 



Ns/NbI for Wino LSP Ng for Wino LSP 




Figure 6.5: The ratio of A'signai/ \/ -^background (left panels) and total number of photons (right 
panels) collected by a generic Cherenkhov array with an effective area of 1.5 x 
10^ cm^ and an energy resolution of 10%, over a range of Af2, after 10^ s (about 
1/3 of an active year) of viewing time, and for both a heavy Wino (top panels) and 
Higgsino LSP (bottom panels) [61] . See figure 16.41 caption for the halo model key. 
The threshold for 5cr discovery has been included for reference in the significance 
graphs, and contours corresponding to 10, 25, and 100 events have been included 
in the event count graphs. It can be seen here that there are real prospects for 
detection with such a Cherenkhov detector, provided that the galactic CDM halo 
density resembles the NFW or Moore et al. profiles. 

Thus far, we have examined the potential for detecting energetic photons from heavy neu- 



tralino annihilation at current and planned facilities (or hypothetical ones only slightly more 
advanced than these), but it is also interesting to take a slightly different point of view and ask 
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exactly what attributes would a detector need in order to conclusively register a discovery. We 
have seen that the effective area A^jj and angular acceptance AQ (which may be tuned to any 
value between the detector's angular resolution and its field of view) are the most important 
of the detector attributes in the detection of the LSP annihilation signal, so these variables 
will become our primary focus. In figure 16.61 we plot contours from both the significance con- 
straint fl6.14p and the event count constraint (16.151) in A^s-AQ space. Since the significance 
contours are slightly different for space telescopes and ACTs, which "see" different backgrounds, 
we present a separate plot for each case (an ACT in the top left panel; a space telescope in 
the bottom left). The results shown are for 10^ s of viewing time, and the AE/E values used 
in computing the significance limits are those given in table 16.21 for the GLAST telescope and 
the VERITAS array. The horizontal bars correspond to GLAST and VERITAS themselves, 
and the endpoints thereof are given by each detector's respective angular resolution and field 
of view. The criterion for discovery at either of these facilities is that for a given halo model, 
both the 5(7 and A^^ > 25 contours for that model lie below the corresponding bar. From this 
figure, it is once again evident that GLAST, primarily due to its small effective area, would be 
unable to detect CDM from PeV-scale split supersymmetry at all, while VERITAS would have 
far better hopes for detection. 

Perhaps more importantly, however, it is also evident from figure 16.61 what improvements 
in detector area, angular resolution, and field of view would be necessary in order to detect a 
signal in even the least sharply peaked of halo models surveyed. A space telescope would need 
about a factor of 5 increase in A^jj in order to register the requisite number of signal events to 
claim a discovery even in the Moore et al. profile case, and this would require significant feats of 
engineering. The best way to increase the effective area of an ACT is to add further telescopes 



113 

to the detector array, but the increase in A^jj from each such addition is merely additive. We 
see from figure 16.61 that in order to register a discovery when the dark matter profile is one of 
the those included in our survey less sharply peaked than the NFW profile (the Burkert profile 
and the two isothermals) , a factor of ten increase in A^ff is required (one could theoretically 
also widen the field of view by several orders of magnitude, but in practice this is even more 
difficult). This is in principle possible, of course, if one adds enough telescopes to a given array, 
but this is an expensive and somewhat impractical proposition. 

In addition to the discrete lines arising from to NiNi 77 and NiNi 7Z, the 7- 
ray spectrum from heavy neutralino annihilation has a continuum component which arises 
primarily from pion decay. The requirement that this continuum contribution not exceed the 
flux observed by EGRET [IISJ provides another phenomenological check on the model. It 
was shown in |90] that the continuum photon flux produced by the decay of a heavy Wino or 
Higgsino LSP will be low enough that no conflict arises between the predictions of the dark 
matter scenario model and EGRET data. This is one advantage of heavy Wino or Higgsino 
dark matter in AMSB scenarios: this constraint is not a trivial one, and may be relevant in 
many models with a lighter LSP. 

So far we have considered only the contribution to the 7-ray flux arising from annihilations 
in our own galactic halo, but it is also interesting to consider the effects of dark matter annihi- 
lation in other galaxies to the total observed photon background. The 7-ray flux arising from 
extragalactic WIMP annihilation has been investigated by several authors |118[ 1119] 11201 [T2T] . 
The results depend significantly on the choice of halo profile and on other astrophysical inputs 
about which there is substantial uncertainty, including the number density distribution (often 
also called the cosmological mass function) ^(M, z) of halos in the universe as a function of 
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halo mass M and redshift z and the correlation between halo mass and dark matter concen- 
tration. The situation is further complicated by the possibility that a nontrivial fraction of 
dark matter could be bound in small subhalos within larger, virialized halos [122[ I123j . for in 
this case the contribution to the 7-ray flux from extragalactic dark matter annihilation can 
be substantially increased. Owing to these uncertainties, the best way to proceed is to choose 
a reasonable astrophysical model that is not too conservative (in terms of the resulting 7-ray 
flux) and ascertain whether the continuum spectrum from LSP annihilation is consistent with 
EGRET data. Results corresponding to the choice of a Press-Schechter |124] cosmological mass 
function, an NEW halo profile without subhalos, and a Wino or Higgsino LSP with a mass of 
180 GeV were derived in |118j . It was shown that for the relatively cuspy NEW profile, the 
continuum component of this flux was not large enough to conflict with EGRET data, and 
one would expect a heavier LSP of the same sort would also be compatible with the observed 
7-ray spectrum at energies i?^ < 100 GeV. Eor a cuspier distribution, such as the Moore et. al 
profile, or a model where a substantial fraction of the dark matter is concentrated in subhalos, 
however, the continuum flux could potentially come into conflict with these measurements. 

In addition to the continuum flux, the extragalactic photon spectrum from LSP annihilation 
also has a "discrete" component — the NiNi 77 and NiNi 7Z lines, which will be smeared 
out somewhat due to absorption and redshifting effects. If the diffuse 7-ray background in the 
TeV range indeed obeys the same power-law relation fl6.13p that it does at lower energies, 
the 77 and Z7 peaks will likely be difficult to detect in particle physics models of this sort. 
However, if the contribution to the 7-ray spectrum from all other sources drops significantly 
when > 100 GeV, as is the case in certain models where blazar emissions dominate this 
spectrum at high energies |118j . these peaks could dominate over the rest of the background 
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by several orders of magnitude. In practice, however, it is unlikely that the signal from the 
extragalactic annihilation of a TeV-scale Wino or Higgsino LSP offers much in the way of 
additional prospects for indirect detection. It was illustrated in figure [616] that GLAST, like all 
space telescopes, is hindered in its ability to meet the event count requirement (16.151) necessary 
for signal detection by its small effective area. Even when AQ is adjusted to near the telescope's 
field of view, the extragalactic photon flux from LSP annihilation will likely be insufficient to 
meet this requirement, and hence there is good reason to think that GLAST will fare little 
better in detecting extragalactic LSP annihilation than it will in detecting LSP annihilation 
at the center of our own galaxy, no matter what the shape of the diffuse photon background 
is at high energies. AGTs are also unlikely to detect this signal, since the background they 
see (16.161) - (16.171) is an order of magnitude larger than that seen by space telescopes and does 
not depend on the shape of the diffuse 7-ray background. 

The message here is that while the 7-ray spectrum from extragalactic LSP annihilation 
is interesting both because it provides a consistency check on a given particle physics model 
(that the continuum spectrum does not conflict with EGRET data) and because it contains 
distinctive structures (the "monoenergetic" peaks) at energies as yet unprobed, the former 
consideration is not damning for the case of heavy Wino or Higgsino dark matter, and the 
latter is not likely to yield any new possibilities for indirect detection. As with the case with 
the signal from LSP annihilations in our own galactic halo, the difficulty in detecting the signal 
from extragalactic LSP annihilation stems from the limitations of satellite detectors due to 
their small collection areas. 

To summarize the results of this chapter, although energy resolution considerations rule out 
discrimination between the 77 and Zj lines from heavy LSP decay, the detection of a single 
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monoenergetic photon "line" is still a possibility at atmospheric Cherenkhov detectors like HESS 
and VERITAS, which have a large enough effective collection area to register a statistically 
significant discovery. Space telescopes like GLAST, on the other hand, have far smaller effective 
areas and hence will be unlikely ever to observe the signal in question. Detection at Cherenkhov 
facilities is far from guaranteed, however, and prospects for it depend significantly on the 
dark matter distribution in our galaxy. Continuum photons are also produced during LSP 
annihilations, both in our galactic halo and in the halos of other galaxies, but these are not 
likely to run afoul of current experimental constraints on the 7-ray background unless the 
universal dark matter profile is extremely cuspy or involves substantial substructure. 
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Figure 6.6: Detection boundary contours in — AQ parameter space for the 7-ray signature 
of a 2.3 TeV Wino, based on the 5a significance requirement from equation 16.141 
for VERITAS (top left panel), with = 1 x 10^ cm^ and AE/E, and GLAST 
(bottom left panel), with A^s = 1.5 x 10^ cm^; and based on the Ns > 25 event 
count requirement from equation 16.151 (right panel), for a variety of halo profiles. 
Bars showing the range of angular acceptances that can be chosen at VERITAS 
and GLAST have also been included. In order to register a discovery at either of 
these facilities for a given halo model, both the 5a and Ng > 25 contours for that 
model must lie below the bar corresponding to that facility 



CHAPTER VII 



CONCLUSIONS 



Leptogenesis is an attractive model, and the prospect that its abihty to relate the smallness 
of neutrino masses to the mechanism responsible for baryogenesis could be realized in a sce- 
nario where neutrinos are purely Dirac is an interesting one for many reasons. We have shown 
that thermal Dirac leptogenesis is not only an intriguing theoretical curiosity, but also a gen- 
uinely viable phenomenological model. Not only is it capable of producing a baryon-to-photon 
ratio 1] for the universe that matches that observed by WMAP, but it can also satisfy all rele- 
vant constraints from cosmology, flavor physics, etc. while reproducing the observed neutrino 
spectrum. 

Specifically, from the research presented in this work and that originally presented in the 
publications [Ml 1121 EI] and reviewed herein, the following results and conclusions were ob- 
tained: 

• There is a particularly simple, theoretically-motivated realization of Dirac 
leptogenesis, constrained hierarchical Dirac leptogenesis (CHDL), which re- 
produces the observed neutrino spectrum. 
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In its most general form, the Dirac leptogenesis superpotential ( 12. ip contains a large 
number of new free parameters, including the heavy particle masses M$-, the scalar 
VEV X the elements of the general, complex trilinear coupling matrices A and h. 
It is thus not terribly surprising that some configuration of parameters is capable of 
satisfying all experimental bounds on neutrino masses and mixings, but it is of interest 
that this set of constraints can be satisfied in a simple, theoretically-motivated model that 
contains only five free parameters (four, if one does not include the highly constrained 
coupling matrix entry 63). In [38], it was shown that a particularly simple model of Dirac 
leptogenesis, which was dubbed constrained hierarchical Dirac leptogenesis (CHDL), in 
which the trilinear coupling matrices A and h are antisymmetric, the mass matrix M$ 
is real and diagonal, and the small ratio 5 = M$^/M$^ dictates the textures of the 
neutrino mass matrix, the experimental constraints are readily satisfied. The model is 
simple, robust and well- motivated by flavon physics (see section 13721) . Furthermore, as was 
shown in [3H] and reviewed in chapter 14. 2^ this model is compatible with the additional 
constraints arising from baryogenesis, astrophysics, and cosmology. 

• Dirac leptogenesis strongly prefers a heavy gravitino with 7713/2 ^ 10^ GeV and 
is compatible with split supersymmetry. 

Dirac leptogenesis must respect a battery of astrophysical constraints in order to be con- 
sidered a phenomenologically viable theory, and the most stringent of these are related to 
gravitino cosmology. The connection between this subject and leptogenesis occurs through 
the reheating temperature Tr. On the one hand, Tr < is required for thermal lepto- 
genesis; on the other, for certain values of the gravitino mass m3/2, considerations related 
to BBN and nonthermal LSP production place stringent limits on Tr. In [38], it was 
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shown that baryogenesis considerations require > 10^° GeV, unless some additional 
mechanism — such as resonant leptogenesis — is invoked. Thus, when 7713/2 ^ 10^ GeV, 
nucleosynthesis bounds on Tr from late gravitino decay render Dirac leptogenesis in its 
simplest form essentially unworkable; when 10^ GeV < ^,3/2 ^ 10* GeV things are better, 
but constraints related to nonthermal LSP creation also cause problems when the split- 
ting between m3/2 and misp is small. These results are summarized in table mi For this 
reason, Dirac leptogenesis strongly prefers a heayy gravitino with m3/2 ^ 10^ GeV. Since 
this is a natural consequence of split supersymmetry, this scenario becomes one possible 
context for Dirac leptogenesis. This scenario has the added advantage of keeping squark 
and slepton masses heavy to alleviate flavor-violation concerns, while keeping gauginos 
light to provide a (characteristic) dark matter candidate. 

• Dirac leptogenesis can succeed in models with light (m^ = 200 GeV) scalar 
masses. 

In any model where supersymmetric scalar particles are light and where the mass eigen- 
states of these particles are not diagonal in the flavor basis of the Standard Model quarks 
and leptons, loop diagrams involving such particles can potentially give rise to unac- 
ceptable levels of flavor violation. In Dirac leptogenesis, the renormalization group run- 
ning of slepton soft masses will tend to generate off-diagonal terms in the sneutrino and 
charged slepton mass matrices. It was shown in [13j that flavor violation constraints do 
not rule out CHDL: the model is successful for a universal scalar mass as low as 
200 GeV, respecting all lepton-flavor- violation constraints from processes such as /i ^ 67 
and r — > yU7. Dirac leptogenesis therefore does not require split supersymmetry and will 
be phenomenologically viable in a wide range of contexts. More general AMSB [581 El] 
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models, for example, can yield a splitting between and misp large enough to permit 
a gravitino with a mass slightly above the BBN bound {171^/2 < 10^ GeV) and sparticle 
masses of (9(200 TeV) without running afoul of any of the late gravitino decay constraints. 

• Obtaining the appropriate superpotential for Dirac leptogenesis can be ar- 
ranged in a consistent, anomaly-free manner through the introduction of a 
hidden sector U{1)n symmetry. 

The Dirac leptogenesis superpotential requires some additional symmetry to forbid tree- 
level Dirac and Majorana masses for right-handed neutrinos. It was shown in section 15.31 
that this can be accomplished in a simple, consistent, anomaly-free manner by enlarging 
the hidden sector of the model. It was shown that the necessary VEVs are obtained to 
yield small but nonzero masses and to ensure that all i?-parity odd fields in the hidden 
sector become massive after supersymmetry is broken. 

• Right-handed sneutrino dark matter is a viable possibility in Dirac leptogen- 
esis. 

Unlike in Majorana leptogenesis, in Dirac leptogenesis the lightest right-handed sneutrino 
can be light, and as it is neutral under the Standard Model gauge group, it is a 
potential dark matter candidate if it is the LSP. Since the interactions of this particle are 
extremely weak (as they are suppressed by the powers of the effective neutrino Yukawa 
coupling), its relic abundance will in general be far larger than the upper bound fl3.4p 
on VLcdm obtained from WMAP data. However, Dirac leptogenesis comes equipped with 
a mechanism for alleviating this problem: the theory requires an additional symmetry 
(to forbid Majorana masses for neutrinos) under which z/^^ is necessarily charged. In 
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section [573l it was shown that if this symmetry is a gauged U{1), annihilations through 
the gauge bosons and gauginos associated with this additional symmetry can reduce the 
right-handed sneutrino abundance to an appropriate level. 

• Linking ^Icdm to Vlb via an asymmetry in some conserved quantum number is 
generally quite difficult, at least in supersymmetric theories, and requires a 
hierarchy among annihilation rates. 

The possibility that the cold dark matter abundance flcDM and the relic density Qb of 
baryonic matter in the universe, which are of roughly the same order, might be linked is an 
intriguing one. One way of arranging this is the method of [93], in which the association 
is forged by charging a hidden-sector dark matter candidate under a globally conserved 
quantum number. Then dark matter particles annihilate with their antiparticles until 
only the excess given by the asymmetry is left, and this excess becomes the dark matter 
abundance. This cannot work in the usual neutralino dark matter scenario, since Majo- 
rana neutralinos cannot develop a charge asymmetry, but a right-handed sneutrino LSP in 
Dirac leptogenesis, which carries lepton number charge, has the right properties to allow 
this mechanism to work. Unfortunately, linking the right-handed sneutrino abundance to 
the lepton number L^,^ stored in z//j^ in this manner proves difficult in practice. It was 
shown in section 15.31 that this required the Lj^^j-violating rate Fp,) to be small, so as to 
preserve the lepton number asymmetry at the appropriate level, whereas it required the 
Lj^^-conserving rate Tj^c* to be several orders of magnitude larger, so that right-handed 
sneutrinos could annihilate efficiently and dark matter would not consequently be over- 
produced. However, in a supersymmetric model, s-channel diagrams in which a dark 
matter particle and its antiparticle annihilate through some virtual Z' are accompanied 



123 

by t-channel diagrams (with Majorana gaugino intermediaries) that violate L^^. It is 
worth restating that the evolution of Lj,^ during the CDM freeze-out epoch will not affect 
the baryon number B of the universe whatsoever (the effect of Tj^c is simply to shuffle lep- 
ton number from uji-^ to and other light hidden sector fields, and not to change Lhid), 
and so right-handed sneutrino dark matter is perfectly compatible with baryogenesis in 
Dirac leptogenesis scenarios. However, it was shown in section [5^ that while it is possible 
to obtain the correct ^Icdm with a right-handed sneutrino, the inherent tension between 
Tio) and T^^ makes it very difficult simultaneously to prevent L^^^ from being washed out 
to the extent where it will have no bearing on the dark matter abundance. It was shown 
that there are ways of overcoming this hurdle, such as adding a large number of light 
scalars to the theory, but these methods are somewhat contrived. The prospects may 
be brighter for linking flcDM and fib in other situations where a dark matter candidate 
is charged under some conserved quantum number, but since it is the supersymmetriza- 
tion of annihilation processes that lies at the root of the problem, the issue is a general 
one — especially since supersymmetric models naturally provide a compelling solution to 
the dark matter problem, in the form of the LSP, when i?-parity is conserved. 

• The prospects for detecting the monoenergetic photon signal of heavy Wino 
or Higgsino dark matter reasonably good at atmospheric Chernenkhov tele- 
scopes (ACTs) like HESS and VERITAS, but slim at space-based facilities 
like GLAST. 

Since split supersymmetry emerges as one promising context for CHDL, it is relevant to 
the study of Dirac leptogenesis, and of more general interest as well, to examine the dis- 
covery potential for the dark matter candidate that naturally emerges from this scenario: 
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a predominately Wino or Higgsino LSP with a mass around 1 TeV. The direct detection 
prospects for such a particle are somewhat dim, and the easiest way of discovering such 
a particle is through the observation of the effectively monoenergetic photon signal pro- 
duced by its annihilations at the galactic center. In [6Tj, we show that space telescopes 
such as GLAST are highly unlikely to record the necessary number of signal events to 
claim a statistically significant discovery, due to their small effective collection area Aejj. 
Ground-based Cherenkhov arrays, however, have far larger effective areas and hence far 
better prospects for the discovery of a heavy Wino or Higgsino LSP. 



Gravitino 
Mass Range 
(GeV) 


Maximum 
(GeV) 


Workability 
of CHDL? 


Comments 


"^3/2 < 10^ 


lO" - 10* 


Very Low 


G decay during or after BBN. 
Insufficient rj generated. 


10^ < m3/2 < 10* 


109 - ioi° 

or higher 


Depends on 
'^LSP and the 
ratio 
^LSP /i^^/2 


LSP annihilations ineffective unless 

LSP 1 1 1 small. Tn constrained 
by nonthermal LSP abundance from G 
decay. Loop-split SUSY works only for 
'"^3/2 10 GeV. More general split 
SUSY theories can be successful. 


10* < "13/2 S 10^° 


None 


Excellent 


f2^5P is thermal, since G decays 
before LSP freeze-out. 

hli^^ > 10"'""'" GeV allowed, v spectrum 
requirements compatible with CHDL. 


IQl" < m3/2 < 10^^ 


None 


Questionable 
{depends on 

gluino 
properties) 


V sector and baryogenesis okay, but 
model may have a cosmological gluino 
problem. 


10^^ < "13/2 


None 


Very Low 


Potential gluino problem becomes a 
serious concern. 



Table 7.1: The various gravitino mass regimes for split supersymmetry models and the viability 
of thermal Gonstrained Hierarchical Dirac Leptogenesis (GHDL) in each case. 
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It should be mentioned that there are a number of experimental checks on the viability of 
Dirac leptogenesis in general and on the specific model we have dubbed CHDL. The major 
prediction of Dirac leptogenesis is that neutrinoless double-beta decay will not be observed to 
any degree, for this process relies on the existence of a Majorana mass term for right-handed 
neutrinos. The discovery of such a process experimentally would rule the theory out. In 
addition, forthcoming results from MiniBooNE [125] should either confirm or deny the LSND 
result, which will reveal whether or not the neutrino spectrum produced by CHDL is in fact the 
one present in nature (if MiniBooNE confirms the LSND result, Dirac leptogenesis should still 
be workable, albeit in a more complicated manifestation than CHDL). Finally, CHDL places 
constraints on the neutrino mixing parameter sin ^13 (see figure [231), the value of which will be 
measured in future experiments. 
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APPENDIX A 

Effective Couplings for Flavor- Violating Calculations 



For completeness, we list here the results used in our analysis for lepton flavor violating 
processes. The amplitudes A'" and in equation fl3.44p were computed in ^83j and, with a 
trivial extension to include three right-handed sneutrinos, are given by 



= + A^")'^ and A^ = A^"^^^ + A^")'^, 



where the individual amplitudes A^"^'^, A^")'^, A^"^'^, and A^")'^ 
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Here, the indices A and X respectively label the gaugino (chargino or neutralino) and slepton 
(sneutrino or charged slepton) mass eigenstates, xax = '"^xa/'^^x' ^tix i^iAx) denote 
the effective couplings of charged lepton i to chargino (neutralino) A and sneutrino (charged 
slepton) X. The flavor mixing terms in (I3.38P - fl3.40p enter into the overall rate (I3.44p through 
Cfj^ and A^j^X' which contain elements of the matrices and Ui that diagonalize the mass- 
squared matrices for sneutrinos and charged sleptons, respectively. The slepton masses also 
enter into the partial amplitudes (IA.2tiA.4|) . 
The effective couplings nI'J^ and Cl'J^ are 
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in terms of the chargino mixing matrices {OR)A,i and {OL)A,i the neutralino mixing matrix Uxi, 
and the sneutrino and charged slepton mixing matrices and Uxi- The chargino mixings 
matrices are defined by the relation 
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and M^^"'^ is diagonal. The sneutrino mixing matrix Ux i and the charged slepton mixing matrix 
Ux i are defined by the relations 
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where the matrices mj^ and are given by the sum of the MSSM contribution and the 
respective Dirac leptogenesis contributions in (13.411) . The neutrahno mixing matrix Un is 
defined by the relation 

(m^)'*'^^ = U^mj^Ul, (A.9) 



where 
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APPENDIX B 



Derivation of the Boltzmann Equations for Dirac 

Leptogenesis 



In this appendix, we derive the Boltzmann equations (14 ■2p - (14. 3p for the evolution of baryon 
and lepton number in Dirac leptogenesis, following the methods and notation of |126j . We 
begin by observing that the Boltzmann equation for any particle species a in the early universe 
can be written in terms of the number density of a as 



dn 

— +3Hna 
at a 



d^Pa d^pi d^pj d^pk ,^ .4.. . 
(27r)3 (27r)3 (27r)3 ' ' ' [27,)^^ ' ^ _^ 



n=a,i,j,...k 



J2\M{a...z ^ j...k)\\U..f,) 

int. 



int. 



(B.l) 



where is the phase-space distribution function of particle i, \ J^{a...i are the squared 

matrix elements for particle-number-changing interactions involving a, and the sums are over all 
interaction processes which create or destroy a. The term proportional to the Hubble parameter 
H can be absorbed by rewriting equation ( IB. II) in terms of the ratio Ya = Ua/s, where s is the 
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entropy density of the unverse [42] : 

(IYa 1 f d^Pa d^Pi d^pj d^pk 
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In the high temperature hmit T ^ m^, the ratio Ya is well- approximated by 



Ya = Ua/s ~ - 



1 Qa 



[B.3) 



where Qa and /i^ respectively denote the the number of degrees of freedom and chemical potential 



of species a and g^^ 
at temperature T|l 



ST) represents the total number of relativistic interacting degrees of freedom 
If an asymmetry Aa = Ya — Yac develops between a and its antiparticle a'^, 
the chemical potential /i^ can be expressed in terms of this asymmetry (assuming /x^, ^ 1) 
using the relation 

.2Ma/r_i)^ (B.4) 
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which implies that 



9*s 

9a 



(B.5) 



In deriving Boltzmann equations for the heavy fields, we will begin with those for the scalars 
(f) and (f)'^. To first order in these fields, the leading processes include only decays and inverse 
decays, the Boltzmann equations for the abundances Y^ and Y-r are given by 



dt 
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c zjc\ |2 



(B.6) 



^During the epoch of interest, ff*s(T) remains constant, and thus wo will henceforth drop the T-dependence 
in our notation. 
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and 



^^^^ _ 1.^* 
dt ~ s^'' 



-Ulh \M(J''<Pl - - UJh^ \M{uaHu - , (B.7) 

where A"; ;^ has been used as a shorthand to denote the appropriate phase-space integral and we 
use the notation £ and to denote the usual Lepton and Higgs-up doublets of SU{2). These 
equations can be simplified by noting that energy conservation implies that the inverse decay 
processes of the form {a + i ^ h) appearing in equations flB.6p and flB.7p . in which particle i 
is in chemical equilibrium with other particles in the thermal bath (i.e. /Xj = 0), the previous 
relation allows us to write 

/a/.^/r(l + - — ), (B.8) 

where /^"^ is the equilibrium distribution of a. The assumption of CPT-invariance also simplifies 
these expressions by enforcing the relation 

\M{a = \M{ff a'^)p. (B.9) 

The particle asymmetries with which we will be concerned here are the lepton number 
abundances = Ql- [rii — n1) / s carried by each field i in the theory with lepton number charge 
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Qli- From the charge assignments given in table [2?T] the nonzero Li are given by 
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Using these definitions and relations (IB.Sp and flB.9p . the Boltzmann equations for (p and 
become 



dY^ 

dt 

dY^. 


dt 



(27r)3 
d^p^ 
(27r)3 
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where the interaction rates F^^ = r(0 ^ £ + x) and F^ = F( 
relations 
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z/^ + H^) are defined by the 
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(27r)3 (27r)3' 

with F^, R being the rates of the conjugate processes and F^) = Fj^ + Fr is the total decay rate 
for 0, etc. given in (12. 6p . Because of supersymmetry, the total decay rate of the fermion 
components of the heavy supermultiplets $ and $ will also be F^), with the same F^ and Vr 
as their partial rates. 



134 



It will be more convenient for our purposes to express the Boltzmann equations in terms of 
F^c and . Subtracting ( ]B.17p from ( IB. 161) yields 



dt 



d^P4 



(B.20) 
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-2^^*:(y(r^ + ri) + L„„(r„ + ra)' 

After integrating equations (1B.17|) and (]B.20p over the incoming momentum p^, using rela- 
tions (12. 2p - (12.51) to express the result in terms of the decay asymmetry e, and averaging over 



time-dilation factors |126j we obtain 
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were, the time-dilation-averaged rates (Fj) are given in terms of the naive rates Fj factors |126j : 

(B.23) 



(r.) . r. 



K2{M^/T) 

where Ki{x) and K2{x) are modified Bessel functions. 

The Boltzmann equations for the scalar component of the $ superfield (and its conjugate), 
as well as the ones for the fermion components of $ and are obtained in a similar manner. 



They turn out to be 

dLj,_ 



dt 



dY^_ 
dt 



* 2n7 



eq 



■(F1) + L,,(F^) 



(B.24) 



(B.25) 
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where we have used the fact that Y^'^ = Ff^. Dropping neghgibly small terms proportional 
to eLj, we notice that Boltzmann equations for the combinations Yt == Y"^^ + y^c and 



are redundant. We thus only require three equations to describe the dynamics of the heavy 
field sector. 

We now turn to address the evolution of the Lepton number abundance of the particle 
species i, z//j, i, and i/pt, in which we must take into account the effect of 2 ^ 2 processes 
which transfer lepton number between Le, L^^, Lj, and Lj;^. For the moment, we will not 
concern ourselves with the exact form these interaction terms will take, but will will make one 
important observation: the rates for interactions which shuffle lepton number between i, i, and 
the right-handed charged lepton and sleptons fields and will be much larger than those for 
the interactions which shuffle lepton number between ur and any of these other fields. This is 
because the pr interact only via processes pictured in figure I^TTl which involve a virtual (j)^, 
etc. while all of the other fields either take part in SU (2) and/or f/(l)y gauge interactions. We 
will represent the effects of these rapid equilibration processes by including terms (where 
A is the relevant particle asymmetry) to represent them in the Boltzmann equations. The 
slower 2^2 processes through which right handed neutrinos ur interact with lepton doublets 
i and i and with right handed sneutrino will be included separately as C^^^^, C^^^j and 
Cp^^p^. As discussed in chapter HVl and elsewhere, vr may or may not also take part in rapid 




Y^_ + y^c 




(B.26) 



(B.27) 



(B.28) 
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equilibration processes, depending on the value of (i^x)- 

To simplify further the notation we will define the terms Fa to account for the collective 
contribution from decays (and inverse decays) of the fermionic components of $ and $ (which we 
will not write explicitly, being of similar form to the contribution from the scalar components). 

Expressed in the above notation, the equations for the left-handed lepton field I and its 
conjugate l'^ are 



m 

dt 

dYec 
dt 



1 

— A,* 



'■12 



\M{4>1 



+F^(V':)+s^. 
1 



-A * 



+F£(V^J + Ei 



(B.29) 



(B.30) 



The lepton number abundance is obtained by subtracting the second of these equations from 
the first, and the result is 

dLf 



dt 



+ + L^U] + - + (Sl - Si) 



(B.31) 
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The equations for L^^, Ljr^, and Lj are determined in a similar manner: 



dL 



dt 



(B.32) 



dLj 



-e (r^)(i;. + y;i) + L,^(r,) - y;i^Lj(^-{v,) + le(r^)) 

e {rn){Y,._ + Yp - L, (Tn) - Y;'_^L,, ({Tn) + ^eF^) 



dt 



(B.33) 





(B.34) 



Before going any further we still need to compute the terms C^^^j, C,^^<_>jr^ and C^^^^ 
corresponding to the 2 2 processes mediated by heavy fields. We begin by calculating 
Ci.„^L> which is given by 



Ai' here refers to the amplitude for the specified 2 <-> 2 process to which we have substracted 
the contribution from the resonant intermediate state (RIS) in which a real field 0^ is produced 
and then decayed into the 2 particle final state. The RIS contribution must be substracted since 
we have already counted contributions from decays of real 0^ fields. 

The leading term in the difference between a 2 2 process of the form ah ij involving 
a heavy intermediary k and its conjugate process depends on the contribution of the on-shell 
(resonant) intermediate state: 



M\ab ^ - \M\ij ^ ab)\^ = \MnisitJ ^ ab)\^ - \MTZjsiab ^ tj)\' (B.36) 





{\M'{Ku% ^ ex)\' + \M'{ex - 



(B.35) 



2 
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with 

\Mnxs{ab - ^^^(s - mj)|^(a6 ^ k)\' x \M{k ^ (B.37) 

where k represents the intermediate-state particle, and s is the usual kinematic variable s = 
+ P2^Y- In our case, making use of the equality T^Tr — F^^r^ = eT^ we find 



iMnxsiix - '^^iH^Jl' - \Mnxs{H:v% ^ Ix)? ^ e— — 5(s - mJ)|A^t/ (B.38) 



and substituting this result in eq. ( 1B.35P we finally obtain 



C^^^j = 2eY;i{Tn) + (L^^-hj) n,{va^^^j+vaj^J. (B.39) 

For the other sets of 2 ^ 2 processes, Cy^^^-p^ and Cy^^^, the procedure is essentially the same. 
The rates (F^^^j) = n^{va^^^j+ vai^^^), {Ty^^t) = n^{vayj^^e. + vai^y^^) and (F,,^^^;^) = 
n^ivcTy^^y^ + vay^^y^), assoclatcd with these interactions are calculated in section S21 We 
will denote the total contribution from these processes as (F2«->2)- 

At this point, the full Boltzmann system comprises fifteen individual differential equations: 
four to represent the evolution of (p^, 0_, and their conjugates (in whatever basis we choose); an 
additional four for the fermionic superpartners in the $ and $ supermultiplets; six for the the 
Lepton asymmetries L^, Lj,^, Lj, Ly^, = e^ — e^j, and Lg^ stored in various individual lepton 
and slepton species; and one for the overall baryon number B of the universe, which interacts 
with Li via sphaleron processes of the form given in equation 11.171^ 



We have already noted 



that equations flB.22p and flB.25p are redundant up to terms of 0{eLi), and the assumption of 
unbroken supersymmetry during the leptogenesis epoch allows us to equate the abundances of 
the fermionic and bosonic components of $i and $i, which reduces the number of equations in 



^Keeping track of individual baryon numbers is unnecessary, since QCD gauge interactions among the quarks 
are rapid. 
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our system to ten. As discussed in section 14.11 we can also make use of the fact that many of 
the hght fields in our theory will be brought into chemical equilibrium by rapid SU{2) x U{1)y 
gauge interactions, Yukawa interactions, A-terms, etc. whose collective rates we have denoted 
T,A- In this case, any lepton number stored in i, i, cr, or should be rapidly distributed 
among all of these particles in proportion to the relative number of degrees of freedom for each 
field. Whether is also in equilibrium with these fields depends on the size of the effective 
A-teim (12. lip induced by (Fy.): if (F^) is large, the upt fields will equilibrate chemically with 
the left-handed leptons and sleptons; if (F^) is small (or zero), they will not, but can be seen 
as part of the hidden sector which contains the right-handed neutrino fields. In either case, we 
can define two aggregate lepton numbers Lyis and L^id for the fields in the visible and hidden 
sectors respectively and reduce the number of equations in our Boltzmann system to six. In 
the case where (F^) is large and is in the visible sector, we have 

Lvis = Li + Lyj^ + Li + + Lg^ (B.40) 
Lhid = Ly^, (B-41) 

When vr is in the hidden sector, we have 

Lvis = Li + Ly^ + Lf + Lg^ (B.42) 
Lhid = Ly^ + Le^. (B.43) 

In either case, chemical equilibrium enforces that the lepton number stored in each individual 
particle species be distributed evenly among them in proportion to the number of degrees of 
freedom of each, or in other words 

7 7 

-Le, = -Li = ILy^ = ILe^ = 7Len = Lyis (B.44) 
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in the large- (F-^) scenario, and 

3Le = 3L^ = 6Le^ = Lg^ = Lyis (B.45) 

in the small- (F^) scenario. In the latter case, if chemical equilibrium is established among 
the fields in the hidden sector (for example, via the gauge interactions of the gauged f/(l)y 
discussed in section [5^ . one obtains the similar relation 

= 2^eR = Lhid. (B.46) 

In the large-(F-^) scenario, the Boltzmann equation that describes the evolution of Lms, 
which is obtained by summing equations flB.3ip - flB.34p with the appropriate numerical pref- 
actor from (]B.44p . is 



dLyig 
dt 



-y;j^L„,((r,5) + (r^)) + (l,,, - (r^^^), (bat) 

where small terms proportional to e times Li, L^,^, Lj or Lj>^ have been dropped. Here, we 
have used the fact that the Boltzmann equations for L^.^ and Lg^ are trivial, consisting of only 

terms, which all cancel after taking the sum of all Lepton abundances. 

In the event that annihilation processes of the form 0^0^ X, where X represents some 
final state comprising light fields, are effective in reducing the abundance of (j)^ and its conjugate, 
equation flB.22p must be modified slightly to account for the presence of source and sink terms 



dYA,c I 



-A 



34 



dt s' 



uj^.^ i-M(0,0: - - hf, \M{tj ^ . (b.48 



Here, i and j represent the unspecified final states products of 0^ — annihilation (in which 
2^2 processes will again dominate) and A'^g is the appropriate four-particle phase-space 
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integral. Ignoring negligible CP-violation effects and invoking conservation of momentum, we 
obtain the relation fifj = (/^^)^, which leads to the result 



dt 



H<P,<Pl X)\v\) Yl.^-{Y;if . (B.49) 



Using equation (I4.18p . this can be written in terms of the effective annihilation rate Ta- 

The only task left is to couple the Bolztmann equations for B and L^g via electroweak 
sphaleron interactions. The high-temperature rate for these interactions, F^p/i, is given in equa- 
tion fll.20p . and the proportionality constant between B and L^s is given, in the supersymmetric 
case, by equation (11.351) . and since B is otherwise conserved, the equation for its evolution is 



dB 
'dt 



-{T,,,)(^B + ^L„s) (B.50) 



We are now ready to write down the complete Boltzmann system for the evolution of baryon 
and lepton number in Dirac leptogenesis. Collecting equations (IB. 211) - (]B.22p . (IB. 241) . (lB.32p . 
(]B.47p . and (1B.50P and defining the shorthand expressions 

en en en 

(rD),, = ^^{^d) (Tl),, = ^^{^l) (Tn),, = -^{Tn) (B.51) 
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for the inverse decay rates, we obtain the following result: 



dB 

HI 

dLyiS 

dt 



—'2Lms{(^ d) lo + (^L)ir>^ + {Lhid — -jLns)(^2^2) 

/ 8 



(B.52) 



dL 



hid 



dt 



dYfkc 



dL^ 

dt 

dLs. 
_j 

dt 



2e{Tr>){Y^.^ - Y;1) + L^^ (r^) - 21^^ n) lu 

— {Lhid — jLyis){T2^2) 

-(rD)(Y^l - Y^^) + -Lns{^L)iD + ij^Lhid {^r)i 



- 1 



~(^d)L^^^ + 2Lms(r l) jD + 2Lhid(^R)i 

— (^D)Ltj)_ + 2Lyis{T£))j^, 



(B.53) 
(B.54) 

(B.55) 
(B.56) 
(B.57) 



where once again, negligibly small terms proportional to eLyis or eLhid have been dropped. 
Changing variables from t to 2; = M<^JT, one obtains equations fl4.2p - (14. 7p . 

In the "drift and decay" limit, in which all the rates in (F^) and (r2^2) are assumed to be 
much smaller than the rate of expansion of the universe H, we note that this system simplifies 
considerably. In this case 



dL, 



dt 
dLhid 

dt 
dY2 



dt 



-2e{Vi,){Y^.^-Yl 

2e(r,)(F,. -r-) 

-{Tn){Y,.^-Y;i) 



(B.58) 
(B.59) 

(B.60) 



and since the right sides of these expressions are proportional to one another, the Boltzmann 
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system reduced to a single differential equation, whose solution is 

Lhid(t ^ oo) = -L^s(t ^ oo) = 2eyf^(i = to), (B.61) 

given the boundary conditions Lyis{t = 0) = Lhid{t = 0) = and Y^'^it oc) = and 
approximating F^c by its equilibrium abundance at time to, defined as the time at which 
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APPENDIX C 



Derivation of the Boltzmann Equations for Hidden 

Sector Dark Matter 



Here we present a derivation of the Boltzmann equations for right-handed sneutrino dark 
matter, using the methods of |126j . We begin by writing down the equations in terms of the 
abundances: Yi,^^ = riy^/s and Y^c^ = riyc^/ s: 
dY- 1 

+UJua mu%un - i>Ri>%)\' - U^J,, \M{DIDr (C.l) 

dY 1 

+UJu. mu'^^un - i>Ri>'k)\' - U^J.a \M{i>^^VR ^ u^^un)?] ■ (C.2) 



Here the phase space integral Afj is defined as in equation flB.6p . The energy-conservation 
condition (IB.Sp implies that 

UJuu = e-(^-«i+^-«.)/^(l + ^..L.J (C.3) 

U^J.^,^ = e-(^-«i+^-«.)/^(l-^7*3L.J (C.4) 

UJ.^„ = e-(^-Hi+^--)/^ (C.5) 
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and we also will neglect CP violation (a good assumption at scales T <^ M$^, and define 



\M{v'kVR^T>'kv'k)\^ = \M{vrVr^VrVr)\^ = \Mrr\^ 
\M{vrv%^vrv'^^)\^ = \Mrrc\\ 



(C.6) 
(C.7) 



With these simplifications, the Boltzmann equations become 



dt " 



2 \Mrr\'' [l + g.sL 



n 



dt 



k = 1a34 e-(^^iv«,+i?iv«J/r 



+ \Mrrc\ I i - -mbZmb 
2 \Mrr\' 1 - - 



+ \MrRc\ I i - ^MB^MB 




(C.8) 



(C.9) 



It will be more convenient for our purposes (to relate the dark matter abundance to the hidden- 
sector lepton asymmetry Lhin) to write these equations in terms of Lj>^ = Fjj^ — Y^c^ and 
Ydm — Yi^j^ + Yf^c We also define the effective rates 



r~~ — 



^i'rVr 



1 A 34 p-{ENJ,,+ENJ,^/T\ KA |2 _ ^^fl^fl 



(C.IO) 
(C.ll) 
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where Y^^ = Yj^^{fi = 0), in terms of which the Boltzmann equations become 



dL 



dt 



dV, 



DM 



dt 



2 



^R 



n 



fR 



t^R 



s 



\M 



RR\ 



n 



MB 
i>R 



^R 



n 



MB 



+ \Mrrc\^ I 1 



^UR^i^'k 

^MB^MB 



(T- -* + 2r- - ) i r^"^ - "^^^^ 1 + (T- -2T- - ) 



^•>R 



(C.12) 



(C.13) 



Defining the quantities T^"^^ = Tjj^y^ ± 2ry^y^ and y^'j^^- = 21^^^, these two equations take the 
far more manageable form 



dL 



^R 



dt 
dVoM 



dt 



i'R 



1 

2l 



f(+) f 2 _ y2 

' '-DM 'dm 



(C.14) 
(C.15) 



given in ([ESi-ESTD- 
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APPENDIX D 
Thermally Averaged Cross Sections 

Here we provide a computation of (3"(z>z>) and o"(z>z/*). Let us begin with For a right- 

handed sneutrino LSP, the only relevant process is annihilation into neutrinos via t-channel 
exchange of hidden sector neutralinos x'j shown in figure ID.li It it the Z' component of each 
neutralino that mediates this exchange, hence neutralino mass mixing will introduce factors 
of Uij, where Uij is unitary matrix that diagonalizes (15.191) . The squared amplitude for this 
process is 




Figure D.l: Processes contributing to (3"(i/z/) in the right-handed sneutrino dark matter scenario. 
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in terms of the Mandelstam variables s, t, and u; the U{1)n gauge couphng g and the neutrahno 
masses m^^. Let us assume that m^j -C m^. for alH > 1. Then if we define the the quantities 



^xi ~ ^/™xi' ~ ^/^Ir^ ^ — i'^^n^ " 2a^x^ "~ 1)) ^'^is expression simphfies to 



^ t — m^^ M — Xi"^ 



64g 



^2 _y2 



Using equation (I4.17P and defining r = \/l — Ami /s, "^^ fi^'i 



cr(s 



In 



IT \ A 



4- 



where the s-dependence comes in through both A and r. 



(D.2) 



(D.3) 




Figure D.2: Processes contributing to (t(z/z/*) in the right-handed sneutrino dark matter sce- 
nario. 



149 

The reduced cross-section a^uu*) for L,>-conserving processes is slightly more complicated, 
as we have not only a t-channel contribution from neutralino exchange, but also an s-channel 
contribution mediated by the Z' gauge boson of U{l)j\f, as shown in figure [Dl2l The squared 
amplitude for each of these processes is 

s,s' ^ 

s,s' ^ X ^ 

- V^3J^. (D.5) 

and the interference term between them is 

s,s' ^ 

= —. (D.6) 

^ s- A + Y ^ ^ 

Straightforward calculation yields the reduced cross-section contributions 



= ^ -*'- + A\n(P^'] 1 (D,8) 



the sum of which appears in equations (15.291) - (15.30 
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